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Executive  summary 

During  the  first  year  contract,  we  have  successfully  initiated  and  continued  our 
study  of  defects  in  III-V  and  ZnO  semiconductors  using  ultrasoft  pseudopotential  (USPP) 
method  as  planned.  We  have  also  investigated  and  implemented  a  numerical  approach  to 
calculate  isotropic  hyperfine-parameters  for  isolated  defects  based  on  a  supercell  method. 
We  have  surpassed  our  original  milestones  in  several  aspects  and  already  accomplished 
important  calculations  that  aid  ZnO  and  III-V  research,  a)  Based  on  our  calculations  and 
experimental  work  by  Dr.  D.C.  Look  at  AFRL/WPAFB,  we  together  have  identified  a 
newly  observed  native  defect  donor  in  ZnO  as  a  Zni-N0  defect  complex,  b)  Through  our 
collaboration  with  experimentalists  at  Linkopin  University  (Sweden),  we  identified  the 
recently  observed  EPR  signals  in  diluted  GaPN  to  be  Gallium  interstitial  (Gai) 
complexes,  c)  We  have  investigated  and  identified  two  important  mechanisms  that 
deactivate  nitrogen  doped  ZnO  and  turns  them  into  nitrogen-hydrogen  complex  (No-H) 
and  substitutional  diatomic  molecules  (N2)o-  Our  predicted  infrared  signatures  of  the 
complexes  are  in  agreement  with  the  infrared  measurements  by  our  collaborators  at  the 
National  Renewable  Energy  Laboratory,  d)  We  proposed  a  novel  N-2H  complex  model  in 
GaAsN  that  has  the  IR  signature  matching  the  observed  result. 

Publications: 

1)  Appl.  Phys.  Lett.  86,  211910  (2005). 

“Substitutional  diatomic  molecules  NO,  NC,  CO,  N:,  and  O?:  Their  vibrational  frequencies  and  effects 
on  p  doping  of  ZnO  ” 

2)  Appl.  Phys.  Lett.  86,  151910  (2005). 

“ Resolving  hydrogen  binding  sites  by  pressure  -  A  first-principles  prediction  for  ZnO  ” 

3)  Phys.  Rev.  B  71,  125209  (2005). 

“Properties  of  Ga-interstitial  defects  in  AlxGai.xNyPi.y” 

4)  Phys.  Rev.  Lett.  95,  225  502  (2005). 

“ Evidence  for  native-defect  donors  in  n-type  ZnO ” 

5)  Physica  B  376-377,  686  (2006). 

“ Probing  deactivations  in  Nitrogen  doped  ZnO  by  vibrational  signatures:  A  first  principles  study  ” 

6)  Physica  B  376-377,  583  (2006). 

“Ion  relaxation  and  hydrogen  LVM  in  H-irradiated  GaAsN” 

7)  J.  Cryst.  Growth  287,  94  (2006). 

“Impurity  effects  in  ZnO  and  nitrogen-doped  ZnO  thin  films  fabricated  by  MOCVD” 

(All  publications  are  acknowledged  this  AFOSR/AOARD  funding  in  the  acknowledgement  section.) 
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(Following  work  schedule  and  milestones  are  the  same  as  previously  submitted) 

Work  schedule: 

Month  0: 

o  Start  state-of-the-art  calculations  for  common  native  defects  in  GaP  and  GaAs 
using  ultrasoft  pseudopotential  (USPP)  and  large  supercell. 

o  Investigate  the  first  principles  all  electron  codes  in  detail.  Designing  viable 
numerical  approaches  to  calculate  isotropic  hyperfine-parameters  for  isolated 
defects  using  the  supercell  method.  Develop  the  tools  to  streamline  the  results  from 
USPP  calculations  to  all  electron  calculations. 

o  Study  N0-Znj  complexes  and  other  point  defects  in  ZnO  using  USPP  calculations 
(in  collaboration  with  experimentalists  at  AFRL/WPAFB) . 

Month  6: 

o  Explore  the  formation  of  defect  complexes  in  GaP  and/or  GaAs  using  USPP 
calculations.  Identify  potential  complexes  for  further  study. 

o  Systematically  calculate  the  isotropic  hyperfine  parameters  for  important  native 
defects  in  GaP  and/or  GaAs. 

o  Study  the  stability  of  Frenkel  pairs  in  GaN  and/or  ZnO  using  USPP  calculations  (in 
collaboration  with  experimentalists  at  AFRL/WPAFB). 

Milestones: 

6  month: 

o  Complete  USPP  results  on  common  native  defects  in  GaP  and  GaAs. 

o  Ability  to  calculate  isotropic  hyperfine  parameters  for  point  defects  in 
semiconductors. 

o  Results  on  stability  of  No-Zn,  complexes  in  ZnO  and  preliminary  results  on  their 
electronic  properties. 

12  month: 

o  USPP  results  on  important  defect  complexes  in  GaP  and  GaAs. 

o  Complete  results  on  isotropic  hyperfine  parameters  for  common  native  defects  in 
GaP  and  GaAs. 

o  Results  on  stability  of  Frenkel  pairs  in  GaN  and/or  ZnO  and  preliminary  results 
on  their  electronic  properties. 
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Method  and  Theory 


Computational  method 

We  performed  our  studies  using  a  first-principles  calculation  approach,  based  on  density- 
functional  theory  and  ab  initio  ultrasoft  pseudopotentials  (USPP).  We  used  the  codes 
called  Vienna  Ab-initio  Simulation  Package  (VASP)  [R1  ].  For  calculations  of  hypcrfine 
parameters,  the  electron  density  near  atomic  cores  is  needed.  In  such  cases,  we  use  all 
electrons  (AE)  calculations  based  on  Augmented  Plane  Wave  method  as  implemented  in 
the  WIEN2k  codes  [R2].  Since  AE  calculations  are  quite  demanding  in  tenns  of 
computational  resources,  we  limited  the  use  of  WIEN2k  calculations  to  only  the  final 
stage  after  the  extensive  atomic  relaxation  steps  are  performed  using  the  less  demanding 
yet  accurate  USPP  codes  (VASP).  The  detail  of  the  parameters  and  approximations  used 
in  the  calculations  is  slightly  varied,  depending  on  specific  material.  Further  specific 
computational  detail  for  each  system  under  study  can  be  found  in  each  attached 
publications. 

A  brief  description  of  the  defect  theory 

Calculations  for  defects  in  semiconductor  are  calculated  using  a  supercell 
containing  32  to  96  atoms  and  a  set  of  non-T  k-points  sampling. 

In  supercell  calculations,  the  formation  energy  of  a  defect  X  in  charge  state  q  is 
defined  as 

Ef[X<]  =  EtotSC(Xq)  -  EtotSC(bulk)  +  nafia  +  qEF,  [1] 

where  E^sciX1)  is  the  total  energy  of  a  supercell  containing  the  defect,  Etot,sc(buIk)  is  the 
total  energy  of  the  same  supercell  but  without  any  defect.  na  is  the  number  of  the  atoms 
of  specie  a  being  removed  from  (positive  value)  or  added  to  (negative  value)  a  bulk 
supercell  to  fonn  defect  X.  /ua  is  the  reservoir  energy  of  specie  a  (the  chemical 
potential).  Ep  is  the  Fenni  level,  referenced  to  the  valence-band  maximum  of  bulk  (no 
defect),  Ev. 

The  chemical  potentials  in  Eq.  [1]  depend  on  the  growth  conditions.  For  example, 
(in  the  case  of  defects  in  ZnO)  under  Zn-rich  case,  /.//n  =  formation  energy  of  solid  metal 
Zn.  In  order  for  ZnO  to  be  stable  during  the  growth,  we  also  require  that  //Zn  +  /u o  =  jUzno 
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=  constant  ®.  This  fixes  /uo  =  ®  -  pz n.  The  chemical  potential  of  other  elements  are 
referenced  to  their  natural  phase.  For  example  the  chemical  potential  of  N,  //\  is 
referenced  to  the  N2  gas  precipitation  limit,  i.e.  ju^=  Etot(N2)/2. 

Research  Activities  on  ZnO 

ZnO  is  a  wide-band-gap  semiconductor  (3.4  eV).  Most  of  its  current  applications, 
including  varistors,  phosphors,  piezoelectric  transducers,  and  transparent  conduction 
films,  do  not  directly  benefit  from  its  wide-band-gap  properties.  Large-area  high  quality 
bulk  ZnO  crystal  could  be  produced  by  various  growing  methods.  This  combining  with 
the  recently  reported  p- type  ZnO  (as-grown  ZnO  is  //-type),  opens  up  the  possibility  of 
producing  blue-  and  UV-range  opto-electronics  devices.  Currently  Nitrogen  is  one  of  the 
most  promising  p-type  dopants  in  ZnO  with  limited  success.  To  aid  in  understanding  ip- 
typo)  doping  difficulty,  we  have  investigated  several  donor  defects  and  defect  complexes 
based  on  first  principles  calculations.  We  are  particularly  interested  in  the  process  that 
deactivated  N  acceptors,  i.e.  in  the  substitution  No  form.  To  accomplish  this 
investigation,  we  have  calculated  a  wide  range  of  possible  defects  and  defect  complexes 
involving  No  and  native  defects  or  predominant  impurities  that  could  be  the  end  results  of 
the  deactivation  process.  Up  to  this  point,  we  found  three  types  of  N0  complexes  that  are 
likely  to  fonn  1)  Zni-No  complex,  2)  No-H  complexes,  and  3)  Substitutional  dimer 
molecules  (SDM)  such  as  (N2)o  and  (NC)o-  Appropriate  signature  of  each  complex  is 
calculated  and  compared  with  experimental  measurements.  Since  H  has  been  identified 
as  an  exclusive  donor  that  could  impede  p-doping  of  ZnO,  we  studied  isolated  H  in  ZnO 
as  well.  Based  on  first  principles  calculations,  we  found  that  under  applied  pressure, 
local  vibrational  mode  of  H  at  different  sites  in  ZnO,  namely  antibonding  (AB)  and  bond- 
center  (BC)  sites,  shifts  differently.  We,  therefore,  proposed  that  IR  experiment  under 
high  pressure  could  be  very  beneficial  in  identifying  the  location  of  H  in  ZnO. 

Zni-No  complex 

Zn  interstitial  (Zni)  has  been  predicted  to  have  high  fonnation  energy  in  //-type 
growth  conditions  [R3].  However,  with  its  (2+)-charge  state,  the  formation  energy  goes 
down  twice  as  fast  as  the  Fenni  level  moves  down.  While  Zni  might  not  be  the 
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predominant  cause  of  /2-type  in  undoped  cases,  it  could  exist  in  a  substantial  amount  in  N- 
doped  cases  where  the  Fermi  energy  is  shifted  down  from  the  CBM  considerably.  In 
addition,  to  place  N  on  the  oxygen  site,  an  O-poor  (or  Zn-rich)  growth  condition  is 
generally  used.  This  also  promotes  Zni  formation.  During  a  cool  down  process,  Zni  can 
passivate  N0  acceptor,  forming  Zni-N0  complex  which  is  a  single  donor.  To  investigate 
the  complex  as  well  as  finding  a  binding  energy,  we  separately  studied  N0,  Zni  and  the 
Zni-No  complexes. 


Fig  1.  Calculated  defect  formation  energy  as  a 
function  of  the  Fermi  level  for  N0,  Zn:,  and  the 
ZnrN0  complex  under  the  Zn-rich  and  N2 
precipitation  limits.  The  Fermi  level  is  references 
to  the  top  of  the  VBM.  The  dotted  lines  show  the 
sum  of  the  formation  energy  of  isolated  Zn\  and 
isolated  N0.  The  calculations  show  that  Zni  binds 
with  N0  to  form  the  ZnrN0  with  the  binding 
energy  of  0.9  eV  (the  energy  difference  between 
the  thick  black  line  and  the  thick  brown  line). 


The  N  atom  stays 


Isolated  substitutional  No  has  four  Zn  nearest  neighbors. 


almost  exactly  on  the  Tci  site  with  a  small  contraction  of  Zn-N  bond  comparing  an  an 
average  Zn-0  bond,  i.e.  contracted  by  3.6  %  and  2.6  %  for  the  neutral  and  (l-)-charge 


state,  respectively.  The  formation  energy  of  No  under  Zn-rich  condition  and  N2 
precipitation  limit  is  shown  in  Fig.  1 . 


Zn  interstitial  prefers  a  site  near  an  ideal  octahedral  site  (in  the  middle  of  a 
hexagonal  channel).  However,  a  strong  repulsion  between  the  interstitial  atom  and  lattice 


cations  (Zn)  causes  the  Zni  atom  to  shift  vertically  to  be  roughly  in  the  middle  between 
two  [0001]  Zn  planes.  The  relaxation  occurs  such  that  neighboring  anion  lattice  atoms 


(O)  are  attracted  to  the  interstitial  (due  to  Coulomb  attraction)  and  are  drawn  inward 


whereas  cation  lattice  atoms  (Zn)  are  repelled  by  the  interstitial  and  are  pushed  outward. 


Zni  is  stable  in  (2+)-charge  state  under  p- type  conditions.  Our  calculations  indicate  that 


the  transition  from  2+/0  could  occur  at  EF  =  2.3  eV  [s(2+/0)  =  2.3  eV].  Since  LDA 
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underestimates  ZnO  bulk  band  gap  by  a  rather  large  amount,  with  band  gap  corrections 
the  s(2+/0)  level  would  lie  very  well  above  the  CBM,  implying  that  Zni  is  a  shallow 
donor  in  agreement  with  previous  calculations  [R4]  [R5]. 

Zni  and  No  can  fonn  a  donor-acceptor  pair  and  become  a  single  donor  complex. 
The  complex  is  most  stable  in  a  (l+)-charge  state  and  has  a  calculated  donor  level  s(+/0) 
=  2.7  eV  above  the  VBM.  Since  this  donor  level  is  higher  than  the  calculated  s(2+/0) 
level  of  isolated  Zni,  which  has  been  concluded  to  be  shallow  [R4][R5],  the  complex  can 
therefore  be  argued  to  be  a  shallow  donor  as  well.  The  binding  energy  of  the  complex  is 
~  0.9  eV,  which  can  be  seen  in  Fig.l.  Our  calculated  fonnation  energies  of  No  (solid 
line),  Zni  (solid  line)  and  No-Zni  complex  (thick  solid  line)  as  a  function  of  EF  are 
illustrated.  The  sum  of  the  formation  energies  of  No  and  Zni  is  plot  with  a  dotted  line. 
The  binding  energy  is  simply  a  total  energy  lowered  when  the  complex  is  formed  (the 
difference  between  the  dotted  line  and  the  solid  line). 

A  brief  explanation  of  Zni-N0  complex  study  is  published  in  conjunction  with  the 
photo  luminescence  results  by  Dr.  Look  in  Phys.  Rev.  Lett.  [P4].  The  full  detail  of  the 
theoretical  part  is  under  preparation  for  publication  (tentatively  for  Phys.  Rev.  B). 


Np-H  complex 


Fig  2.  Formation  energy  of  an  interstitial 
H,  N0  (dashed  lines)  and  a  N0-H 
complex  (solid  line)  in  ZnO  as  a 
function  of  electron  Fermi  energy.  The 
dotted  line  shows  the  sum  of  the 
formation  energies  of  an  isolated 
interstitial  hydrogen  and  an  isolated  N0. 
The  Zn-rich  condition,  N2  and  H2  phase 
precipitation  limits  were  assumed.  The 
energy  difference  between  the  dotted 
line  and  the  solid  line  is  the  binding 
energy  of  the  N0-F1  complex. 


Hydrogen,  which  is  present  in  various  growth  techniques,  could  easily  be 
incorporated  in  the  ZnO  film  during  growth.  Hydrogen  is  known  to  be  amphoteric  in 
most  semiconductors,  but  is  identified  as  an  exclusive  donor  in  ZnO  and  can  passivate 
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acceptors  such  as  No  [R7][R8].  Hydrogen  and  nitrogen  can  interact  strongly  to  form  a 
neutral  defect  complex,  thereby  reducing  the  maximum  achievable  hole  concentration. 
Based  on  first  principles  calculations,  we  calculated  the  formation  energies  of  an  isolated 
interstitial  H,  a  nitrogen  acceptor  No,  and  a  No-H  complex  in  ZnO.  An  isolated 
interstitial  H  is  exclusively  a  donor  in  ZnO  and  always  exists  in  the  fonn  of  H+.  For  an 
isolated  H 1  location,  there  are  four  low-energy  sites  surrounding  an  O  atom  with  the  so- 
called  BCj|  site  being  the  site  with  the  lowest  energy  [R9][R10].  An  isolated  No  is  an 
acceptor.  Currently  the  exact  location  of  the  ionization  energy  of  No  is  still  under  debate 
[Rll].  Over  almost  the  entire  Fermi  energy  range,  the  isolated  No  is  stable  in  a  (1-)- 
charge  state  whereas  an  isolated  interstitial  H  atom  is  always  stable  in  a  (l+)-charge  state. 
The  two  defects,  therefore,  have  a  Coulomb  attraction  and  have  a  strong  tendency  to  form 
a  No-H  defect  complex.  The  formation  energy  of  H+  (at  BC||  site)  and  No  are  plotted 
(dashed  lines)  as  a  function  of  Fermi  level  in  Fig.  2  with  the  sum  of  the  two  fonnation 
energies  shown  as  a  dotted  line.  In  the  same  figure,  we  plotted  the  fonnation  energy  of 
the  No-H  complex  (thick  solid  line).  The  binding  energy  of  the  No-H  complex  is  0.95 
eV,  which  is  the  difference  in  energy  between  the  dotted  and  solid  line.  For  the  No-H 
defect  complex,  the  ABN1  configuration  has  the  lowest  formation  energy.  To  compare 

with  experiment,  we  calculated  the  infrared  absorption  signature  of  this  complex  and 
obtained  the  characteristic  N-H  stretch  frequency  of  3070  cm'1.  This  is  in  a  reasonable 
agreement  with  recently  observed  peak  from  an  MOCVD  sample  at  3020  cm’1  [R1 1]. 

More  detail  of  this  work  can  be  found  in  the  attached  manuscript  (published  in 
Physica  B)  [P5]. 

Substitutional  dimer  molecules  (SDM)  such  as  (N2)o  and  (NC)o 

Based  on  total  energy  calculations,  we  found  that  carbon  impurities  and  excess  N 
strongly  prefer  to  passivate  No  and  form  low-energy  SDM  on  the  Oxygen  site,  (NC)o  or 
(No)o,  both  of  which  are  donors  with  several-eV  binding  energies.  These  SDMs  have 
electronic  and  structural  properties  similar  to  free  diatomic  molecules  (Lee  et  al.  has 
previously  illustrated  for  the  case  of  (N2)o  [R12]). 


Fig  3.  Formation  of  (left  panel)  (NC)o  from  isolated  N0  and  (CO)0  and  (right  panel)  (N2)o 
from  N0  and  (NO)o-  The  Zn-rich  condition,  N2  and  diamond  phase  precipitation  limits  were 
assumed. 

Isolated  interstitials  N;  and  C,  have  very  high  fonnation  energies,  they  tend  to 
immediately  bind  with  the  closest  lattice  O,  forming  (NO)o  and  (CO)o.  Figure  3  shows 
the  formation  energies  as  a  function  of  Fermi  energy  for  (CO)o,  (NC)o,  (NO)o,  and  (N2)o- 
However,  in  the  presence  of  No,  (NO)o  and  (CO)o  are  unstable  against  the  formation  of 
(N2)o  and  (NC)o-  Typical  binding  energies  are  several  eV  due  in  part  to  Coulomb 
attraction  between  oppositely  charged  impurities  and  in  part  to  the  combination  of  two 
impurities  into  one.  In  the  right  panel  of  Fig.  3,  the  formation  energy  of  N0  and  (NO)o 
are  shown  as  a  dashed  line  and  a  thin  solid  line  with  the  sum  of  the  two  as  a  dotted  line. 
The  sum  is  to  be  compared  with  the  thick  solid  line,  which  is  the  fonnation  energy  of 
(N2)o  defect.  We  can  immediately  see  that  the  binding  energy,  which  is  the  difference 
between  the  dotted  line  and  the  thick  solid  line,  depends  somewhat  on  the  Fermi  energy, 
but  has  a  typical  value  of  several  eV.  The  binding  energy  is  of  similar  magnitude  for  the 
case  of  (NC)o,  as  shown  in  the  left  panel.  The  formation  of  SDMs  are  enhanced  in  the 
case  of  p- type  ZnO,  for  which  the  binding  energies  are  large  enough  to  lower  the  energy 
of  both  (N2)o  and  (NC)o  below  the  isolated  N0.  It  is  therefore  expected  that  the 
fonnation  of  (N2)o  and  (NC)o  will  compete  with  the  formation  of  No  acceptors.  Since 
these  SDMs  are  donors  in  p- type  samples,  their  fonnation  further  hinders  /7-type  doping 
by  compensating  the  already  lowered  N0  acceptors. 


2+  _i_ 

Calculated  frequencies  for  (N2)o  and  (NC)o  ,  which  are  the  most  stable  charge 
states  in  p- type  ZnO,  are  2108  and  1995  cm'1,  respectively.  These  frequencies  fall  within 
200  cm'1  of  the  experimentally  observed  values  [R13].  In  addition,  recent  XPS 
experiments  have  found  strong  signals  of  N-N  and  C-N  bonds  with  similar  characteristics 
to  those  of  free  diatomic  molecules  [R14]. 

These  results  suggest  that  one  should  be  careful  with  a  C  impurity  that  can  turn  an 
existing  N  acceptor  (N0)  into  an  (NC)o  SDM,  which  is  a  donor.  In  addition,  excessive  N 
incorporation  could  lead  to  an  unwanted  (N2)o  SDM,  which  is  a  double  donor,  instead  of 
the  desired  N0  acceptor. 

More  detail  of  this  work  can  be  found  in  an  attached  manuscripts  [P1][P5]. 

Research  Activities  on  III-V 

III-V  semiconductors  are  used  extensively  in  electronic  and  optoelectronic 
devices,  ranging  from  GaAs  (1.42  eV)  to  wide  band  gap  semiconductor  GaN  (3.4  eV). 
Recently,  there  are  substantial  interests  on  diluted-nitride  alloys,  for  instant  GaAsi_xNx, 
based  on  their  interesting  characteristics.  Due  to  the  large  size  mismatched  between  N 
and  As  (or  P),  a  small  percents  of  N  can  drastically  change  the  electronic  and  optical 
properties  of  the  host  GaAs  or  GaP  samples.  For  example,  the  small  amount  of  N  can 
dramatically  reduce  the  band  gap  of  GaAs.  During  the  course  of  the  program,  we  have 
three  main  themes  of  research  related  to  III-V  semiconductors,  a)  H  is  known  to 
neutralize  the  effects  of  N  in  diluted-nitride  alloys  back  to  nonnal  bulk  properties.  This 
includes  the  recoveries  of  band  gap,  electron  effective  mass,  and  lattice  constant.  We 
have  investigated  several  possible  N-H  complex  configurations  that  are  likely  to  form  in 
GaAs.  Based  on  comparison  with  experiments,  N-2H  complex  emerges  as  the  top 
candidate  that  can  explain  how  H  neutralize  the  effects  of  N.  b)  Native  defects  in 
traditional  semiconductors,  GaP  and  GaAs,  are  not  well  theoretically  studied.  We  have 
systematically  calculated  simple  native  defects  in  GaP  that  should  provide  an  important 
foundation  for  further  study  of  defects  in  these  III-V  materials,  c)  Recent  optically 
detected  magnetic  resonance  (ODMR)  measurements  found  ODMR  lines  that  correspond 
to  Ga,  but  has  much  weaker  Fermi  contact  values.  We  believed  that  the  observed  FC 
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came  from  a  Ga;  which  is  a  part  of  a  defect  complex.  To  aid  in  this  identification,  we 
have  explored  some  of  the  complexes  that  could  explain  the  observed  smaller  FC  values. 

N-2H  complex  in  GaAsN 

Large  size-mismatched  dilute  alloys  such  as 
GaAsN  and  GaPN  are  promising  alloy  semiconductors 
with  unique  physical  properties.  For  example,  with  only 
an  atomic  percent  of  N  incorporated  into  GaAs,  the  band 
gap  decreases  by  several  hundred  meV  [R15].  On  the 

other  hand,  H  irradiation  causes  a  nearly  full  recovery  of 

.  Ball-stick  model  for  N-2H 

the  GaAs  band  gap  and  lattice  parameter  from  those  ot  compiex 

GaAsN  [R16].  The  fonned  H  complexes  after  H  irradiation  should  be  charge  neutral  due 
to  their  large  concentration  (>10“  cir).  While  there  are  existing  complex  models,  such 
as  a  -  FI*  (N)  complex  that  can  explain  the  experimental  observations,  a  recent  infrared 
experiment  [R17]  on  FI-irradiated  samples  did  not  observe  any  Ga-H  modes,  hence, 
contradicting  the  cr-H^N)  model.  Based  on  first-principles  calculations,  we  show  that 

the  hydrogen  configurations  in  GaAsN  are  actually  depended  on  how  hydrogen  is 
introduced  into  the  sample.  Since  proton  and  neutral  H  have  different  ground  states,  the 
proton  injected  into  the  sample  by  H-irradiation  follows  a  unique  energy  pathway  to  form 
a  charged  dihydride,  instead  of  the  charge-neutral  Hj  monohydride.  The  subsequent 

charge  neutralization  causes  the  spontaneous  canting  of  the  dihydrides.  The  resulting 
canted  N-2H  structure  explains  the  recent  puzzling  IR  observation,  the  recoveries  of  the 
GaAs  band  gap  and  lattice  parameter,  and  the  dihydride  symmetry  determined  by  the 
XANES  experiment. 

More  detail  of  this  work  can  be  found  in  the  attached  manuscript  (published  in 
Physica  B)  [P6] . 
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Native  Defects  in  GaP 


We  have  investigated  fundamental  native  defects  in  GaP,  including  Ga,,  P,,  UGa, 
Vp,  Gap,  and  PGa  using  supercell  approach  with  32-atom  supercell  size.  The  fonnation 
energies  of  defects  are  referenced  to  an  fcc-Ga  (plcc.Ga  =  -2.90  eV/atom)  and  a  P2-dimer 
( |Up2-dimer  =  -4.52  eV/atom)  phases.  Our  calculated  GaP  heat  of  fonnation  is  pGap  =  -9.18 
eV/pair. 


Fig  4.  Formation  energies  of 
native  defects  in  GaP  under  Ga- 
rich  conditions.  Solid  lines 
indicate  EPR  active  states  (the 
defect  state  is  half  occupied). 
Under  Ga-rich  condition,  Ga;  has 
the  lowest  energy  in  p-type 
samples  whereas  GaP  has  the 
lowest  energy  in  /?-type  samples. 


Fig  5.  Formation  energies  of 
native  defects  in  GaP  under  P-rich 
conditions.  Solid  lines  indicate 
EPR  active  states  (the  defect  state 
is  half  occupied).  Under  P-rich 
condition,  PGa  has  the  lowest 
energy  in  p-type  samples  whereas 
UGa  has  the  lowest  energy  in  77- 
type  samples. 


Under  Ga-rich  condition,  the  calculated  formation  energies  of  native  defects  are 


shown  in  Fig.4.  As  a  result  of  high  Ga  partial  pressure,  Ga;,  GaP,  and  Vp  become  low- 


energy  defects.  Under  p- type  condition  (where  Fermi  energy,  EF,  is  near  the  VBM),  Ga, 


in  the  (3+)-charge  state  (triple  donor)  has  the  lowest  formation  energy.  This  means  it 
would  serve  as  the  leading  native  donor  as  one  tries  to  dope  GaP  p- type.  Ga  interstitial 
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prefers  to  stay  in  a  tetrahedral  site  Tci  (also  called  a  cage  site)  in  GaP.  In  the  zincblende 
structure,  there  are  two  types  of  Tj  sites,  one  is  surrounded  by  four  Ga  (7yoa)  and  the 
other  is  surrounded  by  four  P  (T^p).  Ga;  at  both  sites  have  similar  fonnation  energy. 
However,  only  the  Ga;  at  Td,v  has  a  stable  (2+)-charge  state  which  is  EPR  active.  As  the 
Fermi  energy  goes  up,  the  Ga  antisite,  Gap,  which  is  an  acceptor  has  the  lowest  formation 
energy  among  the  the  native  defects.  With  its  (2-)-charge  state,  it  will  counteract  the 
attempt  to  dope  GaP  n-type.  Note  that,  our  reference  for  phosphorus  chemical  potential 
is  somewhat  too  high  (we  used  P2  instead  of  black  phosphorus),  the  actual  position  of  GaP 
should  be  around  1  eV  higher  in  relative  to  other  defects  in  the  plot. 

Under  P-rich  condition,  calculated  formation  energies  of  native  defects  are  shown 
in  Fig. 5.  As  a  result  of  low  Ga  partial  pressure  and  high  P  partial  pressure,  P(,a,  U,a,  and 
Pi  become  low-energy  defects.  Under  p- type  conditions,  P(ia  in  the  (2+)-charge  state 
(double  donor)  has  the  lowest  formation  energy.  This  means  it  would  serve  as  the 
leading  native  donor  as  one  tries  to  dope  GaP  p- type.  At  higher  Fermi  energy,  the  Ga 
vacancy  ( foa),  which  is  an  acceptor,  has  the  lowest  formation  energy.  With  its  (3-)- 
charge  state,  it  will  counteract  the  attempt  to  dope  GaP  n-type.  P,  has  rather  high 
fonnation  energy  and  is  unlikely  to  exist  in  a  noticeable  concentration  under  usual  growth 
conditions. 

To  aid  in  identification  of  these  native  defects,  we  plan  to  perfonn  calculations  of 
hyperfine  parameters  for  each  important  defect  explained  above.  At  this  point,  we  have 
finished  calculating  isotropic  hyperfine  parameters  (Fermi  contact,  FC)  of  three  out  of 
four  important  native  defects  in  GaP,  i.e.  P(,a,  GaP,  and  Ga;. 
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Table  I.  Isotropic  hyperfine  parameters  of  selected  native  defects  in  GaP. 


Central 

Nearest  Neighbor 

Defect 

Atom 

P 

AxlO”4 

Atom 

P 

AxlO”4 

(#  atoms) 

spin/A3 

cm”1 

(#  atoms) 

spin/A3 

cm”1 

P+ 

rGa 

P(l) 

9.71 

984 

P(4) 

0.51 

53 

Gap 

Ga(l) 

0.09 

5 

Ga(4) 

1.44 

76 

GarK.oJ 

Gai(l) 

35.40 

1876 

Ga(4) 

1.55 

82 

Gaf+[7;p] 

Gai(l) 

30.18 

1599 

P(4) 

0.39 

35 

Phosphorus  antisite  (PGa)  can  exist  in  three  charge  states,  2+,  1+,  and  neutral 
depending  on  the  Fermi  level  as  shown  in  Fig. 5.  Out  of  the  three  charge  states,  only  the 
(l+)-charge  state  has  an  unpaired  electron.  Our  calculated  FC  for  PGa  at  the  central  P 

atom  and  at  the  four  neighboring  P  atoms  are  984  x  1 0  4  and  53  x  10  4  cm'1,  respectively. 
These  are  in  good  agreement  with  earlier  measurement  by  Kaufmann  et  al.  who  reported 
the  corresponding  values  of  (966  +  13)xl0”4  and  (81.5  +  4)  xlO”4  cm'1,  respectively 
[R18]. 

Gallium  antisite  (Gap)  can  be  stable  in  three  charge  states:  neutral,  1-,  and  2-  as 
shown  in  Fig. 4.  Only  (l-)-charge  state  has  an  unpaired  electron.  We  found  that  the 
unpaired  electron  is  not  localized  on  the  central  Ga  atom  but  is  distributed  out  to  four 
nearest  neighbors  (Ga)  atoms.  The  spin  density  at  the  nucleus  of  the  central  Ga  atom  is 
only  0.09  spin/A3.  This  is  very  small  in  comparison  to  the  spin  density  of  1.44  spin/A3  at 
the  nuclei  of  the  nearest  neighbor  Ga  atoms.  These  spin  densities  lead  to  FC  values  of 
5x10  4  and  76  x  1 0  4  cm'1  for  the  interaction  at  the  central  Ga  atom  and  at  the  four  nearest 
neighbor  Ga  atoms,  respective. 

Gallium  interstitial  (Ga;)  prefers  to  stay  at  a  tetrahedral  site.  While  the  formation 
energies  are  similar  for  Ga;  at  both  Tj, Ga  and  7//.p  sites,  only  the  later  has  a  stable  (2+)- 
charge  state  (which  has  an  unpaired  electron).  The  (2+)-charge  state  for  G^  at  TdyG a  is 
metastable.  We  calculated  the  FC  of  Ga;  at  both  sites  and  found  that  an  unpaired  electron 
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is  strongly  localized*  at  the  central  interstitial  Ga  atom  with  charge  densities  of  35.40  and 
30.18  spin/A3  for  an  interstitial  at  T(/c,a  and  7yP,  respectively.  These  correspond  to  large 
FC  values  of  1876x10  4  and  1599xl0-4  cm4.  Recent  optically  detected  magnetic 
resonance  (ODMR)  measurements  on  AlGaPN  by  our  collaborators  (the  group  of  Prof. 
Weimin  Chen  et  a/.,  Linkoping  University,  Sweden),  found  two  sets  of  pattern  of  ODMR 
lines  that  are  proved  to  come  from  Ga  interstitial  and  labeled  them  Ga,-A  and  Ga,-B. 
However,  the  measured  FC  values  of  both  Ga;  defects  are  much  smaller  than  the 
calculated  values  of  isolated  Ga;,  i.e.  770x1 0-4  cnT1  for  Gaj-A  and  1150  x  1CT4  cnT1  for 
Ga;-B.  We  believed  that  the  observed  FC  came  from  a  Ga;  which  is  a  part  of  a  defect 
complex,  not  an  isolated  Ga;.  This  model  is  supported  by  recent  observations  (both 
theory  and  experiment)  that  Ga;  in  GaN  has  a  low  diffusion  barrier  and  should  diffuse  at 
around  room  temperature.  Since  GaP  is  less  dense  than  GaN,  all  existing  Ga;  should  able 
to  diffuse  easily  and  form  defect  complex  with  other  defects.  We,  therefore,  should  not 
detect  any  isolated  Ga;.  To  aid  in  this  identification,  we  have  explored  some  of  the 
complexes  that  could  explain  the  observed  smaller  FC  values. 

Some  Defect  Complexes  in  GaP 

Our  first  assumption  was  that  an  Al-atom  replacing  lattice  Ga  (Alca)  near  a  Ga; 
might  strongly  reduce  the  FC  value  of  an  isolated  Ga;.  For  this,  we  have  investigated  a 
Ga;  at  the  Tyca  site  (labeled  as  Ga31  \Td  Ga] ).  And  replace  one  nearest  neighbor  Ga  atom 

with  an  A1  atom,  leading  to  a  GarA/Ga  complex.  Since  A1  and  Ga  are  isovalent,  the 
charge  state  of  this  complex  that  results  in  an  unpaired  electron  is  the  same  as  the  charge 
state  of  an  isolated  Ga;,  i.e.,  the  2+  charge  state.  For  a  quick  test,  we  use  the  fully  relaxed 
atomic  structure  of  an  isolated  Ga3  [Td  Cia]  defect  and  just  replace  one  nearest  neighbor 

Ga  atom  with  an  A1  atom.  To  justify  this  assumption,  the  structural  relaxation  test  were 
perfonned  using  a  pseudopotential  codes.  We  indeed  found  only  a  small  atomic 
relaxation.  The  FC  results  show  that  replacing  a  Ga  nearest  neighbor  by  an  A1  atom  leads 
to  a  very  small  reduction  in  the  central  FC  (see  Table  II).  The  FC  value  of  this  complex 


The  impaired  electron’s  density  is  42-49%  of  our  calculated  atomic -Ga  4s  electron.  (  (0)  =  72.7  e"/A3) 
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is  still  way  too  large  in  comparison  to  any  of  the  observed  ones.  Therefore,  we  concluded 
that  Gaj-Aloa  defect  complex  could  not  explain  the  measured  FC. 

Next,  we  aimed  our  attentions  toward  GarGap  complex  since  Ga,  and  GaP  are 
both  have  low  formation  energies  under  Ga-rich  condition.  Ga;  should  serve  as  the 
leading  donor  and  Gap  as  the  leading  acceptor  if  there  are  no  other  lower-energy 
impurities  (see  Fig.  4).  As  the  sample  is  cool  down,  both  complex  are  attracted  to  each 
other  and  can  combine  to  fonn  Ga;-GaP  complex  which  is  a  donor-acceptor  pair. 


Fig  6.  The  formation  of  Gaj-GaP  complex.  The  charge  states  with  an  unpaired  electron  are 
plotted  with  thick  solid  lines  with  boldface  labeling  including  the  charge  state  numbers.  The 
dotted  lines  show  the  summation  of  the  formation  energy  of  the  isolated  defects  forming  the 
complex. 

To  minimize  the  donor-acceptor  pair  distance,  the  Ga;  binds  to  the  GaP  on  the  Td 
site  with  (originally  before  the  Ga  atom  substitute  on  the  P  site)  4  P  atoms.  Since  a  P 
atom  is  substituted  by  a  Ga  atom,  the  Ga;  is  left  with  three  P  and  a  Ga  nearest  neighbor 
atoms.  Substituting  a  neighboring  group  V  atom  with  a  group  III  atom  results  in  a  loss  of 
two  electrons.  Therefore,  the  charge  state  of  the  defect,  that  has  an  unpaired  electron, 
changes  from  2+  in  an  isolated  Ga;  case  to  neutral  in  this  Ga;  -Gap  complex.  From  an 
energetic  viewpoint,  we  find  that  the  complex  is  “bound”  with  a  binding  energy  of  about 
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0.5  eV  (the  energy  different  between  the  dotted  line  and  the  energy  of  Ga;  -Gap 

complex  in  Fig.  6).  The  formation  energy  of  each  constituent  defect  and  the  complex  is 
plotted  as  a  function  of  the  electron  Fenni  level  in  Fig.  6.  Calculated  FC  values  of 
(Ga;  -Gap)°  complex  are  tabulated  in  Table  II.  The  FC  value  at  the  central  Ga;  atom  is 

1 123  x  10  4  cm  1  matched  well  with  an  observed  Ga;-B  of  1 150  x  10~4  cm  '.  Although,  N 
is  not  directly  involve  in  this  complex  formation,  the  growth  conditions  needed  to 
incorporate  N  on  the  P  sites  should  lead  to  a  P-poor  condition  (or  another  words  Ga-rich 
condition)  that  also  promotes  the  formation  of  both  constituents,  i.e.  Ga;  and  Gap,  of  the 
complexes. 


Table  II.  Isotropic  hyperfine  parameters  of  Ga,  and  selected  Gaj-complexes  in  GaP. 


Defect 

Central 

Nearest  neighbor 

Nearest  Neighbor 

atom 

P 

spin/A3 

Axl0~4 

cm-1 

atom 

P 

spin/A3 

Axl0~4 

cm-1 

Atom 

P 

spin/A3 

Axl0~4 

cnT1 

Gar[^GJ 

Ga,(  1 ) 

35.40 

1876 

Ga(4) 

1.55 

82 

(Ga;  -AlGa)2+ 

(unrelaxed) 

Ga,(  1 ) 

34.07 

1805 

AlGa(l) 

0.36 

21 

Ga(3) 

1.58 

84 

Gaf+[r,p] 

Ga,(  1 ) 

30.18 

1599 

P(4) 

0.39 

35 

(Ga;  -Np)2+ 

Ga,(  1 ) 

26.30 

1393 

Np(l) 

0.95 

15 

P(3) 

0.52 

47 

(Ga,  -CP)+ 

Ga,(  1 ) 

25.95 

1375 

Cr(l) 

0.26 

15 

P(3) 

0.53 

48 

(Ga,  -Gap)° 

Ga,(l) 

21.20 

1123 

GaP(l) 

0.02 

1 

P(3) 

0.29 

25 

(Ga,-Alp)° 

Ga,(  1 ) 

20.76 

1101 

A1p(1) 

0.08 

5 

P(3) 

0.30 

27 

(Ga,  -2Gap)2" 

Ga,(  1 ) 

17.72 

939 

GaP(2) 

0.37 

19 

P(2) 

0.26 

18 

(Ga,  -3Gap)4~ 

Ga,(  1 ) 

18.77 

995 

GaP(3) 

0.31 

16 

P(l) 

0.08 

7 

To  study  the  cation  effects  on  the  FC  of  Ga;  -Gap  complex,  we  replace  a  Ga 

antisite  (Gap)  with  an  A1  antisite,  creating  Gctj-Alp  complex.  (Full  atomic  relaxation  has 
been  carried  out  here.)  We  found  that  the  effect  of  replacing  the  Gap  by  Alp  is  small,  i.e., 
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resulting  in  a  reduction  of  the  central  FC  by  only  22  x  10  4  cm  1  (see,  Table  II).  This 
could  qualitatively  explain  the  modest  variation  of  the  FC  (Gai-B)  with  A1  composition. 

Since  C  impurity  is  abundant  in  many  growth  methods,  it  might  form  complexes 
with  Ga,.  We,  therefore,  have  investigated  Gcii-Cp  complex.  Since  here  we  replace  a  P 
atom  (group  V)  with  a  C  atom  (group  IV),  therefore  only  one  electron  is  missing.  In 
order  for  the  complex  to  have  an  unpaired  electron  occupied  Ga;  4s,  we  need  (l+)-charge 
state.  This  complex  can  also  be  considered  as  a  of  donor-acceptor  pair  since  Cp  is  a  deep 
acceptor  in  GaP.  We  have  investigated  the  formation  energy  and  found  that  the  complex 
is  bound  with  the  binding  energy  of  about  1  eV.  The  formation  energy  of  each  constituent 
defect  and  the  complex  is  plotted  as  a  function  of  the  electron  Fermi  level  in  Fig  7.  We 
calculated  the  FC  parameters  of  this  complex  and  tabulated  the  values  in  Table  II.  We 
found  the  reduction  in  the  central  FC  of  only  224x1  (T4  cm  1  from  an  isolated  Ga;  atom, 
leading  to  the  FC  value  of  1375xl0-4  cm  .  This  does  not  fit  well  with  any  observed 
values.  Nevertheless,  this  result  gives  a  very  valuable  trend  for  further  study. 

We  can  see  that  if  one  starts  from  an  isolated  Ga;  on  the  Tj  site  surround  by  4  P, 
one  can  reduce  the  FC  by  replacing  the  neighboring  P  by  other  atom,  preferably  from  a 
different  column.  Replacing  a  nearest  neighbor  P  by  a  group  IV  (C)  atom  leads  to  a 
reduction  in  the  central  FC  of  ~220xlCT4  cm  \  Replace  a  nearest  neighbor  P  by  a  group 
III  (either  Ga  or  Al)  atom  leads  to  reduction  in  the  central  FC  of  ~  476  x  10  4  cm-1,  i.e. 
one  gains  ~220xl  O-4  cm-1  reduction  for  each  column  in  the  periodic  table  one  moves 
away.  If  we  assume  that  this  trend  keeps  carry  on,  then  a  (Ga,  -  Fp  )3  complex  should 

have  a  central  FC  of  around  700  x  10  4  cm  ,  hence,  in  agreement  with  Ga;-A 
measurement.  However,  our  calculations  show  that  the  GarVp  complex  is  not  stable. 
The  Ga;  spontaneously  moves  to  replace  the  vacancy  and  form  a  low-energy  GaP  defect. 
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Fig  7.  The  formation  of  Ga-CP  complex.  The  charge  states  with  an  unpaired  electron  are 
plotted  with  thick  solid  lines  with  boldface  labeling  including  the  charge  state  numbers.  The 
dotted  lines  show  the  summation  of  the  formation  energy  of  the  isolated  defects  forming  the 
complex.  The  chemical  potential  of  a  C  atom  is  referenced  to  the  diamond  phase. 

We  found  that  the  Ga,.  -Gap  can  attract  another  Gap  and  form  a  larger  Ga\-2Gap 

complex.  This  complex  has  a  Ga;  in  the  Tci  site  originally  surrounded  by  4  P  but  has  two 
P  atoms  replaced  by  Ga  atoms.  Substituting  two  of  the  group  V  (P)  by  group  III  (Ga) 
atoms  results  in  a  total  loss  of  4  electrons.  Therefore,  the  charge  state  of  the  defect  (that 
leads  to  an  unpaired  electron  in  the  Ga;  4s  state)  changes  from  2+  in  an  isolated  Ga;  case 
to  2-  in  this  Ga;  -2Gap  complex.  From  the  energetic  viewpoint,  we  find  that  an 

additional  GaP  is  “bound”  to  the  Ga,  -Gap  complex  to  fonn  Ga,  -2Gap  with  a  binding 

energy  of  about  0.5  eV.  The  fonnation  energy  of  each  constituent  defect  and  the  complex 
is  plotted  as  a  function  of  the  electron  Fermi  level  in  Fig.  8. 
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Fig  8.  The  formation  of  Ga;-2GaP  complex.  The  charge  states  with  an  impaired  electron  are 
plotted  with  thick  solid  lines  with  boldface  labeling  including  the  charge  state  numbers.  The 
dotted  lines  show  the  summation  of  the  formation  energy  of  the  isolated  defects  forming  the 
complex. 

The  calculated  FC  values  of  (Ga;  -2Gap)2~  complex  are  tabulated  in  Table  II. 
The  FC  value  of  the  central  Ga;  atom  is  939  x  10  4  cm-1,  which  is  reasonably  close  to  the 
observed  Ga,-A  of  770  x  1CT4  cnT1. 

Since  the  Ga;  -2Gap  complex  does  not  lower  the  FC  down  enough  to  match  well 
with  the  observed  Ga,-A.  We  further  investigate  whether  if  we  convert  one  more  P 
neighbor  to  be  an  antisite  would  help  to  further  lower  the  FC.  It  turns  out  that  a 
Gar 3 Gap  complex  in  a  (4-)-charge  state,  which  is  stable  in  our  calculations  and  has  an 
unpaired  electron  in  the  Ga  As  state,  have  the  FC  higher  than  the  Ga;  -2Gap  complex 

(see,  Table  II).  This  is  not  at  all  surprising  if  we  assume  that  the  FC  of  the  complex  is 
reduced  from  the  isolated  interstitial  because  of  the  result  of  symmetry  breaking.  In  this 
view,  the  Ga;  -3Gap  complex  has  similar  symmetry  as  the  Ga;  -Gap  with  the 

interstitial  atom  shifted  further  away  from  the  center  -  that  is  why  a  FC  for  a  Ga;  -3Gap 
complex  is  slightly  lower  than  a  Ga;  -Gap  but  higher  than  a  Ga;  -2Gap  where  the 
symmetry  is  strongly  broken. 
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Fig  9.  The  formation  of  Ga,-NP  complex.  The  charge  states  with  an  impaired  electron  are 
plotted  with  thick  solid  lines  with  boldface  labeling  including  the  charge  state  numbers.  The 
chemical  potential  of  the  N  atom  is  set  such  that  NP  has  zero  formation  energy. 


Our  experimental  collaborators  noticed  that  Ga,-B  signal  depends  rather  strongly 
with  a  N  concentration.  We  studied  effects  of  replacing  an  isovalent  atom  by  starting 
with  Gaf+[rrf  p]  defect  and  replace  a  nearest  neighbor  P  atom  with  a  N  atom.  At  least,  we 

found  that  replacing  an  anion  isovalent  atom  (i.e.,  replace  a  P  atom  with  a  N  atom)  has  a 
larger  effect  on  the  central  FC  than  replacing  a  cation  isovalent  atom  (Ga  by  Al).  This  is 
probably  due,  at  least  in  part,  to  a  large  structural  relaxation  introduced  by  N  atom.  FC 
values  are  tabulated  in  Table  II. 

However,  from  the  energetic  point  of  view,  we  found  the  Ga;  do  not  bind  with  the 
Np.  The  plot  of  each  defect  as  well  as  the  complex’s  fonnation  energy  as  a  function  of 
electron  Fenni  energy  is  shown  in  Fig.  9.  Here  we  set  the  chemical  potential  of  N  atom 
such  that  the  formation  energy  of  Np  is  zero.  (The  binding  energy  is  an  absolute  quantity 
which  does  not  depend  on  the  chemical  potential  of  N.)  The  total  fonnation  energy  of 
individual  defects  is,  therefore,  just  the  formation  energy  of  Ga;.  It  is  clear  that  the 
complex  has  higher  fonnation  energy  by  almost  1  eV  and  should  not  fonn. 
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In  conclusion  (regarding  the  Ga;  FC),  our  calculated  FC  for  an  isolated  Ga^+  in 

GaP  is  too  large  to  be  consistent  with  the  measured  values.  Our  experiment  collaborators 
(Linkopin  Univ.,  Sweden)  observed  two  signatures  via  ODMR  measurements  and  labeled 
them,  Ga;-A  (Fermi  contact  =  770xl(T4  cm”1)  and  Ga;-B  (1150x10  4  cm-1).  Our 
calculated  values  for  isolated  Ga  interstitials  are  1876x10  4  cm-1  for  a  Ga,  at  a  Tj  site 
surrounded  by  four  Ga  atoms  and  1599xl0”4  cm-1  for  a  Ga;  at  a  Tj  site  surrounded  by 
four  P  atoms.  We  proposed  that  the  observed  hyperfine  spectrum  belongs  to  Ga; 
complexes  and  have  investigated  some  of  the  probable  ones.  Indeed,  we  found  that  the 
value  of  Ga;  FC  can  be  lowered  by  fonning  a  complex.  At  this  point,  we  found  a 
(Ga;  -Gap)°  complex  (FC  =  1123xl0”4  cm”1)  and  a  (Ga;  -2Gap)2”  complex  (FC  = 

939x1  CT4  cm”1)  to  be  the  most  probable  models  for  the  observed  Ga;-B  and  Ga;-A, 
respectively.  We  have  published  our  results  together  with  our  experimental 
collaborators  in  Phys.  Rev.  B  [P3]. 


-22- 


Future  Research  Activities 

Defects  and  defect  complexes  in  III-V 

Summarize  the  theoretical  results  on  native  defects  and  some  complexes  in  GaP  for 
publication.  Complete  the  study  of  similar  defects  in  GaAs.  In  addition,  we  will  study 
some  defect  complexes  in  GaAs  that  are  of  current  interests,  such  as  nitrogen-carbon 
dimer  complex  in  GaAs  that  are  recently  identified  by  an  IR  experiment. 

Frenkel  pairs  and  defect  complexes  in  ZnO 

1)  We  will  compose  a  detailed  theoretical  manuscript  that  describes  in  detail  the 
formation  of  the  Zni-No  complex,  which  is  believed  to  responsible  for  the  observed  PL 
signature  by  Dr.  Look.  2)  We  will  investigate  the  stability,  the  formation,  and  the 
recombination  of  Zn,-Vzn  Frenkel  pair  via  first-principle  method.  This  study  will  involve 
the  calculation  of  diffusion  barrier  of  Zn,  and/or  Vzn  in  different  charge  states.  3)  Since 
0;  is  proposed  to  be  an  important  native  acceptor  that  neutralizes  the  75-meV  donors  in 
the  irradiation  experiment,  we  will  investigate  the  electronic  properties  and  formations  of 
the  complexes  involving  O;  and  substitution  A1  and/or  Ga. 
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First-principles  calculations  show  that  AB  defects  substituting  on  an  O  site  in  ZnO  where  A,  B 
=  N,  O,  or  C  are  an  important  class  of  defects  whose  physical  properties  cannot  be  described  by  the 
usual  split  interstitials  but  rather  by  substitutional  diatomic  molecules.  The  molecular  natures  of  the 
(AB)0  defects  are  reflected  in  their  vibrational  frequencies  which  are  redshifted  from  those  of  the 
corresponding  free  molecules  but  only  by  about  10%.  These  calculated  results  agree  with  the 
frequency  range  recently  observed  by  IR  measurement  on  N-doped  ZnO.  Moreover,  most  (AB)0 
defects  are  donors  in  p-type  samples.  The  (NC)o  and  (N2)o  defects  have  sufficiently  low  energies 
to  convert  substituional  N0  acceptors  into  donors,  thereby  hindering  the  efforts  of  doping  ZnO  p 
type.  ©  2005  American  Institute  of  Physics.  [DOI:  10.1063/1.1931823] 


Recent  experimental  investigations  on  chemical  vapor 
deposition  grown  nitrogen-doped  ZnO1  4  revealed  that  these 
samples  contain  not  only  a  large  amount  of  H,  but  also  C,  as 
well  as  an  excess  of  N  in  unwanted  forms.  This,  combined 
with  the  difficulty  of  achieving  p-type  ZnO,  suggests  that 
these  defects  may  play  important  roles  in  compensating  or 
passivating  the  N  acceptors.  Moreover,  IR  absorption  mea¬ 
surements  on  N-doped  ZnO  showed  several  peaks  in  the 
range  of  1800-2000  cm-1,  which  are  surprisingly  close  to 
the  IR  frequencies  of  free  diatomic  molecules  but  signif- 
cantly  higher  than  the  phonon  frequencies  of  the  ZnO  host. 
Recent  x-ray  photoelectron  spectroscopy  (XPS)  experiments 
found  strong  signals  from  the  C-N  and  N-N  bonds,  which 
are  characteristic  of  those  in  free  diatomic  molecules.  In  a 
recent  theory  work,5  charge  distribution  near  a  substitutional 
N2  also  showed  the  sign  of  a  substitutional  molecule. 

In  this  work,  we  investigate  systematically  first-row  C, 
N,  and  O  impurities  in  ZnO  forming  the  so-called  split  inter¬ 
stitials  such  as  (NC)0  and  (N2)0.  First-principles  calculations 
show  that  these  complexes  have  electronic  and  structural 
properties,  which  are  inconsistent  with  the  traditional  split- 
interstitials  picture6  but  in  good  agreement  with  a  substitu¬ 
tional  diatomic  molecule  (SDM)  picture.  In  particular,  the 
calculated  diatomic  bond  length,  wave  functions,  and  stretch 
frequencies  closely  resemble  those  of  free  molecules,  in 
qualitative  agreement  with  experiments.  The  strong  tendency 
of  forming  such  molecules  at  least  in  the  cases  of  (NC)0  and 
(N2)0  also  causes  the  conversion  of  the  N0  acceptors  into 
unwanted  donors  that  further  compensate  the  remaining  ac¬ 
ceptors  in  ZnO. 

Our  calculations  are  performed  using  the  density  func¬ 
tional  theory  with  the  local  density  approximation  (LDA) 
and  ultrasoft  pseudopotentials,7  as  implemented  in  the  VASP 
codes.8  The  Zn  3d  electrons  are  treated  as  valence  electrons. 
The  cutoff  energy  for  the  plane  wave  basis  set  is  300  eV.  The 
calculated  heat  of  formation  of  3.58  eV  agrees  with  the  ex- 
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perimental  value  of  3.60  eV.  We  use  a  supercell  with  36 
atoms  for  the  defect  study.9  Results  for  (N2)0  are  tested  using 
a  96-atom  cell  to  ensure  the  convergence  to  be  +0.2  eV.  For 
charged  defects,  a  jellium  background  was  used.  Since  LDA 
severely  underestimates  the  band  gap,  we  have  instead  ex¬ 
amined  the  electronic  properties  at  the  2X2X2  Monkhorst- 
Pack  special  k  points,  which  are  used  for  the  Brillouin  zone 
integration.  The  band  gap  at  the  special  k  points  is  2.5  eV. 
All  the  atoms  are  relaxed  by  minimization  of  force  to  less 
than  0.05  eV/A. 

The  defect  formation  energy  is  defined  as10 

AEf=Elol(D,q)  -  Etot{ 0)  +  2  Anx/ux  +  qEF, 

where  Etot(D,q)  is  the  total  energy  of  the  supercell  with  de¬ 
fect  D  in  charge  state  q\  £tot(0)  is  the  total  energy  of  the 
supercell  without  the  defect;  A nx  is  the  number  of  species  X 
(=  Zn,  O,  N,  and  C)  being  removed  from  a  defect-free  cell  to 
its  respective  reservoir  with  chemical  potential  ptx  to  form 
the  defect  cell.  The  upper  limits  for  /xZn,  ptQ,  /rN,  and  ptc  are 
the  calculated  energies  of  metallic  Zn,  gaseous  02  and  N2, 
and  diamond.  For  mathematical  simplicity,  however,  we  ex¬ 
press  here  the  chemical  potentials  relative  to  their  respective 
elemental  natural  phases.  This  redefines  the  ptx  so  that  ptx 
«0.  To  keep  ZnO  thermodynamically  stable,  it  is  also  re¬ 
quired  that  /uZn+/A0=/xZn0=-3.58  eV.  EF  is  the  electron 
Fermi  level  with  respect  to  the  valence  band  maximum 
(VBM)  averaged  at  the  special  k  points,  as  described  in  Refs. 
11  and  12. 

First,  let  us  consider  split  interstitial  in  the  “traditional” 
picture:  namely,  a  complex  embedded  in  a  local  sp 3  environ¬ 
ment.  For  example,  the  C2  split  interstitial  in  diamond  where 
each  carbon  atom  is  threefold  coordinated.6  One  can  expect 
that  the  pair  form  single  bonds  with  surrounding  host  atoms. 
Alternatively,  one  can  view  the  complex  according  to  the 
molecular  orbital  theory,  in  which  the  energy  levels  of  free 
molecules  are  sscr,  sscr*,  ppu  (doublet)  pp it,  ppu  (doublet), 
and  ppu  in  the  order  of  increasing  energy.  When  embedded 
in  ZnO,  the  sscr,  sscr* ,  ppu,  and  ppu  states  are  all  in  the 
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FIG.  1.  (Color  online),  (a)  Atomic  structure  of  an  (AB)Q  defect,  A=N,  B 
=C,  where  the  encircled  AB  dimer  replaces  an  O  atom.  Large  sphere  is  Zn 
and  small  sphere  is  O.  (b)  Charge  density  of  the  ppir *  states  of  the  (NC)0 
defect.  The  two  states  should  be  degenerate  in  a  free  dimer  as  shown  in  (c) 
but  the  degeneracy  is  lifted  by  ZnO  crystal  field. 

valence  band  as  expected  and  are  hence  doubly  occupied.  If 
we  do  an  electron  counting  for  (NC)0,  for  example,  we  have 
5  doubly  occupied  states  with  11  available  electrons,  2  from 
4  Zn  neighbors  in  Fig.  1(a),  5  from  N,  and  4  from  C.  This 
leaves  (1 1  - 10)  =  1  electron  to  occupy  the  higher  pp  v"  (dou¬ 
blet)  states  above  the  VBM  as  shown  in  Fig.  2. 

The  calculated  wave  functions  in  Fig.  1(b)  are  almost 
copies  of  the  pp v*  states  for  the  free  dimer  in  Fig.  1(c).  Most 
strikingly,  the  charge  contour  in  Fig.  1(b)  is  for  the  entire 
supercell,  suggesting  that  these  doublet  states  are  highly  lo¬ 
calized.  Indeed,  the  dispersion  in  k  space  due  to  cell-cell 
interaction  through  the  ZnO  host  is  small,  only  0.2  eV  for 
these  states.  Thus,  there  is  little  coupling  between  the  /; /577- 
states  and  the  surrounding  Zn  atoms.  This  shows  clearly  that 
the  SDM  model  is  a  better  description  of  the  (NC)o  com¬ 
plex.  Due  to  the  crystal  field  of  ZnO,  however,  the  double 
degenerate  ppv  states  of  the  free  molecule  is  replaced  by 
two  nondegenerate,  spatially  orthogonal  states. 

The  calculated  bond  lengths  compare  favorably  with 
those  of  the  corresponding  free  molecules,  also.  Note,  how¬ 
ever,  that  as  a  substitutional  impurity,  each  (AB)0  accepts 
two  electrons  from  its  Zn  neighbors.  Thus,  we  should  com¬ 
pare  neutral  free  molecules  with  SDMs  of  (2  +  )  charge  state 
to  keep  the  same  orbital  occupation.  Table  I  shows  that  for 
(N2)q+,  the  bond  length,  g(n_n=1.14A,  is  only  3%  longer 
than  the  triple  bond  of  free  N2  of  1.1 1  A.  In  contrast,  the  sum 
of  single -bond  atomic  radii  yields  1.50  A,  which  is  32% 
longer.  For  (CO)q+,  c/c-o  is  longer  than  that  of  free  CO  by 
only  2%.  For  (NO)q+,  ^n-o  Is  5%  longer.  An  elongation  is 
generally  expected  due  to  the  screening  of  the  ZnO  host. 

To  calculate  stretch  frequencies,  we  employed  the  ap¬ 
proach  used  in  Ref.  13.  Our  calculated  results  for  the  free 
molecules  are  in  reasonable  agreement  with  the  measured 
values:  1485(cm-1)  versus  1580.2  for  02,  1819  versus 


Conduction  band  'd 


(0,)0  <NO)0  (CO)0  (N,)c  (NC)q 


FIG.  2.  Calculated  energy  levels  of  the  ppir*  states  for  charge  neutral  de¬ 
fects  with  respect  to  the  band  gap  of  ZnO  calculated  by  the  special  /,-poml 
average.  Solid  dots  show  the  electron  occupations. 
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TABLE  I.  Calculate  bond  length,  stretch  frequencies,  and  formation  energy 
of  the  (AB)0  defects  at  EF= 0,  p.Zn= pN= p,c= 0.  Jfree  is  the  bond  length  of 
free  molecules. 


Defect 

(AB)0 

Charge 

<1 

^A-B 

(A) 

^free 

(A) 

CO 

(cm-1) 

AEf 

(eV) 

(02)o 

0 

1.46 

1.23 

1051 

5.1 

(NO)0 

2+ 

1.23 

1462 

2.5 

1+ 

1.29 

1357 

3.1 

0 

1.36 

1173 

4.4 

1- 

1.48 

1.17 

967 

6.1 

(CO)0 

2+ 

1.16 

1894 

1.8 

1+ 

1.22 

1652 

3.5 

0 

1.27 

1.14 

1373 

5.7 

(^2)0 

2+ 

1.14 

2108 

-1.5 

1+ 

1.20 

1852 

0.4 

0 

1.28 

1.11 

1537 

2.5 

(NC)0 

1+ 

1.18 

1.17 

1995 

0.4 

1904.2  for  NO,  2197  versus  2358.6  for  N2,  2132  versus 
2169.8  for  CO,  and  2052  versus  2068.6  for  NC.  To  calculate 
the  frequency  u>  for  the  SDMs  in  ZnO,  in  principle,  one 
should  solve  the  dynamic  matrix  with  coupling  to  at  least  a 
few  shells  of  neighboring  atoms.  However,  we  found  that  the 
coupling  to  the  ZnO  host  is  very  weak,  which  is  a  reminis¬ 
cence  of  the  substitutional  molecules:  for  example  for  (N2)0, 
the  coupling  results  in  a  shift  in  &>  less  than  2  cm-1.  Because 
of  this,  in  the  calculation  we  keep  the  host  atoms  fixed, 
which  allows  us  to  include  also  the  anharmonic  effect.  Table 
I,  column  5  shows  to  scaled  by  the  ratio  between  experiment 
and  theory  for  free  molecules  (given  by  the  previous  results) 
to  reduce  systematic  errors  in  the  results.  The  frequencies  of 
the  SDMs  are  reasonably  close  to  those  of  the  free  mol¬ 
ecules. 

Because  the  ppir'  molecular  orbitals  are  antibonding  or¬ 
bitals,  its  occupation  will  lead  to  an  increase  in  the  AB  bond 
length.  We  can  see  from  Table  I  that  adding  electrons  to  the 
ppir  states  always  causes  an  increase  in  the  AB  bond  length, 
resulting  in  weaker  bonds  and  a  reduction  in  oj.  For  example 
for  (N2)0,  the  bondlength  increases  from  1.14  A  for  c/=2+, 
to  1.20  A  for  </  =  +,  and  to  1.28  A  for  q= 0.  Accordingly,  the 
frequency  decreases  from  2108,  to  1852,  and  eventually  to 
1537  cm-1.  Because  of  the  Coulomb  repulsion  between  elec¬ 
trons,  energy  levels  in  Fig.  2  will  shift  when  occupation 
changes.  For  example,  in  the  case  of  (02)0,  the  removal  of 
any  electron  from  the  gap  state  in  Fig.  2  would  lower  the 
level  to  below  the  VBM,  thus  making  the  neutral  state  the 
only  stable  state. 

The  formation  energies  of  the  SDM  impurities  (A Ej)  are 
tabulated  in  Table  I.  Figure  3  shows  the  EF  dependence  of 
A Ef  for  (CO)0,  (NC)o,  (NO)0,  (N2)0,  and  N0.  Our  results 
for  N0  agree  to  within  0.2  eV  with  those  in  Refs.  5,  14,  and 
15.  For  (N2)0,  our  results  agree  to  within  0.3  eV  with  Ref. 
15.  Isolated  interstitial  N,  and  C;  are  unstable  against  spon¬ 
taneous  formation  of  (NO)o  and  (CO)o  so  they  are  ignored. 
In  the  presence  of  N0,  (NO)o,  and  (CO)o  are  also  unstable 
against  the  formation  of  (N2)Q  and  (NC)0.  On  the  other 
hand,  both  (N2)0  and  (NC)0  have  lower  formation  energy 
than  that  of  isolated  N0  in  p-type  ZnO.  Hence,  their  forma¬ 
tion  will  severely  reduce  the  concentration  of  N0  acceptors. 
AIP  license  or  copyright,  see  http://apl.aip.org/apl/copyright.jsp 
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FIG.  3.  Formation  of  (a)  (NC)0  from  isolated  NO  and  (CO)0  and  (b)  (N2)0 
from  N0  and  (NO)0-  The  chemical  potentials  are  set  at  /rZn=/u,N=/u,c=0. 

Also,  the  SDMs  are  donors,  therefore  compensating  any  re¬ 
maining  N0  acceptors. 

The  situation  could  be  more  severe  when  higher  N  con¬ 
centration,  [N],  than  the  equilibrium  value  is  incorporated 
via  nonequilibrium  growth  methods  because  high  [N]  im¬ 
plies  higher  chemical  potential  /j,s.  Note  that  the  formation 
energy  of  the  (N2)o  complex  decreases  twice  as  fast  as  that 
of  N0  when  one  raises  /%.  However,  because  the  formation 
of  (N2)o  from  N0+(NO)o  involves  the  breakup  of  one 
SDM,  the  diffusion  of  Nj,  and  the  formation  of  another 
SDM,  one  might  be  able  to  control  the  growth  temperature  to 
be  low  enough  to  suppress  (N2)0.  Also,  one  should  be  cau¬ 
tioned  about  carbon  impurities  in  p-type  ZnO.  Carbon  can 
form  not  only  double  donor  C/n,  by  substituting  Zn,  but  also 
single  donor  (NC)0,  which  has  low  formation  energy  when 
EF  is  low.  Hence,  beside  controlling  the  concentration  of 
(N2)0  by  kinetic  means,  one  also  has  to  suppress  CZn  and 
(NC)0  in  order  to  get  high  enough  hole  concentrations. 

Recent  XPS  experiment  has  found  strong  signals  of  N-N 
and  C-N  bonds  with  similar  characteristics  to  those  of  free 
diatomic  molecules.4  These  observations  strongly  suggest  the 
existence  of  (N2)0  and  (NC)0  in  these  samples.  The  calcu¬ 
lated  frequencies  for  (N2)q+  and  (NC)q+,  which  are  the  most 
stable  form  for  these  complexes  in  /7-type  ZnO,  are  2108  and 
1995  cm-1,  respectively.  These  frequencies  fall  within 
200  cm-1  of  the  experimentally  observed  ones."  Also,  the 
calculated  frequency  of  1894  cm'1  for  (CO)q+  is  consistent 


with  experiment  but  the  calculated  formation  energy  is 
higher  than  (NC)0-  Although  we  do  not  have  a  definitive 
answer  why  (CO)Q  should  coexist  with  (NC)0,  we  realize 
that  the  probability  for  an  interstitial  C;  to  find  an  O  site 
should  be  significantly  larger  than  to  find  an  N0  site. 

In  conclusion,  first  principles  total  energy  calculations 
show  that  first-row  impurities,  C,  N,  and  O,  can  form  substi¬ 
tutional  diatomic  molecules  on  the  oxygen  site  in  ZnO.  Both 
the  electronic  and  structural  properties  of  the  SDMs  are  con¬ 
sistent  with  those  of  free  molecules  but  in  contrast  to  those  of 
conventional  split  interstitials.  The  calculated  vibrational  fre¬ 
quencies  are  in  qualitative  agreement  with  recent  IR  experi¬ 
ment.  Our  calculations  further  show  that  all  these  SDMs  are 
donors  in  p-type  ZnO.  In  particular,  (NC)o  and  (N2)Q  could 
have  low  formation  energies.  These  results  should  shed  new 
lights  on  how  to  improve  the  fabrication  of  high  quality 
p- type  ZnO  by  nitrogen. 
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The  binding  sites  and  vibrational  frequencies  w  of  H  in  ZnO  are  studied  by  first-principles 
total-energy  calculations.  In  the  past,  different  experiments  have  observed  different  primary  H 
vibrational  modes,  making  the  comparison  with  theory,  and  hence  the  identification  of  the  most 
favorable  H  site,  difficult.  Here,  we  show  that  by  applying  a  hydrostatic  pressure,  one  should  be  able 
to  make  an  unambiguous  distinction,  in  particular,  between  the  bond  center  sites  and  antibonding 
sites.  This  is  because  a>  should  increase  with  pressure  for  the  former  but  decrease  for  the  latter  with 
the  magnitude  of  calculated  slopes  about  4  cm-  /GPa,  which  should  be  large  enough  to  measure. 
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Recently,  ZnO  has  attracted  a  lot  of  attention  due  to  its 
potential  for  blue  light  optoelectronics.  If  reasonably  good 
p-type  conductivity  can  be  achieved,  ZnO  could  eventually 
emerge  as  a  superior  alternative  to  GaN.  Hydrogen  is  one  of 
the  most  common  impurities  in  semiconductors.  It  is  abun¬ 
dant  in  air,  used  in  the  various  growth  techniques,  and  gen¬ 
erally  has  a  rather  high  diffusivity  inside  semiconductors. 
Usually,  H  is  an  amphoteric  impurity,  namely,  acting  either 
as  a  donor  or  an  acceptor  depending  on  the  Fermi  level. 
However,  ZnO  is  an  exception1  for  which  isolated  interstitial 
is  predicted  to  be  the  predominant  form  of  H  (Ref.  2),  and  H 
acts  only  as  a  shallow  donor. 1-5  Both  first-principles  calcula¬ 
tions  and  local  vibration  measurements  indicate  that  hydro¬ 
gen  prefers  to  bond  with  oxygen.  However,  what  is  the  most 
favorable  geometry  is  still  under  intensive  debate. “  8 

In  a  tetrahedrally  coordinated  semiconductor,  a  hydrogen 
atom  can  bind  to  an  anion  atom  predominantly  in  two  ways: 
in  the  bond  center  (BC)  site  or  in  the  antibonding  (AB)  site. 
Due  to  the  low  symmetry  of  the  ZnO  wurtzite  crystal,  how¬ 
ever,  there  are  two  types  of  orientations  (II  and  _L)  for  each 
site  above,  as  shown  in  Fig.  1(a).  First-principles  calculations 
suggested  that  all  four  sites  have  relatively  low  energies  with 
the  global-minimum  BCy  site  lower  than  the  AB0  sites  by 
less  than  0.2  eV.  The  occupation  of  the  BC  sites  by  an  H 
atom  requires  large  displacements  of  the  host  atoms,  which 
may  rarely  happen  due  to  their  relatively  large  masses  with 
respect  to  that  of  H.  Therefore,  from  a  kinetic  point  of  view, 
the  H  atom  might  occupy  the  AB0  sites,  instead,  despite  their 
slightly  higher  energies.  In  cases  like  this,  often,  a  calcula¬ 
tion  of  the  local  vibrational  modes  (LVMs)  could  differenti¬ 
ate  among  the  structures.  Indeed,  the  calculated  H  stretch 
frequencies  for  the  BC  and  AB0  sites  differ  by  about 
300  cm-1.  However,  the  measured  primary  frequencies  for 
interstitial  H  by  different  groups  also  differ  by  this  amount. 
This  makes  the  assignment  difficult  because  assigning  the 
modes  to  different  sites  would  contradict  energetic  argument 
in  favor  of  kinetic  argument.  Therefore,  it  is  highly  desirable 
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to  have  another  experimental  means  to  independently  con¬ 
firm  the  assignment. 

In  this  letter,  we  predict,  based  on  first-principles  calcu¬ 
lations,  that  an  experiment  on  the  LVMs  (e.g.,  IR,  Raman,  or 
inelastic  neutron  scattering)  under  pressure  should  indepen¬ 
dently  provide  the  critical  test  of  the  H  binding  sites,  differ¬ 
entiating  between  BC  and  AB0.  The  reason  is  because  the 
O-H  bond  length  for  the  BC  sites  decreases  under  pressure, 
resulting  in  an  increase  in  the  LVM  frequencies.  In  contrast, 
the  O-H  bond  length  for  the  AB0  sites  increases  under  pres¬ 
sure,  resulting  in  a  decrease  in  the  LVM  frequencies.  Fre¬ 
quency  shifts  (Aw),  as  much  as  ±35  cm-1,  are  predicted  over 
a  pressure  range  of  8  GPa,  at  which  phase  transition  takes 
place.9  These  shifts  should  be  large  enough  to  be  readily 
detected  by  experiments. 

We  used  the  density-functional  theory,  within  the  local- 
density  approximation  (LDA),  and  the  Vanderbilt- type  ultra- 
soft  pseudopotentials,10  as  implemented  in  the  VASP  code.11 
Zinc  3d  states  are  treated  as  valence  states.  The  cutoff  energy 
for  the  plane-wave  expansion  is  400  eV,  with  additional  tests 
carried  out  at  515  eV.  The  calculated  equilibrium  lattice  con¬ 
stant  of  ZnO  is  only  1 .4%  shorter  than  experiment.  We  used 
a  supercell  approach  with  a  Monkhorst-Pack  A: -point  mesh 


FIG.  1.  (a)  Schematic  representation  of  the  four  H  sites  in  ZnO  and  (b) 
calculated  atomic  positions  for  H+  at  the  BCy  site,  in  the  (11-20)  plane.  The 
dashed  circles  and  dashed  lines  indicate  bulk  atomic  positions  and  bonds, 
respectively.  Distances  are  given  in  unit  of  A  and  percentage  changes  are 
with  respect  to  the  bond  length  of  bulk  ZnO. 
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TABLE  I.  Calculated  formation  energy  with  respect  to  that  of  BC||  and  the 
O-H  bond  length  for  H+  in  ZnO. 


Site 

A E  (eV) 

^O-H 

(A) 

Present 

Ref.  2 

Ref.  6 

Present 

Ref.  6 

BC, 

0.00 

0.00 

0.0 

0.986 

0.990 

bca 

0.15 

-0.02 

0.1 

0.982 

ABqii 

0.17 

0.23 

0.1 

1.003 

1.006 

AB01 

0.14 

0.04 

0.1 

1.003 

for  Brillouin-zone  integration  (2X2X2  for  the  36-atom  cell 
and  1  X  1  X  1  for  the  96-atom  cell).  For  charged  H,  we  use  a 
jellium  background.  All  atoms  in  the  supercell  are  allowed  to 
relax  in  the  structural  optimizations.  A  test  using  the  96-atom 
cell  shows  that  the  energy  difference  between  the  two  BC 
configurations  in  the  36-atom  calculations  are  converged  to 
within  0.01  eV.  The  smaller  supercell  is  used  for  extracting 
the  energy  difference  needed  for  LVM  calculation. 

We  simulate  the  stretch  frequency  by  tracing  the  poten¬ 
tial  energy  curve  of  the  vibration.  This  has  been  carried  out 
by  calculating  the  energy  as  a  function  of  H  displacement  in 
both  compression  (toward  O)  and  extension  (away  from  O) 
directions.  Because  (1)  H  is  much  lighter  than  the  host  atoms 
and  (2)  the  LVM  of  the  O-H  bond  (around  3000  cm-1)  is 
much  higher  than  the  crystal  phonon  modes,  the  displace¬ 
ment  of  the  host  atoms  during  the  vibration  can  thus  be  ne- 
glected.  A  previous  study  “  has  explicitly  compared  the 
LVMs  obtained  by  freezing  host  atoms  and  the  LVMs  ob¬ 
tained  from  dynamic  matrix  approach.  A  good  agreement 
was  found  provided  that  the  reduced  mass  of  the  anion-H 
pair  is  used.  Due  to  the  light  mass  of  the  H  atom,  anharmonic 
contribution  ( wAC)  could  also  be  important.13  We  include 
«ac  in  our  calculation  by  fitting  the  calculated  energy  curve 
as  a  function  of  the  H  displacement  with  a  fourth-degree 
polynomial,  i.e., 

k 

V{x)=~x2+ax3  +  /3x4.  (1) 

To  solve  the  Schrodinger  equation  with  such  a  potential,  we 
used  a  perturbation  theory14  to  obtain  an  approximate  ana¬ 
lytical  solution,  which  reads 

—  _  A/1  ^/5a2  f}\ 

w-wHar+wAC-  \  “3  2  -  I.  (2) 

y  /r  fL  \  Zk~  k  / 


Table  I  shows  the  calculated  energy  for  the  four  H  sites 
and  the  corresponding  O-H  bond  lengths.  It  suggests  that 
BCy  is  the  most  stable  H  configuration,  in  agreement  with  the 
recent  results  by  Lavrov  et  al.6  Figure  1(b)  shows  the  atomic 
geometry  of  the  fully  relaxed  BC||.  We  see  that,  with  respect 
to  ideal  ZnO,  the  Zn  atom  adjacent  to  the  H  atom  has  relaxed 
away  from  the  H  through  the  basal  plain  defined  by  the  three 
remaining  oxygen  atoms  by  as  much  as  41.1%.  As  a  result, 
the  three  remaining  Zn-O  bonds  are  compressed  by  3.3%. 
Our  calculated  O-H  bondlengths  (c/0-h)  agree  with  those  in 
Ref.  6  to  within  0.005  A. 

Table  II  shows  the  calculated  LVM  frequencies, 
wHar,  coAC,  and  <5,  as  defined  in  Eq.  (2).  The  anharmonic 
contribution  wAC,  for  example,  -282  cm-1  for  BCy  and 
-352  cm-1  for  AB0y,  is  typically  8-10  %  of  the  correspond¬ 
ing  harmonic  contribution  <wHar,  3659,  and  3370  cm-1,  re¬ 
spectively,  and  the  two  have  opposite  signs.  It  appears  that 
the  difference  in  u>  between  the  II  and  ±  modes  are  relatively 
small,  typically  40  cm-1,  with  the  latter  being  slightly  larger. 
We  may  expect  some  errors  in  our  calculated  Li’s,  as  u>  for 
the  free  water  molecule  in  Table  II  is  lower  than 
experiment15  by  ojhr=98  cm-1.  The  error  is  most  likely  sys¬ 
tematic  because  the  calculated  O-H  bond  lengths  are  consis¬ 
tently  longer  than  experiment  by  about  2%.  We  have  there¬ 
fore  corrected  the  errors  for  all  the  w’s  by  adding  ojkr  to 
obtain  oj.  The  same  procedure,  but  with  ojkr  =  5  1  cm-1,  is 
applied  to  the  calculated  results  in  Ref.  6.  After  the  correc¬ 
tion,  the  average  values  for  the  BC  and  AB0  modes  are 
oj(BC)  =  3497  and  w(AB0)  =  3135  cm-1,  respectively.  These 
can  be  compared  with  the  corrected  values  in  Table  II  for 
Ref.  6,  w(BC)  =  3551  and  «(AB0)  =  3259  cm'1.  Thus,  theo¬ 
ries  agree  to  within  3551-3497  =  54  cm-1  for  <w(BC)  and 
3259-3135=  124  cm-1  for  oj(AB0).  These  differences,  54 
and  124  cm-1,  are  significantly  smaller  than  the  calculated 
splitting  between  the  BC  and  AB0  modes:  362  (present)  and 
292  cm-1  (Ref.  6).  It  comes  to  our  attention  that  the  mea¬ 
sured  LVM  frequencies  in  Refs.  6,8  differ  by  285  cm-1,  sug¬ 
gesting  that  one  of  them  could  be  the  BC  mode  whereas  the 
other  is  the  AB0  mode. 

More  definitive  identification  of  the  experimental  LVMs 
may  be  obtained  by  pressure  measurement.  This  is  possible 
because  the  neighboring  atoms  of  the  BC  and  AB0  sites  are 
different,  which  could  result  in  different  responses  to  exter¬ 
nal  pressure.  To  simulate  the  pressure  effect,  we  need  the 
total  energy  as  a  function  of  volume,  from  which  the 


TABLE  II.  Stretch  modes  for  interstitial  H+  in  ZnO,  in  units  of  cm  1 .  cuHar,  <oAC,  and  m  are  described  in  Eq.  (2), 
whereas  the  calculated  w’s  (columns  4  and  5)  include  the  systematic  error  correction,  as  described  in  the  text. 


Site 

"Har 

Present 

“AC 

co(a>) 

Ref.  6 
a>(a>) 

Measured 

CO 

BC, 

3659 

-282 

3475  (3377) 

3551  (3500) 

361  la 

BCl 

3706 

-285 

3519  (3421) 

AB0|| 

3370 

-352 

3116  (3018) 

3259  (3208) 

3326b 

AB01 

3379 

-323 

3154  (3056) 

H20 

3815 

-257 

3657  (3559) 

3657  (3606d) 

3657c 

JReference  6. 
h Reference  S. 
Reference  15 


^Calculated  using  FHI98  codes. 
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FIG.  2.  The  shifts  in  the  stretch  mode  frequencies  of  the  O-H  bonds  for  the 
four  H  sites  in  Fig.  1,  as  a  function  of  the  external  pressure.  Two  possible 
orientations  of  the  O-H  bonds  in  the  crystal  are  coded  with  solid  (II)  and 
open  (X)  symbols.  The  lines  are  least-square  fits  to  the  calculated  data. 

pressure  P{V)  as  a  function  of  V  is  obtained  through 
P(V)=-dE(V)/dV.  Next,  we  need  the  frequency  as  a  func¬ 
tion  of  volume,  oj( V) .  From  co( V)  and  P( V),  we  can  deduce 
( o(P ).  The  above  calculation  should  be  straightforward  if  we 
use  very  large  supercells.  In  practice,  we  use  a  finite-size 
36-atom  supercell,  which  might  introduce  errors  both  in 
P{V)  and  oj( V) .  To  minimize  the  error  in  P{V),  we  calculate 
instead  dE(V)/dV  using  a  primitive  cell  for  the  reasons  that 
P  is  an  intensive  quantity  and  that  in  the  dilute  limit  the 
effect  of  the  defect  should  diminish.  To  minimize  the  error  in 
co(V),  we  did  selectively  96-atom  cell  calculations  to  adjust 
the  36-atom  cell  results.  The  correction  on  co  was  found  to  be 
rather  small,  whose  main  role  was  to  eliminate  the  artificial 
small  splitting  between  the  II  and  X  modes  due  to  the  specific 
shape  of  the  36-atom  cell.  Figure  2  shows  the  calculated  co 
for  the  four  H  binding  sites,  as  a  function  of  P.  We  see 
clearly  the  qualitative  difference  between  the  LVMs  for  the 
BC  and  AB0  sites:  while  oj(BC)  increases  with  pressure, 
w(AB0)  decreases,  instead.  In  spite  of  that,  w(BC)  and 
w(AB0)  have,  fortuitously,  a  similar  magnitude  of  slope 
around  4.2  ctrr'/GPa. 

The  distinct  pressure  dependences  in  Fig.  2  may  be 
qualitatively  understood  by  the  change  of  the  O-H  bond 
length  (c/0-h)  under  pressure,  namely,  a  decrease  in  af0_H 
would  correspond  to  an  increase  in  co,  as  demonstrated  in 
Ref.  13.  Here,  dQ_ H  for  AB0  monotonically  increases  with  P, 
which  explains  why  co  decreases  in  Fig.  2,  e.g.,  at  P 
=  5.4  GPa,Aii0.H  =  0.0013  A  whereas  Aw=-23  cm-1  for 
AB0||.  On  the  other  hand,  d0_H  for  BC  monotonically  de¬ 
creases  with  P,  which  explains  why  co  increases  in  Fig.  2:  at 
P= 5.2  GPa,A<io-H= -0.0021  A  whereas  Aw=+22  cm-1  for 
BC||.  We  speculate  that  the  change  in  d0_H  for  AB0  under 
pressure  is  because  of  the  compression  of  the  neighboring 
Zn-O  bonds.  To  test  this  hypothesis,  we  have  calculated  an 
acetylene  molecule  (H-C-C-H)  and  found  that  indeed  if  we 
compress  the  C-C  bond  by  0.02  A,  both  of  the  C-H  bonds 
will  increase,  as  expected,  by  about  0.001  A.  Concerning  the 
change  in  c/0_H  for  BC  under  pressure,  we  note  that  the  H 
atom  sits  between  an  O  and  a  Zn  atom.  When  the  O-Zn 
distance  is  reduced  due  to  the  external  pressure,  the  O-H 
distance  (c/0_H)  is  also  reduced. 

A  recent  experiment  showed  that  the  3326  cm-1  fre¬ 
quency  corresponds  to  an  H  mode  in  the  X  direction,  whose 
frequency  changes  with  pressure  by  a  negative  slope. 


-1.5  cnT'/GPa  (Ref.  16).  This  mode  has  also  been  assigned 
to  hydrogen  on  the  AB0 ,  sites.  In  fact,  this  assignment 
agrees  with  our  present  calculation  where  the  AB01  site 
should  have  a  negative  pressure  slope,  as  well  as  a  slightly 
lower  formation  energy  than  that  of  the  AB0n  site  (see  Table 
I).  Even  the  measured  slope,  -1.5  cm_1/GPa,  is  within  a 
factor  of  3  of  the  calculated  one,  -4.2  cm_1/GPa.  Note  that 
to  calculate  the  small  changes  in  co  accurate  to  within 
10  cm-1  requires  a  very  high  convergence  to  at  least  within 
0.001  A  in  the  atomic  positions  and  1  meV  in  energy.  Hence, 
such  a  deviation  from  experiment  is  not  at  all  surprising. 
Also,  while  both  the  measured  3326  and  3611  cm-1  modes  in 
Table  II  have  been  previously  assigned  to  interstitial  H,  there 
is  always  the  concern  whether  they  belong  to  some 
hydrogen/(unknown)-defect  complexes,  instead.  We  believe 
that  the  pressure  dependence  of  the  LVMs  should  not  depend 
very  sensitively  on  the  complex  formation.  In  other  words, 
the  pressure  coefficient  should  still  be  positive  if  the  H  in  a 
complex  can  be  described  as  a  bond-center  H,  but  negative  if 
the  H  can  be  described  as  an  antibonding  H. 

In  conclusion,  first-principles  studies  predict  that 
hydrostatic -pressure  experiments  could  be  an  important  tool 
for  making  unambiguous  identification  of  interstitial  H  in 
ZnO;  in  particular,  showing  a  clear,  qualitative  difference 
between  the  bond  center  and  antibonding  sites.  The  fre¬ 
quency  co  for  the  former  should  increase  with  pressure,  but  co 
for  the  latter  should  decrease,  with  a  slope  of  4.2  cnT'/GPa. 
The  magnitude  of  the  calculated  slopes  qualitatively  agrees 
with  the  only  known  experiment. 

We  thank  M.  D.  McCluskey  and  S.-H.  Wei  for  helpful 
discussions.  This  work  was  supported  by  the  Thailand  Re¬ 
search  Fund  under  Contract  No.  BRG4680003  and  by  the 
AFOSR/AOARD  under  Contract  No.  FA5209-05-P-0309. 
The  work  in  the  United  States  was  supported  by  the  U.S. 
DOE/BES  and  DOE/EERE  under  Contract  No.  DE-AC36- 
99GO10337  and  by  the  NERSC  for  MPP  time. 

'C.  G.  Van  de  Walle  and  J.  Neugebauer,  Nature  (London)  423,  626  (2003). 
2C.  G.  Van  de  Walle.  Phys.  Rev.  Lett.  85,  1012  (2000). 

3S.  F.  J.  Cox,  E.  A.  Davis,  S.  P.  Cottrell,  P.  J.  C.  King,  J.  S.  Lord,  J.  M.  Gil, 
H.  V.  Alberto,  R.  C.  Vilao,  J.  Piroto  Duarte,  N.  Ayres  de  Campos,  A. 
Weidinger,  R.  L.  Lichti,  and  S.  J.  C.  Irvine,  Phys.  Rev.  Lett.  86,  2601 
(2001). 

4D.  M.  Hofmann,  A.  Hofstaetter,  F.  Leiter,  H.  Zhou,  F.  Henecker,  B.  K. 
Meyer,  S.  B.  Orlinskii,  J.  Schmidt,  and  P.  G.  Baranov,  Phys.  Rev.  Lett.  88, 
045504  (2002). 

K.  Shimomura,  K.  Nishiyama,  and  R.  Kadono,  Phys.  Rev.  Lett.  89. 
255505  (2002). 

6E.  V.  Lavrov,  J.  Weber,  F.  Borrnert,  C.  G.  Van  de  Walle,  and  R.  Helbig, 
Phys.  Rev.  B  66,  165205  (2002). 

7C.  H.  Seager  and  S.  M.  Myers,  J.  Appl.  Phys.  94,  2888  (2003). 

SM.  D.  McCluskey,  S.  J.  Jokela,  K.  K.  Zhuravlev,  P.  J.  Simpson,  and  K.  G. 
Lynn,  Appl.  Phys.  Lett.  81,  3807  (2002). 

9S.  Limpijumnong  and  S.  Jungthawan,  Phys.  Rev.  B  70,  054104  (2004). 
10D.  Vanderbilt,  Phys.  Rev.  B  41,  7892  (1990). 

UG.  Kresse  and  J.  Furthmuller,  Comput.  Mater.  Sci.  6,  15  (1996). 

12S.  Limpijumnong  and  C.  G.  Van  de  Walle,  Mater.  Res.  Soc.  Symp.  Proc. 
639,  G4.3  (2001). 

13S.  Limpijumnong,  J.  E.  Northrup,  and  C.  G.  Van  de  Walle.  Phys.  Rev.  B 
68,  075206  (2003). 

14L.  D.  Landau  and  E.  M.  Lifshitz,  Quantum  Mechanics ,  3rd  ed.  (Pergamon, 
Oxford,  1977),  p.  136. 

15E.  McCartney,  Absorption  and  Emission  by  Atmospheric  Gases:  The 
Physical  Processes  (Wiley,  New  York,  1983). 

16M.  D.  McCluskey  and  S.  J.  Jokela,  in  Proceedings  of  the  NATO  Advanced 
Research  Workshop  on  Zinc  Oxide,  edited  by  N.  Nickel  and  E.  Terukov, 
NATO  Science  Series  II  (2005),  Vol.  194. 


Downloaded  07  Apr  2005  to  203.147.33.2.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://apl.aip.org/apl/copyright.jsp 


PHYSICAL  REVIEW  B  71,  125209  (2005) 


Properties  of  Ga-interstitial  defects  in  Al^Ga^Nj,?!^ 

N.  Q.  Thinh,  I.  P.  Vorona,*  I.  A.  Buyanova,  and  W.  M.  Chen1 

Department  of  Physics  and  Measurement  Technology,  Linkoping  University,  58183  Linkoping,  Sweden 

Sukit  Limpijumnong 

School  of  Physics,  Institute  of  Science,  Suranaree  University  of  Technology,  Nakhon  Ratchasima  30000,  Thailand 

S.  B.  Zhang 

National  Renewable  Energy  Laboratory,  Golden,  Colorado  80401,  USA 

Y.  G.  Hong,  H.  P.  Xin,  and  C.  W.  Tu 

Department  of  Electrical  and  Computer  Engineering,  University  of  California,  La  Jolla,  California  92093-0407,  USA 

A.  Utsumi,  Y.  Furukawa,  S.  Moon,  A.  Wakahara,  and  H.  Yonezu 
Department  of  Electrical  and  Electronic  Engineering,  Toyohashi  University  of  Technology,  Toyohashi,  Aichi,  441-8580,  Japan 

(Received  14  October  2004;  published  18  March  2005) 


A  detailed  account  of  the  experimental  results  from  optically  detected  magnetic  resonance  (ODMR)  studies 
of  grown-in  defects  in  (Al)GaNP  alloys,  prepared  by  molecular  beam  epitaxy,  is  presented.  The  experimental 
procedure  and  an  in-depth  analysis  by  a  spin  Hamiltonian  lead  to  the  identification  of  two  Ga;  defects 
(Ga,-A  and  Ga ,-B).  New  information  on  the  electronic  properties  of  these  defects  and  the  recombination 
processes  leading  to  the  observation  of  the  ODMR  signals  will  be  provided.  These  defects  are  deep-level 
defects.  In  conditions  when  the  defect  is  directly  involved  in  radiative  recombination  of  the  near-infrared 
photoluminescence  band,  the  energy  level  of  the  Ga rB  defect  was  estimated  to  be  deeper  than  ~1.2  eV  from 
either  the  conduction  or  valence  band  edge.  In  most  cases,  however,  these  defects  act  as  nonradiative  recom¬ 
bination  centers,  reducing  the  efficiency  of  light  emission  from  the  alloys.  They  can  thus  undermine  the 
performance  of  potential  photonic  devices.  High  thermal  stability  is  observed  for  these  defects. 
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I.  INTRODUCTION 

One  of  the  most  important  aspects  in  the  physics  and  tech¬ 
nology  of  semiconductors  is  to  understand  and  control  de¬ 
fects,  which  are  known  to  often  produce  deep  levels  in  the 
forbidden  energy  gap.  Intrinsic  defects  such  as  self¬ 
interstitials,  vacancies,  and  antisites  have  been  found  to  be 
commonly  occurring  defects  in  semiconductor  crystals 
grown  under  nonoptimized,  nonequilibrium  conditions,  or 
subject  to  high-energy  particle  bombardment.  They  are  fun¬ 
damental  building  blocks  for  defect  complexes  and  in  some 
cases  nucleation  sites  for  the  formation  of  clusters  and  ex¬ 
tended  defects.  Intrinsic  defects  and  their  complexes  are 
known  to  play  a  crucial  role  in  determining  the  electronic 
and  optical  properties  of  semiconductors.  In  some  cases 
when  their  properties  are  well  understood  and  controlled, 
they  can  be  used  to  tailor  the  physical  properties  of  semicon¬ 
ductors  such  as  carrier  lifetime,  Fermi  level  position,  and 
doping  [e.g.,  the  antisite-related  defects  in  GaAs  (Refs.  1  and 
2)  and  PIn  antisite  in  InP  (Refs.  3  and  4)].  In  most  cases, 
however,  defects  are  known  to  severely  degrade  device  per¬ 
formance  and  can  hinder  a  semiconductor  from  practical  ap¬ 
plications.  In  any  case,  a  good  understanding  of  the  defect 
properties  is  vital,  especially  concerning  a  positive  identifi¬ 
cation  of  their  chemical  nature  and  their  role  in  altering  ma¬ 
terial  properties.  Up  to  now  only  very  few  experimental  tech¬ 
niques  have  been  proven  successful  in  meeting  the  demand. 
Among  them,  magnetic  resonance  techniques  have  been 


shown  to  be  the  most  powerful  of  the  ones  that  have  success¬ 
fully  identified  many  intrinsic  defects  in  semiconductors.5"7 
When  the  conventional  electron  spin  resonance  (ESR)  tech¬ 
nique  fails  to  provide  the  required  sensitivity  to  reveal  de¬ 
fects  in  semiconductor  thin  films  and  quantum  structures 
grown  by  modern  epitaxy  techniques,  its  advanced  variant — 
namely,  the  optically  detected  magnetic  resonance  (ODMR) 
technique — has  been  shown  to  be  successful.8"11 

The  aim  of  the  present  work  is  to  reveal  important 
grown-in  defects  in  (Al)GaNP  alloys  by  employing  the 
ODMR  technique.  (Al)GaNP  represents  a  class  of  newly 
emerging  dilute  nitrides,  which  have  recently  attracted  great 
attention.  These  materials  exhibit  unusual  and  fascinating 
new  physical  properties,  such  as  a  giant  band-gap  bowing 
that  allows  widely  extended  band-structure  engineering.12"19 
One  can  also  vary  the  lattice  constant  of  the  alloy  materials 
over  a  wide  range  by  varying  the  N  content.  For  example, 
adding  about  1.8%  of  N  into  GaP  can  make  the  parent  ma¬ 
terial  lattice  matched  to  Si,  opening  the  possibility  for  opto¬ 
electronic  integrated  circuits19,20  and  multijunction  III-V  so¬ 
lar  cells  fabricated  on  Si  wafers.  A  crossover  from  an  indirect 
to  a  direct  band  gap  induced  by  incorporation  of  >0.4%  N 
(Refs.  15-18)  promises  high  radiative  efficiency.  Due  to  the 
large  mismatch  in  size  and  electronegativity  between  N 
and  the  replaced  P  atom,  however,  successful  growth  of 
(Al)GaNP  has  to  be  undertaken  under  nonequilibrium 
conditions.21,22  The  nonequilibrium  growth  conditions,  to¬ 
gether  with  the  disparity  between  N  and  P  atoms  resulting  in 
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a  strong  local  potential  and  strain  field,  lead  to  the  formation 
of  various  point  defects  and  clusters.  The  efficiency  of  vis¬ 
ible  light  emissions  in  (Al)GaNP  has  been  shown  to  deterio¬ 
rate  with  increasing  N  and  Al  compositions.  Up  to  now, 
band-to-band  transitions  have  not  been  observed  in  this  ma¬ 
terial  despite  the  crossover  to  a  direct  band  gap,  due  pre¬ 
dominantly  to  nonradiative  recombination  at  the  defects. 
Currently,  there  is  a  great  need  in  identifying  important 
grown-in  defects  in  this  dilute  nitride  and  in  assessing  their 
exact  roles  in  carrier  recombination. 

In  this  work,  we  were  able  to  identify  Ga,  interstitial  de¬ 
fects  by  ODMR  as  the  dominant  grown-in  intrinsic  defects  in 
(Al)GaNP,  prepared  by  molecular  beam  epitaxy  (MBE).  We 
obtained  detailed  information  on  the  electronic  properties  of 
the  defects  and  their  role  in  carrier  recombination.  Prelimi¬ 
nary  results  were  briefly  presented  in  a  recent  Rapid 
Communication.23  In  this  paper,  we  shall  present  in  details 
our  experimental  procedure  and  an  in-depth  analysis  by  a 
spin  Hamiltonian  that  led  to  the  identification  of  the  two  Ga, 
defects.  New  information  on  the  electronic  properties  of 
these  defects  and  the  recombination  processes  leading  to  the 
observation  of  the  ODMR  signals  will  be  provided.  A  discus¬ 
sion  about  the  formation  mechanism  of  the  defects  will  not 
be  included  here,  but  can  be  found  elsewhere.24  The  paper  is 
organized  as  follows:  In  Sec.  II,  we  shall  describe  the 
samples  under  study,  the  details  of  the  ODMR,  and  related 
photoluminescence  (PL)  experiments.  In  Sec.  Ill,  we  shall 
present  a  detailed  account  of  the  ODMR  results  including 
temperature,  angular,  and  spectral  dependences.  An  analysis 
of  the  ODMR  results  by  a  spin  Hamiltonian  and  a  discussion 
of  the  microscopic  structures  of  the  defects  and  the  involved 
recombination  processes  will  be  presented  in  Sec.  IV.  The 
main  conclusions  of  the  present  study  will  be  summarized  in 
Sec.  V. 

II.  EXPERIMENT 
A.  Samples 

Three  types  of  (Al)GaNP  samples  were  studied  in  this 
work,  varying  by  their  chemical  compositions  and  substrate 
materials.  The  first  two  types  were  grown  at  520  °C  on  GaP 
substrates  by  gas-source  molecular  beam  epitaxy  (GS-MBE) 
and  have  a  typical  thickness  between  0.75  and  0.9  /um.  They 
are  (1)  quaternary  Alj-Ga^NyP^.,  with  a  fixed  value 
v=0.012  and  varied  x  values,  x=0.00,  0.01,  0.02,  and  0.30; 
(2)  Al-free  GaNvP!_v  with  varied  y  values,  y=0.013,  0.023, 
and  0.031.  They  will  be  referred  below  as  AlGaNP/GaP  and 
GaNP/GaP,  respectively.  The  third  type  (denoted  as 
GaNP/Si  below)  is  GaN00i8P0982  with  a  thickness  of 
0.15  /jl m  grown  at  590  °C  on  a  Si  (100)  substrate,  by  solid- 
source  MBE.  Rapid  thermal  annealing  (RTA)  was  performed 
for  30  sec  at  the  temperature  range  700-900  °C  with  halo¬ 
gen  lamps  in  a  flowing  N2  ambient. 

B.  Experimental  procedure 

The  ODMR  measurements  were  performed  at  two  micro- 
wave  frequencies  on  a  modified  Bruker  ER-200  D  X-band 
(—9.3  GHz)  setup  with  a  flow  cryostat  and  a  R'-band 


Photon  Energy  (eV)  Magnetic  Field  (ml) 

FIG.  1.  Representative  PL  and  ODMR  spectra  from  three  types 
of  (Al)GaNP  samples  studied  in  this  work.  The  PL  spectra  were 
measured  at  2  K.  The  ODMR  spectra  were  measured  at  5  K  and  at 
the  W  band  (94.8  GHz),  by  monitoring  the  NIR  PL  emissions.  The 
sign  of  the  ODMR  signals  is  negative  for  GaNP/GaP  and 
AlGaNP/GaP,  but  positive  for  GaNP/Si.  They  are  all  shown  as 
positive  for  easy  viewing. 

(—95  GHz)  setup  with  an  Oxford  superconducting  split-coil 
magnet  (0-5  T)  cryostat.  A  solid-state  laser  with  532-nm 
wavelength  was  used  as  an  excitation  source  in  the  ODMR 
measurements.  The  PL  signal  was  detected  through  proper 
optical  filters  by  a  cooled  Ge  detector  and  a  GaAs  photodi¬ 
ode  for  the  near-infrared  (NIR)  and  visible  spectral  range, 
respectively.  The  ODMR  signals  were  detected  by  the 
lock-in  technique  as  spin-resonance-induced  changes  of  PL 
intensity,  in  phase  with  an  amplitude  modulation 
(0.5-20  kHz)  of  the  microwave  radiation.  The  typical  mi¬ 
crowave  power  employed  was  170-200  mW.  The  measure¬ 
ments  were  performed  in  the  temperature  range  of  2-40  K. 

III.  RESULTS 
A.  PL  and  ODMR 

Upon  above-band-gap  optical  excitation,  all  three  types  of 
the  samples  studied  exhibit  PL  emissions  in  the  visible  spec¬ 
tral  range  (see  Lig.  1)  that  are  known  to  arise  from  N-related 
localized  states.25  They  also  give  rise  to  defect-related  deep 
PL  emissions  within  the  NIR  spectral  range  whose  origins 
are  still  largely  unknown  (see  the  PL  bands  at  the  lower 
energies  in  Lig.  1).  Judging  from  their  energy  positions  and 
their  shifts  with  N  and  Al  compositions,  as  well  as  different 
signs  of  the  ODMR  signals  (to  be  presented  below),  we  be¬ 
lieve  that  the  deep  emissions  from  the  three  types  of  samples 
are  of  different  origins.  Representative  PL  spectra  from  the 
three  types  of  samples  are  given  in  Lig.  1,  together  with 
typical  ODMR  spectra  obtained  by  monitoring  these  PL 
emissions.  The  ODMR  spectra  show  the  following  two  dis¬ 
tinct  features  originated  from  different  defects.  The  first  one 
is  a  single  strong  line  situated  in  the  middle  of  the  ODMR 
spectra  from  the  AlGaNP/GaP  and  GaNP/GaP  samples, 
with  a  g  value  close  to  2.  Due  to  a  lack  of  hyperfine  (HP) 
structure,  unfortunately,  the  chemical  nature  of  the  corre- 
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Magnetic  Field  (mT) 

FIG.  2.  ODMR  spectra  taken  at  two  different  temperatures, 
5  and  40  K,  from  GaN003iPo.969/GaP  and 
Al0.02Ga0.98N0.0i8P0.982/CiaP.  The  Y-band  microwave  frequency 
was  9.31  GHz. 
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FIG.  3.  ODMR  spectra  of  the  three  types  of  (Al)GaNP  samples 
obtained  at  5  K  and  when  the  external  magnetic  field  was 
oriented  along  the  three  main  crystallographic  axes.  The 
Al0.02Ga098N0.0i8P0.982/GaP  sample  was  measured  at  the  Y  band, 
while  the  other  two  samples  were  measured  at  the  W  band. 


sponding  defect  cannot  be  identified.  Below  we  shall  simply 
refer  to  it  as  the  “g  =  2”  defect.  The  second  feature  of  the 
ODMR  spectra  consists  of  a  complicated  pattern  of  lines 
spreading  over  a  wide  field  range.  Such  multiple  ODMR 
lines  can  arise  from  a  high-electron-spin  state  exhibiting  a 
zero-field  splitting  due  to  a  defect  crystal  field,  from  an  HF 
interaction  involving  a  high  nuclear  spin,  or  from  several 
overlapping  ODMR  signals  due  to  different  defects.  To  re¬ 
solve  this  issue,  the  temperature  and  angular  dependences  of 
the  ODMR  spectra  were  investigated. 

B.  Temperature  dependence  of  the  ODMR 

To  examine  if  the  second  feature  of  the  ODMR  spectra 
originates  from  a  single  or  more  than  one  defect,  temperature 
dependence  studies  were  performed.  This  is  because  the  tem¬ 
perature  dependence  of  spin-lattice  relaxation  is  often  differ¬ 
ent,  such  that  relative  intensity  of  the  ODMR  signals  from 
different  defects  may  be  altered  by  the  measurement  tem¬ 
perature.  ODMR  spectra  from  AlGaNP/GaP  and 
GaNP/GaP  at  5  K  and  40  K  (Fig.  2)  clearly  show  that  about 
half  of  the  multiple  ODMR  lines  undergo  thermal  quenching 
in  intensity.  This  finding  provides  strong  evidence  that  the 
multiline  feature  of  the  ODMR  spectra  in  fact  consists  of  two 
groups  of  lines  from  two  different  defects  (denoted  as 
Gar,4  and  Ga -B),  of  which  Ga,-/i  loses  signal  strength  at 
elevated  temperatures  and  vanishes  at  T  3*40  K  (see  Fig.  2). 

C.  Angular  dependence  of  the  ODMR 

Both  Ga,-A  and  Ga -B  give  rise  to  multiple  ODMR  lines, 
which  can  arise  either  from  a  high-electron-spin  state  exhib¬ 
iting  a  zero-field  splitting  caused  by  a  defect  crystal  field  or 
from  an  HF  interaction  involving  a  high  nuclear  spin.  The 
former  should  give  rise  to  an  anisotropic  ODMR  spectrum 
reflecting  the  symmetry  of  the  defect,  whereas  the  latter  may 
or  may  not  do  so.  To  identify  the  origin  of  the  multiple 
ODMR  lines,  angular  dependence  studies  were  carried  out 
by  rotating  the  external  magnetic  field  in  a  {110}  crystallo¬ 
graphic  plane.  Typical  ODMR  spectra  are  displayed  in  Fig.  3 
when  the  direction  of  the  external  magnetic  field  is  oriented 
along  the  three  main  crystallographic  axes,  as  examples.  The 
ODMR  spectra  are  shown  to  be  isotropic,  within  the  experi¬ 


mental  error.  This  finding  supports  the  assignment  that  the 
multiple  ODMR  lines  for  each  defect  should  be  due  to  an  HF 
structure,  which  will  be  shown  below  to  arise  from  a  strong 
interaction  between  an  unpaired  electron  spin  (5=1/2)  and 
the  nuclear  spin  (7=3/2)  of  a  Ga  atom. 

D.  Spectral  dependence  of  the  ODMR 

In  the  AlGaNP/GaP  and  GaNP/GaP  samples,  all 
defects — i.e.,  the  “g  =  2”  defect,  the  Ga,-A,  and  the 
Ga rB — give  rise  to  negative  ODMR  signals  that  correspond 
to  a  decrease  upon  spin  resonance  in  the  intensity  of  the 
monitored  PL  emissions  over  the  entire  visible  and  NIR 
spectral  range.  (The  ODMR  signals  are  displayed  as  positive 
ones  throughout  the  paper  merely  for  easy  viewing.)  The 
only  exception  occurs  for  the  Al^Gai^NyP^  sample  series 
when  x=0.00  and  y=0.012,  where  the  sign  of  the  ODMR 
signals  turns  to  positive  by  monitoring  the  NIR  PL  emis¬ 
sions.  A  possible  reason  for  this  sign  change  will  be  dis¬ 
cussed  in  Sec.  IV  E.  In  the  GaNP/Si  samples,  however,  the 
ODMR  signals  (shown  in  Figs.  1  and  3)  correspond  to  a 
spin-resonance-induced  increase  in  the  monitored  NIR  PL 
emission — i.e.,  a  positive  sign  of  ODMR.  On  the  other  hand, 
in  the  visible  spectral  range,  the  ODMR  signal  strength  was 
too  weak  to  be  detected. 

E.  Effect  of  RTA  on  the  ODMR 

In  Fig.  4,  we  show  representative  ODMR  spectra  from  the 
AlGaNP  alloys  before  and  after  post-growth  RTA  treatments. 
All  the  defects  under  study  exhibit  high  thermal  stability  up 
to  900  °C  and  can  only  be  partially  removed  by  RTA.  It 
should  be  noted  that  the  RTA  causes  a  larger  reduction  of  the 
“g  =  2”  defect  than  the  two  Ga,  defects. 

IV.  DISCUSSION 
A.  Spin  Hamiltonian  analysis 

The  left  panels  of  Figs.  5-7  show  that  the  ODMR  spectra 
at  the  X  band  are  rather  irregular,  which  makes  it  very  diffi¬ 
cult  to  assign  them  to  a  specific  nuclear  spin  state.  For  this 
reason,  the  ODMR  was  also  performed  at  the  W  band  in  a 
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FIG.  4.  ODMR  spectra  from  the  as-grown  and  RTA  treated 
Alo.o2Gao.98No.oi8Po.982/GaP  samples,  obtained  at  5  K  for  the  X 
band.  The  RTA  treatment  was  done  at  900  °C. 

higher  magnetic  field,  which  usually  avoids  a  strong  nonlin¬ 
ear  Zeeman  effect  due  to  level  repelling  at  low  fields.  The 
ODMR  spectra  at  the  W  band  show  a  clear  pattern  with  two 
groups  of  lines  from  GarA  and  Ga rB  with  distinctly  differ¬ 
ent  intensities.  Each  defect  shows  grossly  four  lines  with 
additional  features  (a  shoulder  or  a  splitting).  A  close  inspec¬ 
tion  reveals  that  each  defect  in  fact  gives  rise  to  two  sets  of 
four  lines  with  a  fixed  intensity  ratio  of  40/60  and  a  fixed 
ratio  of  line  spacing=  1.3.  The  most  plausible  explanation  for 
these  experimental  findings  is  that  the  ODMR  lines  from 
each  defect  should  originate  from  the  HF  structure  related  to 
an  atom  with  the  following  properties:  (i)  it  has  two  isotopes, 
of  which  both  have  a  nuclear  spin  7=3/2  to  account  for  the 
two  sets  of  four  lines;  (ii)  the  ratio  of  the  natural  abundances 
of  the  two  isotopes  should  be  near  40/60  to  account  for  the 
ratio  of  the  corresponding  ODMR  intensity;  (iii)  the  ratio  of 
the  nuclear  magnetic  moments  of  the  two  isotopes  should  be 
near  1.3  to  account  for  the  ratio  of  the  line  spacing. 

To  confirm  the  assignment  and  to  obtain  information  on 
the  physical  properties  of  the  defects,  a  detailed  analysis  of 
the  experimental  data  was  performed  using  the  following 
spin  Hamiltonian: 


FIG.  5.  (a)  Experimental  ODMR  spectra  from 

Alo  (uGao  98^0  oi 8Po  982/GaP,  taken  at  5  K  for  both  the  X 
(9.31  GHz)  and  W  (94.8  GHz)  bands,  (b)  The  simulated  ODMR 
spectra  including  contributions  from  (c)-(e)  below,  (c)  The  simu¬ 
lated  ODMR  spectrum  for  the  “g= 2”  defect.  A  fitting  value  of  g 
=  2.017  was  used,  (d)  and  (e)  The  simulated  ODMR  spectra  for 
Ga,-A  and  Ga ,-fi,  respectively,  by  using  the  spin  Hamiltonian  pa¬ 
rameters  given  in  Table  I. 


Magnetic  Field  (mT) 

FIG.  6.  (a)  Experimental  ODMR  spectra  from 

GaN0031P0969/GaP,  taken  at  5  K  at  both  the  X  (9.31  GHz)  and  IV 
(94.8  GHz)  bands,  (b)-(e)  have  the  same  meaning  as  in  Fig.  5. 

77=  hbB  •  g  •  S  +  S  •  A  •  I.  (1) 

Here  the  first  and  second  terms  describe  electron  Zeeman 
and  central  HF  interactions,  respectively;  /jlb  is  the  Bohr 
magneton,  B  is  the  external  magnetic  field,  g  is  the  Zeeman 
splitting  tensor,  and  A  is  the  central  HF  interaction  tensor  for 
each  isotope.  The  effective  electron  spin  is  5=1/2  and  the 
nuclear  spin  is  7=3/2.  By  fitting  the  spin  Hamiltonian  to  the 
experimental  data,  the  spin  Hamiltonian  parameters  were  de¬ 
termined  for  both  Ga,-A  and  Ga,-//  in  all  samples  and  are 
given  in  Table  I.  Both  g  and  A  tensors  are  deduced  to  be 
isotropic,  within  the  experimental  accuracy,  and  are  reduced 
to  the  scalars  g  and  A.  This  finding,  together  with  the  rather 
strong  HF  interaction,  points  to  A ,  symmetry  (5  like)  for  the 
electron  wave  function  localized  at  the  defects.  For  example. 
Figs.  5-7  show  the  simulated  ODMR  curves  for 
Alo.o2Gao.98N0o12Po.988/GaP,  GaN003iPo.969/GaP,  and 
GaN0oisPo.982/Si,  assuming  a  Gaussian  line  shape  for  each 
ODMR  line.  It  is  apparent  that  the  agreement  between  the 
simulations  and  the  experimental  results  is  very  good  for 
both  X-  and  IV-band  spectra. 

The  intensity  of  the  low-field  A-band  ODMR  lines  ap¬ 
pears  to  be  weaker  in  experiment  than  in  simulation.  This 
can  be  attributed  to  the  fact  that  the  simulations  only  take 
into  account  magnetic-dipole-allowed  electron  spin  reso¬ 
nance  transitions  but  not  the  difference  in  recombination 
rates  of  the  spin  sublevels.1011  At  lower  fields  when  the  elec¬ 
tronic  Zeeman  term  becomes  comparable  or  weaker  than  the 
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FIG.  7.  (a)  Experimental  ODMR  spectra  from 

GaN0018P0982/Si,  taken  at  5  K  at  both  the  X  (9.31  GHz)  and  IV 
(94.8  GHz)  bands,  (b)  The  simulated  ODMR  spectra  for  Ga,-B. 
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TABLE  I.  Spin  Hamiltonian  parameters  determined  by  fitting  Eq.  (1)  to  experimental  data  for  the  Ga,-A  and  Ga,-B  defects  in  the  three 
types  of  samples.  Spin  density  7f  at  Ga,-  is  also  shown  and  is  deduced  from  our  calculated  4 s  electron  density,  1 04^(0) |2  =  72.7 
X  1024  cm-3.  Note  that  for  the  GaN0018P09g2/Si  and  A I  r  G  a  |  _  r  Nt  P  j  _y  /  G  a  P  (x = 0  and  y  =  0.012)  samples,  the  ODMR  signal  of  Ga,-A,  if 
existing,  is  obscured  by  the  presence  of  the  much  stronger  ODMR  signals  from  other  defects  and  thus  its  HF  interaction  parameter  cannot 
be  reliably  determined.  No  ODMR  signals  from  Ga,-A  and  Ga ,-B  were  observed  in  the  GaN00i^Po.987^GaP  sample  within  the  detection  limit 
of  our  ODMR  experiments. 


Samples 

GaNyP^ 

yl  GaP 

GaNyPj-y/Si 

AkGa, 

-,-NyPl 

_v/GaP 

X 

0 

0 

0 

0.01 

0.02 

0.3 

y 

0.031 

0.023 

0.018 

0.012 

Defect 

Ga,-A 

Ga  ,-fi 

Ga,-A 

Ga  ,-B 

Ga,-B 

Ga  ,-B 

Ga,-A 

Ga  ,-B 

Ga-A 

Ga,-B 

Ga,-A 

Ga  ,-B 

S 

1/2 

1/2 

1/2 

1/2 

1/2 

1/2 

1/2 

1/2 

1/2 

1/2 

1/2 

HI 

1 

3/2 

3/2 

3/2 

3/2 

3/2 

3/2 

3/2 

3/2 

3/2 

3/2 

3/2 

3/2 

g 

2.001 

2.003 

2.001 

2.003 

2.001 

2.005 

2.01 

2.005 

2.01 

2.005 

2.01 

2.005 

A(69  Ga)  X  10“4  cm-1 

770 

1150 

770 

1150 

1145 

1150 

490 

1030 

490 

1030 

450 

980 

A(7lGa)  X  10_4cm_1 

1000 

1480 

1000 

1480 

1477 

1480 

630 
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0.30 

0.30 

0.30 

0.13 

0.27 

0.13 

0.27 

0.12 

0.25 

HF  interaction  [see  Eq.  (1)],  the  spin  sublevels  undergo 
strong  admixing  due  to  the  HF  interaction,  which  reduces  the 
difference  in  the  recombination  rates.  As  the  ODMR  inten¬ 
sity  is  proportional1011  to  the  product  of  the  microwave- 
induced  spin  transition  probability  and  the  difference  in  the 
recombination  rates,  a  reduction  of  the  latter  should  thus  lead 
to  a  decrease  in  the  actual  ODMR  intensity. 

It  should  also  be  noted  that  the  fitting  of  the  ODMR  spec¬ 
tra  from  the  GaNP/Si  samples  in  Fig.  7  is  less  satisfactory 
by  considering  Ga rB  alone,  due  to  a  contribution  from  a 
broad  background  close  to  the  middle  of  the  spectra.  Though 
a  contribution  from  Ga,-A  to  the  background  is  possible,  it 
cannot  be  definitely  concluded  since  the  characteristic  HF 
structure  was  not  resolved  in  this  case.  Even  assuming  the 
contribution  from  Ga,-A,  additional  contributions  from  other 
unknown  defect(s)  are  still  apparent.  A  similar  situation  oc¬ 
curs  in  the  Al-free  sample  in  the  AlGaNP/GaP  series — i.e., 
the  only  other  sample  that  also  exhibits  positive  ODMR  sig¬ 
nals.  In  both  cases,  the  strong  ODMR  line  related  to  the 
“g  =  2”  defect,  which  dominates  the  ODMR  spectra  in  the 
rest  of  the  samples,  is  absent. 

B.  Defect  identification 

The  above  analysis  proves  that  both  Ga,-A  and  Ga -B  in¬ 
volve  an  atom  with  well-defined  natural  abundances  of  the 
two  isotopes  and  their  respective  nuclear  spin  numbers  and 
magnetic  moments.  Among  all  elements,  only  Ga  meets  all 
these  specifications.  Ga  has  two  naturally  abundant  isotopes: 
69Ga  (60.4%  abundant)  and  71Ga  (39.6%  abundant).  The  ra¬ 
tio  of  their  nuclear  magnetic  moments  is  /r(71  Ga)//r(69Ga) 
=  1.27.  Our  ODMR  results  thus  provide  compelling  evidence 
for  the  involvement  of  a  Ga  atom  in  the  core  of  both  Ga,-A 
and  Ga -B  defects. 

Now  let  us  turn  to  the  location  of  the  Ga  atom  in  Ga,-A 
and  Ga -B  that  gives  rise  to  the  observed  HF  structure.  A  Ga 
atom,  instead  of  residing  in  its  designated  lattice  site,  can 
reside  either  on  a  site  of  the  group-V  sublattice  forming  a 
cation  antisite  or  on  an  interstitial  position  leading  to  Ga 


self-interstitials.  A  Ga-related  cation  antisite  was  predicted 
by  earlier26  and  our  present  theoretical  calculations  to  pos¬ 
sess  an  electron  wave  function  of  7'2  symmetry,  which  will 
yield  a  very  weak  and  anisotropic  HF  interaction  at  the  de¬ 
fect  center.  This  is  clearly  inconsistent  with  our  experimental 
observation  of  the  strong  isotropic  HF  interaction  and  can 
therefore  be  ruled  out  as  a  candidate  for  the  model  of 
Ga,-A  and  Ga rB  under  study.  Our  theoretical  calculations 
also  estimate  that  the  HF  interaction  of  a  neutral  Ga,  should 
be  weaker  than  100  X  10-4  cm-1,  which  cannot  explain  the 
observed  strong  HF  interaction.  A  Ga,  in  its  Ga2+  charge 
state,  on  the  other  hand,  is  predicted  to  possess  a  strongly 
localized  electron  wave  function  of  symmetry  in  agree¬ 
ment  with  the  experimental  results.  Therefore  we  can  con¬ 
clude  that  the  Ga  atom  in  the  Ga,-A  and  Ga ,-B  defects  is  a 
Ga2+  interstitial. 

The  next  question  is  whether  both  Ga,-A  and  Ga rB  are 
just  simple,  isolated  Ga,  residing  on  different  Ga,  sites  or 
Ga, -related  complexes.  For  that  we  calculated  for  each  mi¬ 
croscopic  configuration  the  HF  interaction  A  from  the  inter¬ 
action  between  the  spin  density  of  the  electron  and  the 
nuclear  spins.  Using  state-of-the-art  first-principles  supercell 
calculations  as  implemented  in  the  WlEN2k  code,27  we  calcu¬ 
lated  the  electron  spin  distribution  topology  that  was  then 
used  to  calculate  the  HF  tensor.  The  calculations  were  based 
on  generalized  gradient  approximations  (GGA’s)  within  the 
framework  of  density  functional  theory  (DFT),  using  the  all¬ 
electron  full-potential  linearlized  augmented-plane-wave 
(FPLAPW)  method.  The  relativistic  effects  were  included 
through  a  scalar  relativistic  treatment  for  valence  electrons. 
All  atoms  in  the  supercell  are  allowed  to  relax.  The  atomic 
geometry  optimization  has  been  performed  using  the  first- 
principles  ultrasoft  pseudopotential  plane-wave  method  with 
GGA  as  implemented  in  the  VASP  code.28  Our  test  calcula¬ 
tions  of  a  well-established  signature,  Zn,  in  ZnSe,29  yielded 
the  central  HF  parameter  in  good  agreement  with  the  experi¬ 
ment  to  within  10%.  A  supercell  containing  32  atoms,  which 
was  sufficient  to  provide  converged  HF  parameters,30  was 
used.  The  electron  density  of  a  free  atom  was  calculated  by 
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placing  the  atom  in  a  large  cubic  supercell  (10X10 
X  10  A3)  to  reduce  the  intercell  interactions.  Our  theoretical 
calculations  predict  strong  HF  interactions  for  isolated  Ga,  in 
GaP:  i.e.,  A(69Ga)=  1876  X  10-4  cm-1  for  Ga,  at  the  T d  site 
surrounded  by  four  Ga  atoms  and  A(69Ga)=1599 
X  10-4  cm-1  for  Ga,  at  the  Td  site  surrounded  by  four  P  at¬ 
oms.  These  values  are  much  larger  (>30%)  than  the  experi¬ 
mental  values  determined  for  Ga,-A  and  Ga rB.  A  reduced 
HF  interaction  is  commonly  regarded  as  evidence  for  a  com¬ 
plex  due  to  charge  transfer  or  redistribution  of  the  electron 
wave  function  from  Ga,  to  the  other  partner(s)  of  the  com¬ 
plex.  This  conclusion  is  supported  by  recent  experimental 
ODMR  studies  in  wurtzite  GaN  that  have  convincingly 
shown  that  isolated  Ga  self-interstitials  are  unstable  at  tem¬ 
peratures  well  below  room  temperature  and  readily  form 
complexes.3 1-33  Considering  the  more  dilute  crystal  lattice  of 
the  GaNP  and  AlGaNP  alloys  that  makes  Ga,  easier  to  mi¬ 
grate,  the  isolated  Ga  self-interstitials  in  these  alloys  are 
more  likely  to  be  unstable  at  the  growth  temperatures,  lead¬ 
ing  to  the  formation  of  complexes.  Such  complexes  usually 
possess  larger  binding  energies,  which  are  consistent  with  the 
observed  high  thermal  stability  of  the  Ga,-A  and  Ga rB  de¬ 
fects. 

C.  Localization  of  electron  wave  functions  at  the  defects 

The  degree  of  localization  of  the  electron  wave  functions 
at  the  defect  can  be  estimated  using  a  one-electron  linear 
combination  of  atomic  orbital  (LCAO)  scheme.34  The  wave 
function  for  the  unpaired  electron  can  be  constructed  as  a 
LCAO  centered  on  the  atoms  surrounding  the  defect: 

rp  =  2  rjjtf/j.  (2) 

i 

The  wave  function  on  the  ;th  site  can  be  approximated  as  a 
hybrid  of  the  ns,  np  orbitals: 

4'i  ~  ,) !  "1  (3) 

where  n  =  4  for  Ga.  The  normalization,  ignoring  the  overlap, 
requires 

ar  +  /3,2=l,  2  Vi  =  1  •  (4) 

i 

Generally,  with  an  axially  symmetric  A  tensor,  the  principal 
values  can  be  described  as 

Ay  =  a i  +  2  bh  (5) 

A'±  =  a j  -  bh  (6) 

where 

a,  =  TTgogN^BM-Na,2  V2\4aM\2i  0) 

2 

b,  =  -gQgNV-BV-NP,2  y,2(r3p-3)r  (8) 

Here  g0  and  gN  denote  electron  and  nuclear  g  values.  In  our 
case,  the  HF  tensor  of  the  Ga,  site  is  shown  to  be  isotropic — 


i.e.,  A°a=A(jia  and  thus  a^a=  1  and  /3(2ia=0.  The  wave  func¬ 
tion  of  the  unpaired  electron  at  the  Ga,-A  and  Ga rB  defects 
can  then  be  approximated  as  being  constructed  from  the  s 
orbital  only. 

For  a  free  neutral  Ga  atom,  we  calculated  the  charge  den¬ 
sity  of  the  4 s  electron  using  the  relativistic  GGA-FP-LAPW 
method  described  in  Sec.  IV  B  and  found  |f/qs(0)|2=72.7 
X  1024  cm-3.  Based  on  this,  we  estimate  the  localization  of 
the  unpaired  electron  at  the  Ga,  to  be  in  the  range  of  12%- 
30%  (Table  I).  The  rather  strong  localization  confirms  that 
these  are  deep-level  defects. 

D.  Local  surrounding  of  the  defects 

In  a  tetrahedral  semiconductor  like  GaNP  or  AlGaNP, 
self-  interstitials  can  occupy  three  high-symmetry  positions 
in  which  two  of  them  have  Td  symmetry  surrounded  by 
group-III  or  group-V  atoms,  respectively,  and  the  third  has 
hexagonal  symmetry  (D3d)  surrounded  by  both  group-III  and 
group-V  sublattices.  In  order  to  provide  information  on  the 
specific  Ga,-  site  and  the  sublattice  immediately  surrounding 
the  Ga,  in  Ga,-A  and  Ga -B,  the  chemical  compositions  of 
each  sublattice  were  varied  experimentally  by  varying  Al  and 
N  compositions.  The  results  are  summarized  in  Table  I.  It 
can  be  seen  that  a  rather  strong  reduction  (about  10%  for 
Ga rB)  of  the  HF  interaction  was  observed  when  the  Al  com¬ 
position  is  raised  from  v=( ).()()  to  as  low  as  .v = ().()  I .  With  a 
further  increase  in  Al  composition,  however,  the  rate  of  de¬ 
crease  in  the  HF  interaction  is  drastically  reduced  to  about 
5%  from  x=0.01  to  x=0.30.  This  observation  indicates  that 
Al  directly  affects  the  Ga,-A  and  Ga rB  defects.  The  nonlin¬ 
ear  dependence  also  indicates  that  the  presence  of  Al  near  the 
defects  may  not  necessarily  follow  the  statistic  distribution 
expected  for  random  alloy  formation.24  This  suggestion  is 
further  supported  by  the  fact  that  the  ODMR  intensity  in¬ 
creases  significantly  upon  incorporation  of  Al  in  the  alloy. 
We  should  point  out  that  the  incorporation  of  Al  in  GaAs  was 
also  required  to  introduce  the  Ga,  defect  in  AlGaAs.35-36 
Moreover,  the  strength  of  the  HF  interaction  in  the 
(Al,Ga)As-based  structures35-38  is  very  close  to  that  of 
Ga,-A  in  the  Al-containing  AlGaNP  studied  in  this  work. 
Assuming  that  a  similar  Ga,  defect  complex  is  involved,  the 
correlation  between  these  two  alloy  systems  seems  to 
suggest  that  the  Ga,  in  Ga,-A  is  surrounded  by  the  group-III 
sublattice  as  the  group-V  sublattice  is  made  of  completely 
different  atoms  between  AlGaAs  and  AlGaNP.  In  the 
previous  studies  of  Ga,  in  AlGaAs,35  it  was  suggested  that 
Ga,  is  not  very  near  the  group-III  sublattice  based  on  the 
experimental  observation  that  the  ODMR  spectrum  is  insen¬ 
sitive  to  the  change  of  Al  composition  from  20%  to  50%. 
Our  finding  of  the  nonlinear  dependence  of  the  HF  interac¬ 
tion  on  Al  compositions,  which  exhibits  a  much-reduced 
variation  for  high  Al  compositions,  questions  this  earlier 
suggestion. 

In  contrast  to  the  dependence  on  Al  composition,  the  HF 
interaction  of  both  Ga,-A  and  Ga -B  defects  in  GaNP 
was  found  in  this  work  to  be  insensitive  to  N  compositions 
up  to  3%.  On  the  other  hand,  the  incorporation  of  N  is  es¬ 
sential  for  the  observation  of  the  Ga,-A  and  Ga ,-B  defects  in 
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GaNP.  This  indicates  that  either  N  is  directly  involved  as  part 
of  the  defect  or  the  required  nonequilibrium  growth  condi¬ 
tions  or  possible  ion  bombardment  during  the  incorporation 
of  N  in  GaNP  facilitates  the  formation  of  the  Ga,-A  and 
Ga rB  defects.  It  is  interesting  to  note  that  the  strength  of  the 
HF  interaction  of  Gar,4  in  GaNP,  A(69  Ga)  =  770 
X  10-4  cm-1,  is  very  close  to  741  X  10-4  cm-1  observed  for 
the  Ga,-  in  its  parent  GaP.39  The  same  applies  to  Ga,-A  in 
AlGaNP;  i.e.,  A(69Ga)  =  (450-490)  X  1CT4  cm-1  is  very  close 
to  500  X  10-4  cm-1  observed  in  the  AlGaAs-based 
alloys.35"38  In  contrast,  Ga,-/i  was  never  observed  in  N-free 
GaP  and  AlGaAs-based  alloys.  This  may  indicate  that 
N  could  be  directly  involved  in  the  structure  of  the  Ga -B 
defect. 


E.  Recombination  processes 

The  ODMR  signals  were  detected  as  changes  of  the  PL 
intensity  induced  by  the  spin  resonance.  The  sign  of  the 
ODMR  signals  can  thus  provide  us  with  useful  information 
about  the  role  of  the  involved  defects  in  carrier  recombina¬ 
tion  processes. 

In  the  GaNP/GaP  series  of  samples,  in  Figs.  1  and  6, 
the  ODMR  signals  correspond  to  a  decrease  upon  spin 
resonance  in  the  intensity  of  all  monitored  PL  emissions 
over  the  entire  spectral  range.  Such  a  negative  and  emission- 
insensitive  ODMR  signal  has  often  been  taken  as  a  signature 
of  a  nonradiative  defect,10,39  which  is  not  directly  related 
to  the  monitored  PL  emissions.  The  reason  why  a  nonradia¬ 
tive  defect  can  be  detected  via  a  radiative  recombination 
process  lies  on  their  competition  in  recombination  of 
photogenerated  free  carriers.  Figure  8(a)  shows  schemati¬ 
cally  the  recombination  processes  leading  to  the  observation 
of  the  negative  ODMR  signals.  Spin-resonance-enhanced 
recombination  via  a  nonradiative  defect  leads  to  a  reduction 
of  free  carriers  available  for  radiative  recombination  and 
thus  to  a  decrease  in  the  intensity  of  PL  emissions.10,11,25,40 
Such  nonradiative  recombination  can  occur  via  a  single 
defect  or  via  intercenter  charge  transfer41  between  two 
defects — e.g.,  between  the  “g  =  2”  defect  and  Ga,-A 
(or  Ga -B).  The  latter  scenario  is  quite  likely  in  the  present 
case,  judging  from  the  simultaneous  appearance  of  these 
two  defects  for  each  sample  under  various  measurement 
conditions.  Therefore  it  can  be  concluded  that  Ga,-A,  Ga ,-B, 
and  the  “g  =  2”  defect  act  as  nonradiative  recombination 
channels,  reducing  the  radiative  efficiency  of  the  alloys. 
This  is  confirmed  by  an  anticorrelation  between  the 
intensities  of  the  visible  PL  emissions  and  the  ODMR 
signals,  in  Fig.  9.  For  example,  a  significant  increase  in  the 
ODMR  intensity  of  the  Ga,-A  and  Ga -B  defects  (and  also 
the  “g  =  2”  defect)  can  be  observed  with  increasing  N 
composition,  which  is  accompanied  by  a  substantial  decrease 
of  the  visible  PL  emissions.  These  defects  can  undermine 
performance  of  potential  photonic  and  optoelectronic 
devices  and  therefore  should  be  addressed  adequately. 
The  observed  increase  of  the  NIR  PL  in  Fig.  9,on  the 
other  hand,  can  be  explained  by  the  fact  that  the  NIR  PL 
originates  from  defects  that  were  introduced  by  the  incorpo¬ 
ration  of  N. 
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defect  (Ga.-B)  defect 


FIG.  8.  (a)  A  schematic  diagram  of  the  carrier  recombination 
processes  leading  to  the  observation  of  negative  ODMR  signals  in 
GaNP/GaP.  The  spin  resonance  induces  an  increase  in  the  nonra¬ 
diative  recombination  (indicated  by  the  dashed  arrows)  via  the  “g 
=  2”  and  Ga,-A  (or  Ga  ,-B)  defects  and  a  decrease  in  the  radiative 
recombination  process  monitored  by  the  PL  in  both  the  visible  and 
NIR  spectral  ranges,  (b)  A  schematic  diagram  of  the  carrier  recom¬ 
bination  process  leading  to  the  observation  of  positive  ODMR  sig¬ 
nals  in  GaNP/Si;  i.e.,  the  spin-resonance-induced  increase  in  the 
radiative  recombination  via  the  unknown  defect  and  Ga,-A  (or  Ga, 
-B)  is  directly  monitored  by  the  NIR  PL.  (c)  A  schematic  diagram  of 
the  competitive  carrier  recombination  processes  leading  to  the  ob¬ 
servation  of  positive  and  negative  ODMR  signals  in  the 
AlGaNP/GaP  series  of  samples. 

In  the  GaN0018P09g2/Si  samples,  the  ODMR  signals 
were  found  to  be  positive  by  monitoring  the  NIR  PL 
emissions  (Figs.  1  and  7).  This  can  be  explained  as  follows: 
The  monitored  PL  emission  is  directly  related  to  the 
defects — e.g.,  intercenter  charge  transfer  between  an  un¬ 
known  defect  and  the  Ga,-A  (or  Ga,--B)  defect  [see  Fig.  8(b)]. 
The  spin-resonance-enhanced  radiative  recombination 
between  them  leads  to  an  intensity  increase  of  the 
corresponding  PL.  It  should  be  pointed  out  that  the  “g= 2” 
defect  is  absent  in  this  case.  Instead,  a  broad  background 
ODMR  was  detected,  as  shown  in  Fig.  7,  of  which  the  origin 
of  the  corresponding  defect  is  still  unknown.  It  is  denoted 
as  “unknown  defect”  in  Fig.  8(b).  The  energy  of  the  related 
PL  emission  can  usually  be  employed  to  estimate  the  energy 
position  of  the  defect  level  if  the  recombination  process 
is  fully  understood.  Since  the  spin-dependent  recombination 
processes  involving  the  Ga,  defects  are  likely  to  be  of 
an  intercenter  charge  transfer  origin  and  the  energy  position 
of  the  other  partner  is  unknown,  at  present  we  are  not  able 
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FIG.  9.  Relative  intensities  of  the  PL  emissions  and  the  ODMR 
signals  as  a  function  of  the  N  composition,  obtained  from  the 
GaNP/GaP  series  of  samples.  The  open  and  solid  circles  denote  the 
intensities  of  the  visible  and  NIR  PL  emissions,  respectively.  The 
solid  and  open  triangles  represent  the  intensities  of  the  Ga,-A  and 
Ga-B  ODMR  signals,  respectively.  The  solid  squares  denote  the 
ODMR  intensity  of  the  “g  =  2”  defect. 

to  precisely  deduce  the  energy  level  positions  of  the 
Ga,  defects.  However,  a  rough  estimate  of  the  energy 
position  for  the  Ga -B  defect  can  be  made  to  be  at  least 
>1.2  eV  away  from  either  the  conduction  or  valence  band 
edge. 

In  the  AlGaNP/GaP  series  of  samples,  on  the  other  hand, 
both  positive  and  negative  ODMR  signals  have  been  ob¬ 
served,  depending  on  the  Al  compositions  and  the  monitored 
PL.  In  the  Al-free  sample,  positive  ODMR  signals  similar  to 
that  from  the  GaN0.oi8Po.982/Si  sample  were  detected  by 
monitoring  the  NIR  PL.  The  recombination  process  can  thus 
be  described  in  a  similar  fashion  as  discussed  above  and 
shown  in  Fig.  8(b).  In  the  Al-containing  AlGaNP,  however, 
the  ODMR  signals  switch  sign  to  become  negative.  This  is 
accompanied  by  the  appearance  of  the  “g= 2”  defect  and  a 
simultaneous  disappearance  of  the  “unknown  defect”  in  the 
ODMR  spectra,  as  shown  in  Fig.  5.  This  sign  change  can  be 
understood,  as  schematically  shown  in  Fig.  8(c),  by  the  fol¬ 
lowing  competing  intercenter  charge  transfer  and  recombina¬ 
tion  processes  between  these  defects:  (1)  the  nonradiative 
recombination  between  the  “g  =  2”  and  Ga-A  (or  Ga -B)  de¬ 
fects,  which  gives  rise  to  negative  ODMR  signals,  and  (2) 
the  radiative  recombination  between  the  “unknown 
defect”  and  the  Ga-A  (or  Ga rB)  defect,  which  gives  rise  to 
positive  ODMR  signals.  In  principle,  both  processes  coexist 
and  the  sign  of  the  combined  ODMR  signals  depends  on 
their  relative  importance.  The  change  of  the  signal  sign  from 
positive  in  the  Al-free  sample  to  negative  in  the  Al- 
containing  samples  as  soon  as  Al  is  incorporated  suggests 
that  process  (1)  becomes  dominant  over  process  (2). 
With  increasing  Al  composition,  the  intensity  of  the  negative 
ODMR  signals  in  Fig.  10  increases,  signifying  the  impor¬ 
tance  of  the  nonradiative  recombination  between  the  “g= 2” 
and  Ga,-A  (or  Ga-/i)  defects.  With  a  further  increase  of 
the  Al  composition  to  30%,  however,  the  ODMR  intensity 
decreases  for  all  defects.  This  is  a  sign  that  other  nonradia¬ 
tive  defects  become  important  and  overshadow  the  above 
defects. 

It  should  be  pointed  out  that,  although  the  intercenter 
charge  transfer  is  shown  in  Fig.  8  only  from  the  “g= 2”  de¬ 


FIG.  10.  Relative  intensities  of  the  PL  emissions  and  the  ODMR 
signals  as  a  function  of  the  Al  composition,  obtained  from 
the  AlGaNP/GaP  series  of  samples.  The  legends  are  the  same  as  in 
Fig.  9. 

feet  (or  the  “unknown  defect”)  to  the  Ga-A  (or  Ga -B)  de¬ 
fect,  the  reverse  process  is  also  possible  if  their  energy  order 
is  reversed.  This  possibility  holds  as  long  as  the  paramag¬ 
netic  charge  state  of  Ga,  in  the  Ga-A  (or  Ga-/i)  defect,  be¬ 
fore  the  charge  transfer,  is  2+. 

V.  CONCLUSIONS 

In  summary,  we  have  studied  and  identified  two  different 
Ga,  defects  in  dilute  nitride  AlGaNP  by  ODMR.  The  assign¬ 
ment  of  Ga,  to  both  defects  is  concluded  from  its  unique  HF 
interaction.  Both  defects  are  complexes  involving  a  Ga,2+. 
The  redistribution  of  the  electron  wave  function  leads  to  a 
reduction  of  the  HF  interaction.  The  HF  interaction  also 
strongly  depends  on  the  presence  of  Al,  indicating  that  the 
defects  may  be  surrounded  by  group-HI  atoms.  The  incorpo¬ 
ration  of  both  Al  and  N  is  shown  to  play  an  important  role  in 
the  formation  of  the  defects,  indicative  of  the  involvement  of 
N  in  Ga -B.  The  defects  are  found  to  introduce  energy  levels 
deep  inside  the  band  gap — e.g.,  —1.2  eV  from  the  conduc¬ 
tion  (or  valence)  band  edge  for  the  Ga ~B  defect.  In  most 
cases,  these  defects  act  as  nonradiative  recombination  centers 
reducing  the  efficiency  of  light  emissions  from  the  alloys. 
They  can  thus  undermine  the  performance  of  photonic  de¬ 
vices  based  on  these  materials.  In  view  of  the  high  thermal 
stability  (at  up  to  900  °C)  of  the  defects,  post-growth  RTA 
was  found  to  be  inadequate  to  remove  the  defects.  Future 
studies  are  needed  to  control  their  formation  during  the 
growth,  in  order  to  improve  radiative  efficiency.  The  present 
studies  also  shed  new  light  on  Ga  self-interstitials  in  III-V’s 
in  general. 
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Recent  theory  has  found  that  native  defects  such  as  the  O  vacancy  V0  and  Zn  interstitial  Zn;  have  high 
formation  energies  in  n-type  ZnO  and,  thus,  are  not  important  donors,  especially  in  comparison  to 
impurities  such  as  H.  In  contrast,  we  use  both  theory  and  experiment  to  show  that,  under  N  ambient,  the 
complex  Zn/-N0  is  a  stronger  candidate  than  H  or  any  other  known  impurity  for  a  30  meV  donor 
commonly  found  in  bulk  ZnO  grown  from  the  vapor  phase.  Since  the  Zn  vacancy  is  also  the  dominant 
acceptor  in  such  material,  we  must  conclude  that  native  defects  are  important  donors  and  acceptors  in 
ZnO. 


DOI:  10. 1103/PhysRevLett.95. 225502 

Semiconducting  ZnO  has  generated  great  interest  in  the 
past  decade  because  of  advances  in  bulk  and  epitaxial 
growth  that  have  opened  the  door  for  new  photonic  and 
electronic  applications,  such  as  UV  light  emitting  diodes 
and  transparent  transistors  [1,2].  Previous  research  had 
established  that  ZnO  was  always  n-type  and  that  the  domi¬ 
nant  donors  were  usually  shallow  with  activation  energies 
of  between  30  and  60  meV  [3,4].  Because  it  was  also 
known  that  the  crystal  growth  was  typically  Zn-rich,  the 
dominant  donor  was  almost  always  identified  as  either  the 
O  vacancy  V()  or  the  Zn  interstitial  Zn7  [5,6].  This  model 
was  strongly  challenged  in  the  year  2000  when  Kohan  et  al. 
showed  theoretically  that  both  VQ  and  Zn /  have  high  for¬ 
mation  energies  in  n-type  ZnO  and  that,  furthermore,  both 
are  deep,  not  shallow,  donors  [7].  More  recent  theory  has 
concluded  that  Zn7  is  actually  a  shallow  donor,  rather  than 
deep  [8,9],  as  has  also  been  suggested  by  electron- 
irradiation  experiments  [6];  however,  its  high  formation 
energy  would  still  limit  its  participation  in  the  conductivity 
of  re-type  material.  Also  in  the  year  2000,  the  defect-donor 
model  was  further  challenged  by  Van  de  Walle’s  theoreti¬ 
cal  result  that  H  is  always  a  donor  in  ZnO,  that  it  is  easily 
ionized,  and  that  it  has  a  low  enough  formation  energy  to 
be  abundant;  thus.  Van  de  Walle  suggested  that  it  was 
likely  to  be  a  dominant  background  donor  in  ZnO  materials 
that  were  exposed  to  H  during  growth  [10].  This  proposal 
has  been  amenable  to  testing,  because  H-containing,  high- 
quality,  bulk  ZnO,  grown  by  a  seeded  chemical  vapor 
transport  (SCVT)  technology,  has  been  commercially 
available  for  the  past  few  years.  For  the  most  part,  these 
tests  have  confirmed  that  a  shallow  donor  due  to  H  exists  in 
SCVT  ZnO  and  can  contribute  significantly  to  the  con¬ 
ductivity  [11-16].  This  fact,  coupled  with  the  theoretical 
evidence  of  high  formation  energies  for  the  native  donors 
[7],  has  led  to  a  prevailing  opinion  that  native  donors  do  not 
play  a  significant  role  in  the  conductivity  of  as-grown  ZnO. 
In  contrast,  we  will  offer  evidence  here  that  native  donors 


PACS  numbers:  61.72.Ji,  71.55.Gs,  72.20.Fr,  78.55.Et 

can  contribute  significantly  to  conduction  in  ZnO  but  as 
complexes,  rather  than  isolated  elements. 

The  main  sample  used  in  this  study  was  a  5  mm  X 
5  mm  X  0.42  mm  piece  cut  from  a  wafer  grown  by  the 
SCVT  technique  at  ZN  Technology,  Inc.  [17].  Material  of 
this  type  is  of  very  high  quality,  with  peak  electron  mobil¬ 
ity  >2000  cm2/V  s,  300  K  carrier  concentration  in  the 
1016  cm-3  range,  and  photoluminescence  (PL)  linewidths 
under  0.5  meV,  as  shown  in  Figs.  1  and  2  [4].  The  sample 
was  first  annealed  at  7 15  °  C,  in  order  to  release  most  of  the 
H  [14,15],  then  was  irradiated  with  high-energy  electrons 
to  create  point  defects,  and  finally  was  annealed  again  at 
temperatures  from  200-500  °C,  in  order  to  investigate  the 
annihilation  of  the  point  defects.  As  seen  in  Fig.  1,  several 
sharp  PL  lines  appear  in  the  region  3.357-3.365  eV,  and 
these  spectral  features  are  usually  assigned  to  transitions  of 
excitons  bound  to  neutral  donors  (D{)X  transitions),  in 
which  the  donor  remains  in  its  ground  state  (re  =  1)  during 
the  transition.  Some  of  these  lines  have  been  tentatively 
identified;  for  example,  the  lines  at  3.35964  eV  (/8  or  ISa  in 


3.34  3.35  3.36 

E(eV) 


FIG.  1  (color  online).  4-K  photoluminescence  spectra  for  ZnO 
sample.  The  inset  shows  the  D°X  lines  in  greater  detail. 
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the  literature)  and  3.36042  eV  (/6  or  /6a)  have  been  as¬ 
signed  to  GaZn  and  AlGn,  respectively.  (See  Ref.  [18]  for  an 
excellent  review  of  PL  in  ZnO.)  All  of  the  spectra  pre¬ 
sented  in  this  work  are  normalized  to  the  /8  line,  because  it 
is  relatively  isolated  and  also  not  expected  to  change  sig¬ 
nificantly  as  a  result  of  annealing  or  irradiation  treatments. 
However,  the  most  dominant  PL  line  in  this  particular 
sample  before  annealing  is  the  D°X  line  at  3.36270  eV 
(74),  now  almost  universally  assigned  to  interstitial  H 
[16,18].  This  identification  results  at  least  partially  from 
annealing  experiments,  because  /4  disappears  for  anneals 
above  600  °C,  in  good  correlation  with  the  effusion  of  H 
from  the  sample  [14-16,19].  Our  sample  behaves  in  the 
same  way,  as  evidenced  by  the  strong  reduction  of  /4  after 
an  anneal  of  715  °C  (curve  2  in  Fig.  1).  Not  shown  in  Fig.  1 
is  a  two-electron  satellite  (TES)  replica  of  /4,  appearing  at 
3.32961  eV.  A  TES  transition  is  one  in  which  the  donor  is 
left  in  an  excited  n  =  2  state  after  the  collapse  of  the 
exciton.  If  the  donor  associated  with  I4  follows  a  hydro- 
genic  model,  then  the  ground-state  ( n  =  1)  energy  of  this 
donor  should  be  given  by  4/3(3.36270  —  3.32961)  = 
44.1  meV.  Other  TES  transitions,  also  not  shown,  are 
seen  near  3.32  eV  and  may  be  associated  with  I6  and  /8 
[18].  If  so,  then  the  hydrogenic  model  would  predict  - 
ground-state  energies  of  53-55  meV  for  these  donors, 
presumably  associated  with  GaZn  and  AlZn.  Thus,  the 
donors  identified  from  PL  are  H,  at  44  meV,  and  GaZn 
and  AlZn,  at  about  55  meV. 

To  create  point  defects,  we  have  used  the  Van  de  Graaff 
electron  accelerator  at  Wright  State  University.  The  sample 
was  irradiated  on  the  Zn  face  (0001)  with  1  MeV  electrons 
three  separate  times,  making  a  total  fluence  of  3  X 
1017  cm~2.  The  PL  spectrum  after  the  third  irradiation  is 
presented  as  curve  3  in  Fig.  1 .  A  new  sharp  defect  line  at 
3.36070  eV,  which  we  will  designate  as  ID,  is  generated  by 
the  irradiation.  Concomitantly,  a  triplet  feature,  comprised 
of  energies  3.33711,  3.33793,  and  3.33840  eV,  is  also 
generated.  These  we  will  designate  ID iTES1,  Id.tesi,  and 
ID tes3-  This  triplet  has  an  intensity  about  100  times  less 
than  that  of  Id,  and  this  reduction  is  about  the  same  as  that 
observed  for  the  TES  line  of  H  compared  with  its  parent 
line,  /4.  Thus,  the  defect  triplet  at  3.338  eV  clearly  is 
related  to  ID,  and,  if  the  hydrogenic  model  holds,  the 
associated  energy  is  4/3(3.3607  —  3.3379)  ~  30.4  meV. 
This  value  is  less  than  that  predicted  from  an  empirical 
version  of  Haynes’  rule  presented  in  Ref.  [18].  However, 
that  version  was  developed  for  simple,  isolated  impurities, 
Al,  Ga,  In,  and  H,  and  there  is  no  reason  to  expect  that  it 
should  hold  for  a  defect-related  complex.  Both  Id  and  its 
TES  triplet  are  greatly  reduced  for  anneals  greater  than 
500  °C  (not  shown),  evidently  due  to  defect  recombination. 
It  is  important  to  note  at  this  point  that  the  as-grown 
sample,  represented  by  curve  1  in  Fig.  1,  has  two  small 
lines  at  the  energies  of  /onS|  and  ID XES2,  respectively,  and 
possibly  also  a  line  in  the  ID  region,  although  obscured  by 
I6.  This  is  our  first  indication  that  the  as- grown  sample 
contains  defect-related  donors. 


We  next  discuss  the  temperature-dependent  mobility  j± 
and  carrier  concentration  n,  shown  in  Fig.  2.  To  avoid 
clutter,  we  show  curves  representing  only  four  stages  in 
the  evolution  of  this  sample:  (1)  as-grown;  (2)  annealed  at 
715  °C  for  1/2  hour  in  flowing  Ni  gas;  (3)  irradiated  with 
1  MeV  electrons,  in  three  equal  stages  up  to  a  total  fluence 
of  3  X  1017  cm-2;  and  (4)  annealed  at  400  °C,  following 
previous  anneals  beginning  at  200  °C.  It  should  be  noted 
that  a  final  anneal  at  500  °C  (not  shown)  eliminated  almost 
all  of  the  irradiation  damage  and  basically  returned  the  /ul 
and  n  curves  to  those  of  stage  2.  [Note  that  a  return  to  the  /jl 
and  n  curves  of  stage  1  (as-grown)  is  of  course  impossible, 
because  the  H  content  was  lost  in  stage  2,  the  first  715  °C 
anneal.] 

The  mobility  data  (inset  in  Fig.  2)  were  fitted  to  an 
accurate  charge-carrier  scattering  theory,  described  else¬ 
where  [20],  and  the  only  fitting  parameter  was  the  acceptor 
concentration  NA.  The  carrier  concentration  data  of  Fig.  2 
were  then  fitted  to  the  charge -balance  equation: 


n  +  N 4  — 


Ft 


Nd, 


+  «/ <t>Di  ’ 


(1) 


where  the  subscript  i  denotes  a  particular  donor  and  where 
4>  D,  is  a  function  of  T  and  EDi,  the  donor  energy  [cf.  Eq.  (8) 
of  Ref.  [20]].  The  fitting  parameters  in  Eq.  (1)  are  the  donor 
concentrations  (NDi’s),  donor  energies  (EDj’ s),  and  accep¬ 
tor  concentration  NA.  Excellent  fits  to  the  n  vs  l/T  data  in 
Fig.  2  were  obtained  by  also  including  a  degenerate  surface 
layer  in  the  analysis  [20,21]. 

The  PL  results  discussed  above  determined  the  energies 
of  three  different  donors,  calculated  from  their  respective 
TES  lines:  55  meV,  possibly  associated  with  GaZn  and/or 
AlZn;  44  meV,  associated  with  H;  and  30  meV,  produced  by 
irradiation  and  thus  associated  with  a  defect.  Indeed,  the  30 
and  44  meV  energies  turn  out  to  be  good  fitting  parameters 
for  two  of  the  three  donors  required  to  fit  our  Hall-effect 
data,  but  the  third  donor  is  best  fitted  with  about  75  meV, 
rather  than  55  meV.  The  carrier  concentration  fits  are 


103/T  (K'1) 

FIG.  2  (color  online).  Temperature-dependent  carrier  concen¬ 
tration  for  ZnO  sample.  The  solid  lines  are  theoretical  fits.  The 
inset  shows  the  experimental  mobility  curves  after  the  same  four 
treatments. 
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shown  as  solid  lines  in  Fig.  2,  and  the  fitting  parameters  are 
given  in  Table  I.  Two  major  conclusions  are  evident  from 
the  results:  (1)  the  44  meV  H  level  essentially  disappears 
during  the  715  °C  anneal,  in  agreement  with  the  PL  data; 
and  (2)  a  30  meV  level  is  produced  by  irradiation  but  also 
exists  in  the  as-grown  sample  and  the  sample  annealed  at 
715  °C.  It  should  be  emphasized  that  the  choice  of  30  meV 
as  a  Hall  fitting  parameter  is  not  dependent  on  the  fact  that 
the  PL  analysis  also  found  a  donor  at  30  meV.  Indeed,  a 
donor  of  this  energy  has  been  found  previously  by  a 
number  of  groups  to  give  good  fits  to  ZnO  Hall-effect 
data  [3,4,6].  To  help  identify  this  30  meV  donor,  we  appeal 
to  theory. 

First  of  all,  we  calculate  the  expected  1  MeV-electron- 
bombardment  production  rates  of  Zn  and  O  Frenkel 
pairs  from  molecular  dynamics  simulations  [22].  The 
threshold  energies  were  calculated  by  giving  a  randomly 
chosen  O  or  Zn  atom  a  recoil  energy  in  a  random  direc¬ 
tion  in  an  experimentally  controlled  angular  window  of 
15  degrees  around  the  desired  (0001)  direction.  The  inter¬ 
atomic  interaction  model  for  the  ZnO  system  will  be 
published  elsewhere  [23],  but  the  potential  develop¬ 
ment  principles  are  described  in  Ref.  [24].  To  account  for 
the  experimental  situation  with  a  beam  acceptance  angle, 
the  threshold  was  determined  as  the  average  over  the 
direction-specific  thresholds  obtained  within  the  15°  an¬ 
gular  window. 

For  Zn-face  (0001)  irradiation  at  300  K,  the  threshold  for 
O  displacement  is  44  eV,  and  that  for  Zn  displacement, 
34  eV.  We  then  apply  the  McKinley-Feshbach  relativistic 
cross-section  formula  [25]  to  give  effective  production 
rates  of  0.18  cm-1  for  O  displacement  and  0.30  cm-1  for 
Zn  displacement.  Thus,  we  would  expect  that  our  fluence 
of  3  X  1017  cm-2  would  produce  0/  and  VQ  concentra¬ 
tions  of  about  5  X  1016  cm-3  and  Zn /  and  VZn  concentra¬ 
tions  of  about  9  X  1016  cm-3.  Either  of  these  numbers  is 
consistent  with  the  observed  30  meV  donor  concentration 
of  about  2  X  1017  cm-3  in  the  irradiated  sample  (see 
Table  I),  because  the  value  2  X  1017  cm“3  is  actually  an 
upper  limit.  That  is,  the  mobility  after  such  a  heavy  irra¬ 
diation  is  almost  certainly  reduced  by  electrical  inhomo¬ 
geneity,  and,  thus,  it  is  artificially  low,  leading  to  artificially 
high  donor  and  acceptor  concentrations.  However,  the 
concentration  of  Zn;  should  still  be  about  twice  that  of 
V0,  and,  moreover,  there  is  abundant  evidence  that  VQ  is  a 
deep,  not  shallow,  donor  [7-9,26].  Thus,  Zn 7  is  a  much 
better  candidate  for  the  irradiation  donor  than  VQ. 
However,  there  also  is  evidence  that  isolated  Zn7  is  mobile 


at  room  temperature,  so  that  it  likely  has  to  form  a  complex 
to  be  stable  [27,28]. 

A  search  for  potential  complexing  partners  for  Zn7 
immediately  suggests  N,  which  easily  substitutes  for  O  in 
the  ZnO  lattice  and  which  indeed  has  a  concentration  of 
about  1017  cm-3  in  ZnO  material  of  the  type  we  are  using 
[29].  Thus,  we  have  examined  the  formation  energy,  bind¬ 
ing  energy,  and  (0/  +  )  transition  energy  of  the  complex 
Zn/-N0  using  first  principles  calculations.  We  applied 
density  functional  theory  (DFT)  within  the  local  density 
approximation  (LDA)  and  used  Vanderbilt-type  ultrasoft 
pseudopotentials,  as  implemented  in  the  VASP  code  [30]. 
To  obtain  defect  formation  energies,  defined  elsewhere 
[31,32],  a  supercell  approach  was  used,  with  a  ZnO  super- 
cell  size  of  96  atoms.  The  main  result  of  the  calculation  is 
that  the  binding  energy  of  the  Zn/-N0  complex  is  about 
0.9  eV;  thus,  the  complex  should  be  stable  at  room  tem¬ 
perature.  The  formation  energies  of  the  N0,  Zip,  and 
Zrp-N0  species  in  ZnO  depend  on  the  Fermi  energy  and 
the  partial  pressures  of  the  elements  or,  in  other  words,  the 
chemical  potentials  of  the  elements,  during  growth.  The 
secondary-ion  mass-spectroscopy  measurements  [29] 
found  a  N  concentration  [N]  of  about  1  X  1017  cm-3. 
Based  on  our  calculated  formation  energies  of  N0,  Zip, 
and  Zrp-N0,  together  with  the  measured  N  concentration 
([N]  =  [N0]  +  [Zn/-N0]),  we  estimate  a  N  chemical  po¬ 
tential  of  0.92  eV  below  the  N2  precipitation  limit  (assum¬ 
ing  Zn-rich  and  charge-neutral  growth  conditions  and  a 
growth  temperature  of  950  °C).  The  calculated  formation 
energy  [33]  shows  that,  during  growth,  N0“  acts  as  the 
dominant  acceptor  and  Zrp2+  as  the  dominant  donor,  with 
the  Fermi  energy  pinned  at  about  1 .0  eV  above  the  valence 
band  where  the  two  defects  have  approximately  the  same 
energy,  about  1.4  eV.  The  corresponding  defect  concen¬ 
trations  are  on  the  order  of  1017  cm-3.  The  concentration 
of  the  Zip-N0  complex  is  about  2  orders  of  magnitude 
lower,  as  determined  from  the  reaction  N0  +  Zn7  — » 
Zrp-N0  +  0.9  eV  and  the  associated  detailed-balance 
relationship  [NoJZn/J/NsitgfZnz-No]  =  exp(—Eb/kT). 
Here  we  use  Nsite(ZnO)  =  4.28  X  1022  cm-3  and  Eh  = 
0.9  eV.  At  T  =  950  °C,  we  get  [ZnrN0]/[Zn/]  =  0.012, 
which  means  that  only  about  1/100  of  the  Zn /  ions  take 
part  in  formation  of  the  complex.  However,  during  cool¬ 
down,  it  is  expected  that  the  formation  of  the  complex 
will  accelerate,  assuming  a  sufficiently  low  Zn7  diffusion 
barrier.  As  the  temperature  cools  to  about  500  °  C,  the 
quantity  [Ztp-N0]/[Zip]  approaches  unity,  which  means 
that  about  half  of  the  Zn7  are  already  bound  in  the  com- 


TABLE  I.  Fitting  parameters  for  mobility  and  carrier  concentration  data. 


Em  (meV) 

NDi  ( 1 0 1 6  cm  3) 

Em  (meV) 

Nm  (1016cm“3) 

Em  (meV) 

VD3  (1016cm~3) 

Na  (1016cm-3) 

As-grown 

30 

0.45 

44 

3.0 

75 

2.0 

0.13 

715  °C 

30 

0.74 

75 

3.2 

0.7 

Irradiated 

30 

—20 

75 

0.13 

-20 

400  °C 

30 

1.35 

75 

0.4 

1.2 
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plexes.  At  room  temperature,  [Zn7-N0]/[Zn7]  ~  109,  so 
that  nearly  100%  of  the  Zn7  are  in  complexes.  With  other 
impurities  (such  as  H7)  in  the  sample,  the  above  balance 
conditions  would  vary  somewhat,  but  the  main  conclusions 
should  remain  the  same. 

The  electronic  transition  energy  of  ZnrN0  is  somewhat 
uncertain,  due  to  the  well-known  LDA  gap  error.  However, 
our  calculations  show  that  the  donor  level  of  Zn7-NQ  is 
shallower  than  that  of  isolated  Zn7,  which  itself  is  argued  to 
be  a  shallow  donor  [8,9].  In  addition,  a  wave  function 
analysis  of  Zip  shows  delocalization,  which  is  character¬ 
istic  of  a  shallow  level.  In  short,  although  the  (0/  +  ) 
transition  energy  of  Zip-No  cannot  be  accurately  calcu¬ 
lated,  it  is  entirely  consistent  with  the  observed  value  of 
30  meV.  Further  details  of  these  calculations  will  be  pub¬ 
lished  elsewhere. 

It  is  also  interesting  to  compare  the  0.9  eV  binding 
energy  with  the  activation  energy  for  irradiation-defect 
annealing  in  ZnO,  measured  previously  as  1.7  eV  [34]. 
The  difference  between  these  two  energies,  0.8  eV,  should 
be  the  motional  energy  for  Zrp  and  is  a  reasonable  value  for 
an  interstitial. 

Other  evidence  for  a  Zrp-N()  complex  comes  from  opti¬ 
cally  detected  magnetic  resonance  experiments  in  epitaxial 
N-doped  ZnO  [35].  An  observed  spin- 1/2  center  was  con¬ 
sistent  with  a  Zn  interstitial,  possibly  in  association  with  an 
N  atom,  since  it  was  not  observed  in  a  sample  with  lower  N 
content.  Besides  Zn7,  other  donors,  such  as  H,  are  also 
believed  to  associate  with  N  [14,36].  However,  H-N  would 
be  neutral,  whereas  Zn7-N  has  a  shallow-donor  level. 

Finally,  other  Zn7-acceptor  complexes  that  behave  as 
shallow  donors  are  also  predicted  to  be  stable.  For  ex¬ 
ample,  very  recently  Wardle  et  al.  [37]  have  used  DFT  to 
show  that  Zn7-LiZn  is  bound  by  0.7  eV  and  has  a  shallow- 
donor  transition.  This  result  further  strengthens  and  gen¬ 
eralizes  our  assertion  that  Zn7-related  shallow  donors  can 
exist  and  be  important  in  as-grown  ZnO. 

In  summary,  we  have  carried  out  extensive  low- 
temperature  photoluminescence  and  temperature- 
dependent  Hall-effect  measurements  on  irradiated  ZnO 
and  have  identified  a  defect-related  donor  at  about 
30  meV,  which  also  exists  in  as-grown  ZnO.  We  have 
further  carried  out  molecular  dynamics  simulations,  to 
show  that  the  expected  production  rate  of  Zn7  is  consistent 
with  the  concentrations  of  the  30  meV  donor  after  irradia¬ 
tion,  and  density  functional  calculations  to  show  that  the 
Zn rNQ  defect  complex  is  a  shallow  donor  with  a  sufficient 
binding  energy  to  explain  the  annealing  data.  Thus,  the 
conclusion  is  that  native-defect-related  donors  can  exist  in 
//-type  ZnO  and  contribute  to  its  conductance. 
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Abstract 

Based  on  first  principles  calculations,  we  investigate  two  probable  types  of  deactivation  mechanisms  that  hinder  current  efforts  of 
doping  ZnO  p-type.  (i)  Passivation  by  Hydrogen.  H  prefers  to  bind  with  N0  at  the  anti-bonding  site  and  form  N0-H  complexes  with  a 
binding  energy  of  about  1  eV.  (ii)  Passivation  by  the  formation  of  substitutional  diatomic  molecules  (SDM).  Carbon  impurities  and 
excess  N  strongly  prefer  to  passivate  N0  and  form  low-energy  SDM  on  the  Oxygen  site,  (NC)0  or  (N2)o,  both  of  which  are  donors  with 
several-eV  binding  energies.  Our  calculated  vibrational  frequencies  of  N0-H  complexes  and  SDMs  are  consistent  with  the  frequencies 
recently  observed  by  IR  measurement  on  N-doped  ZnO,  which  is  not  p-type. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Recent  SIMs  results  on  nitrogen-doped  ZnO  grown  by 
chemical  vapor  deposition  technique  show  a  substantial 
amount  of  H,  C  and  N  in  unwanted  forms  [1-4].  These 
impurities  are  likely  to  play  important  roles  in  compensat¬ 
ing  or  passivating  N  acceptors,  making  it  difficult  to  dope 
ZnO  p-type.  Recent  IR  measurements  on  N-doped  ZnO, 
which  is  not  p-type,  revealed  several  peaks  associated  with 
H,  including  a  peak  at  3020  cm-1  which  is  consistent  with 
the  N-H  stretch  vibration  frequency.  In  addition,  the 
spectrum  in  the  lower  frequency  range  contains  several 
peaks  in  the  range  of  1800-2000  cm-1,  which  are  close  to 
the  frequencies  of  diatomic  molecules. 
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Based  on  first  principles  calculations,  we  found  two  types 
of  defect  complexes  that  are  likely  to  form.  While  isolated 
H  is  predicted  to  exclusively  be  a  donor  in  ZnO  [5,6],  it  can 
also  passivate  N0  acceptors  and  form  an  electrically 
inactive  N0-H  complex  (Fig.  1)  [3].  For  the  second  type 
of  complex,  C  impurities  and  excess  N  atoms  prefer  to 
passivate  N0  acceptors  and  form  a  “substitutional  diatomic 
molecule”  (SDM),  which  is  a  strongly  bound  NC  or  N2 
molecule  substituted  on  an  oxygen  site  [(NC)0  or  (N2)o] 
(Fig.  2)  [7,8],  The  (NC)0  or  (N2)Q  SDMs  are  single  and 
double  donors,  respectively.  These  SDMs  have  electronic  and 
structural  properties  similar  to  free  diatomic  molecules.  The 
calculated  local  vibrational  modes  (LVM)  of  these  defects  are 
in  good  agreement  with  the  peaks  recently  observed  by  FTIR 
measurement  on  N-doped  ZnO.  These  results  suggest  that,  in 
addition  to  H,  one  should  be  careful  with  a  C  impurity  that 
can  turn  an  existing  N  acceptor  (N0)  into  an  (NC)0  SDM, 
which  is  a  donor.  In  addition,  excessive  N  incorporation 
could  lead  to  an  unwanted  (N2)0  SDM,  which  is  a  double 
donor,  instead  of  the  desired  N0. 
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Fig.  1.  Atomic  structure  of  the  N0-H  complex.  Large  sphere  is  Zn  and 
small  sphere  is  O.  Impurities  N  and  H  are  labeled. 


Fig.  2.  Atomic  structure  of  the  substitutional  diatomic  molecule,  where 
the  encircled  AB  dimer  replaces  an  O  atom  in  ZnO.  For  (Nolo  SDM, 
j4  =  1?  =  N.  For  (NC)o  SDM,  ^4  =  N  and  B  —  C.  The  large  sphere  is  Zn 
and  the  small  sphere  is  O. 


2.  Theoretical  approach 

2.1.  Total-energy  and  formation  energy  calculations 

Our  first  principles  calculations  are  based  on  the 
pseudopotential  planewave  method  with  density-functional 
theory  in  the  local-density  approximation.  The  Vanderbilt- 
type  ultra-soft  pseudopotential  [9],  as  implemented  in  the 
VASP  codes  [10],  is  used.  The  use  of  ultrasoft  versions  of 
the  pseudopotentials  allows  relatively  low  energy  cutoffs 
(300  eV)  to  be  used  for  the  planewave  basis  set.  The  Zn  3d 
electrons  are  treated  as  valence  electrons.  The  calculated 
heat  of  formation  for  ZnO  of  3.58  eV  agrees  with  the 
experimental  value  of  3.60  eV. 

To  study  defects,  we  use  a  supercell  approach  with  a 
supercell  containing  at  least  36  atoms.  Test  calculations  for 
the  (N2)o  and  N0-H  complexes  are  performed  with  a 
larger  96-atom  cell  to  ensure  that  the  convergence  is  within 
±0.2eV.  For  charged  defects,  a  jellium  background  was 
used.  Since  LDA  substantially  underestimates  the  ZnO 
band  gap,  we  have  instead  examined  the  electronic 
properties  at  the  2x2x2  Monkhorst-Pack  special  k 


points,  which  are  also  used  for  the  Brillouin  zone 
integration.  The  band  gap  at  the  special  k  points  is 
2.5  eV.  All  atoms  in  the  cell  are  allowed  to  relax  until  the 
magnitude  of  the  residue  force  on  each  atom  is  <0.05eV/A. 
The  defect  formation  energy  (All r)  which  has  been 
described  in  detail  elsewhere  (for  example,  Refs.  [11,12]), 
can  be  defined  as 

AH[  =  Etot(D,  q)  -  is  tot  (0)  +  A  «z„/tZn  +  An0g0 

+  Anxgx  +  </£’f,  0) 

where  Elol(D,q)  and  £tot(0)  are  the  total  energies  from  first 
principles  calculations  of  the  supercell  with  and  without  the 
defect,  /),  and  A(  =  H,  N,  C)  represents  impurity  elements. 
Quantities  A nA  and  gA  are  the  number  of  species  A  ( =  Zn, 
O,  H,  N,  C)  removed  from  a  defect-free  cell  to  its  respective 
reservoir  to  form  the  defect  cell  and  the  corresponding 
reservoir  chemical  potential.  The  chemical  potentials  /i/n, 
/to,  /(h,  13 n  and  /tc  have  upper  limits  at  their  respective 
natural  phases,  i.e.  the  energies  of  metallic  Zn,  gaseous  02, 
FI2,  N2  and  solid-state  C,  which  are  offset  to  zero  in  the 
present  study.  To  keep  the  ZnO  thermodynamically  stable, 
it  is  also  required  that  /tZn  +  /t0  =  /tzno.caic  =  -3.58  eV.  This 
imposes  an  additional  constraint  that  g0  be  in  the  range 
— 3.58 eV^/io^O  and  redefines  gZn  =  -3.58eV-/i0-  Unless 
noted  otherwise,  all  calculations  in  this  paper  are  done  in 
the  Zn-rich  limit,  i.e.  gZn  =  0.  In  Eq.  (1),  q  and  EF  are  the 
defect  charge  state  and  Fermi  level  with  respect  to  the 
VBM. 

2.2.  Vibration  frequency  calculations 

We  use  the  so-called  frozen  phonon  calculation  ap¬ 
proach.  After  the  defect  is  relaxed  to  the  lowest  energy 
configuration,  we  dislocate  the  vibrating  atom  (for  example 
FI  atom  in  the  N-H  FVM  case)  along  the  vibrating 
direction  in  a  small  step  at  a  time  (typical  step  is  0.05  A). 
For  each  step,  the  calculated  total  energy  is  recorded.  The 
potential  energy  curve  is  constructed  from  the  plot  of  total 
energy  versus  the  vibrating  distance. 

In  the  case  of  N0-H  complex,  it  is  a  good  approximation 
to  dislocate  only  the  H  atom  to  construct  the  potential 
energy  curve.  This  is  because  H  is  much  lighter  than  other 
atoms.  However,  when  calculating  the  N-H  stretch 
frequency  the  reduced  mass  of  N  and  H  atoms  must  be 
used.  Due  to  its  light  mass,  the  H  atom  vibrates  with  a 
large  amplitude,  i.e.  ~10%  of  the  bond  distance  in  the 
ground  state  and  20%  in  the  first  excited  state.  Therefore,  it 
is  important  to  include  anharmonic  effects  in  the  FVM 
calculation  [13,14], 

In  the  case  of  SDMs,  the  situation  is  somewhat  different. 
N  and  C  have  a  comparable  mass  that  is  over  ten  times 
heavier  than  H,  leading  to  a  much  smaller  zero-point 
vibration  amplitude.  Since  the  amplitude  of  the  vibration  is 
small,  the  anharmonicity  in  N-C  and  N-N  vibrations  are 
much  smaller  than  that  of  the  H-N  vibration.  For  example, 
including  anharmonicity  in  the  calculation  only  leads  to  a 
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change  in  the  LVM  of  a  free  N2  molecule  by  <  20  cm- 1  out 
of  a  total  frequency  of  about  2000  cm-1  [14].  In  calculating 
the  LVM  of  SDMs,  both  atoms  in  the  dimer  are 
simultaneously  dislocated  with  the  same  amplitude  in 
opposite  directions.  Other  atoms  in  the  supercell  are  fixed. 

3.  Results  and  discussions 

3.1.  Passivation  by  hydrogen 

Using  the  definition  of  the  formation  energy  given 
above,  we  calculated  the  formation  energies  of  an  isolated 
interstitial  H,  a  nitrogen  acceptor  N0,  and  a  N0-H 
complex  in  ZnO.  An  isolated  interstitial  H  is  exclusively 
a  donor  in  ZnO  and  always  exists  in  the  form  of  H+.  For 
an  isolated  H+  location,  there  are  four  low-energy  sites 
surrounding  an  O  atom  with  the  so-called  BCy  site  being 
the  site  with  the  lowest  energy  [6,15].  The  other  three  sites 
have  slightly  higher  energies  (approximately  within  0.2  eV). 
An  isolated  N0  is  an  acceptor.  Currently  the  exact  location 
of  the  ionization  energy  of  N0  is  still  under  debate  [3].  Over 
almost  the  entire  Fermi  energy  range,  the  isolated  N0  is 
stable  in  a  (1— )-charge  state  whereas  an  isolated  interstitial 
H  atom  is  always  stable  in  a  (1 +)-charge  state.  The  two 
defects,  therefore,  have  a  Coulomb  attraction  and  have  a 
strong  tendency  to  form  a  N0-FI  defect  complex.  The 
formation  energy  of  H+  (at  BCy  site)  and  N0  are  plotted 
(dashed  lines)  as  a  function  of  Fermi  level  in  Fig.  3  with  the 
sum  of  the  two  formation  energies  shown  as  dotted  lines.  In 
the  same  figure,  we  also  plot  the  formation  energy  of  the 
N0-H  complex  (thick  solid  line).  The  binding  energy  of  the 
N0-H  complex  is  0.95  eV,  i.e.  the  difference  in  energy 
between  the  dotted  and  solid  line.  For  the  N0-H  defect 
complex,  the  ABNX  configuration  (Fig.  1)  has  the  lowest 


Fig.  3.  Formation  energy  of  an  interstitial  FI,  N0  (dashed  lines)  and  a 
No-FI  complex  (solid  line)  in  ZnO  as  a  function  of  electron  Fermi  energy. 
The  dotted  line  shows  the  sum  of  the  formation  energies  of  an  isolated 
interstitial  hydrogen  and  an  isolated  No-  The  Zn-rich  condition.  N2  and 
H2  phase  precipitation  limits  were  assumed.  The  energy  difference  between 
the  dotted  line  and  the  solid  line  is  the  binding  energy  of  the  No-FI 
complex. 


formation  energy.  The  ABN  configurations  are  energeti¬ 
cally  more  stable  than  the  BC  by  about  0.2-0. 3  eV,  making 
the  likelihood  of  forming  BC  configurations  negligible.  The 
difference  in  formation  energy  between  the  two  ABN 
configurations  is  within  the  computational  error  bar.  Both 
ABn  configurations  have  a  similar  N-H  stretch  mode 
frequency,  with  an  average  value  of  2927  cm-1.  The 
anharmonic  contributions  are  rather  large,  i.e.  approxi¬ 
mately  8%.  Both  BC  configurations  also  have  a  similar 
N-FI  stretch  mode  frequency,  with  an  average  value  of 
3319  cm-1,  which  is  clearly  higher  than  those  of  the  ABN 
configurations.  The  anharmonic  contributions  are  less  than 
half  of  those  for  the  ABN  cases.  The  smaller  anharmonic 
contribution  in  BC  in  comparison  with  ABN  agrees  with 
the  usual  trend  found  for  an  isolated  FI  in  ZnO  [6]  as  well 
as  in  the  GaN  system  [13],  To  estimate  the  systematic  error 
of  our  calculation,  we  calculated  the  known  stretch-mode 
frequency  of  a  free  ammonia  (NH3)  molecule.  We  obtained 
the  symmetric  stretch  LVM  of  3194  cm-1  which  is 
143  cm-1  smaller  than  the  known  experimental  value 
(3337  cm-1).  Assuming  that  our  calculated  LVM  of  the 
N-FI  bond  in  ZnO  is  underestimated  by  similar  amount, 
our  adjusted  value  for  ABN  of  2927  +  143  =  3070  cm-1  is  in 
reasonable  agreement  with  the  recently  observed  peak  at 
3020  cm-1  [3], 

3.2.  Passivation  by  the  formation  of  substitutional  diatomic 
molecules 

The  diatomic  molecule  is  one  of  the  most  common  forms 
for  first  row  elements,  namely,  C,  N,  and  O,  because  these 
small  atoms  tend  to  form  strong  triple  bonds  with  each 
other.  Recently,  similar  complexes,  such  as  N2  and  NC 
inclusions  in  GaAs  have  been  studied  both  theoretically 
and  experimentally  [16,17].  In  the  case  of  ZnO,  the  (N2)Q 
and  (NC)0  SDM  are  found  to  have  bonding  and  anti¬ 
bonding  characteristics  very  similar  to  their  free  molecules 
confirmed  by  the  detailed  investigation  of  the  electronic 
wavefunctions  [7,8].  Flere,  we  will  focus  only  on  the 
energetic  and  the  LVM  of  the  complexes. 

Fig.  4  shows  the  EF  dependence  of  the  formation  energy 
(AH{  )  in  the  Zn-rich,  N-rich,  and  C-rich  conditions 
(dzn  =  A<n  =  He  =  0)  for  (CO)o,  (NC)o,  (NO)o,  and  (N2)o- 
Since  the  isolated  interstitials  Nj  and  C;  have  very  high 
formation  energies,  they  tend  to  immediately  bind  with  the 
closest  lattice  O,  forming  (NO)0  and  (CO)0-  However,  in 
the  presence  of  N0,  (NO)o  and  (CO)o  are  unstable  against 
the  formation  of  (N2)Q  and  (NC)0-  The  binding  energies 
are  typically  several  eV  due  in  part  to  Coulomb  attraction 
between  oppositely  charged  impurities  and  in  part  to  the 
combination  of  two  impurities  into  one.  This  can 
significantly  reduce  the  formation  energy  because,  for 
example,  when  N0  +  (NO)0  -»■  (N2)Q  +  00,  one  less  oxy¬ 
gen  site  is  disturbed  by  the  impurities  after  the  reaction.  In 
the  right  panel  of  Fig.  4,  the  formation  energy  of  N0  and 
(NO)o  are  shown  as  the  dashed  line  and  the  thin  solid  line 
with  the  sum  of  the  two  as  the  dotted  line.  The  sum  is  to  be 
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Fig.  4.  Formation  of  (left  panel)  (NC)o  from  isolated  N0  and  (CO)o  and 
(right  panel)  (N2)o  from  N0  and  (NO)o-  The  Zn-rich  condition,  N2  and 
diamond  phase  precipitation  limits  were  assumed. 


compared  with  the  thick  solid  line,  which  is  the  formation 
energy  of  (N2)0  defect.  We  can  immediately  see  that  the 
binding  energy,  which  is  the  difference  between  the  dotted 
line  and  the  thick  solid  line,  depends  somewhat  on  the 
Fermi  energy,  but  has  a  typical  value  of  several  eV.  The 
binding  energy  is  of  similar  magnitude  for  the  case  of 
(NC)o,  as  shown  in  the  left  panel.  This  is  especially  true  in 
the  case  of  p- type  ZnO,  for  which  the  binding  energies  are 
large  enough  to  make  both  (N2)o  and  (NC)0  lower  in 
energy  than  isolated  N0.  It  is  therefore  expected  that  the 
formation  of  (N2)0  and  (NC)0  will  compete  with  the 
formation  of  N0  acceptors.  Since  these  SDMs  are  donors 
in  p-type  samples,  their  formation  further  hinders  p-type 
doping  by  compensating  the  already  lowered  N0  acceptors. 

The  situation  could  be  worse  when  one  tries  to  increase 
N  concentration,  [N],  above  the  equilibrium  value  via  non¬ 
equilibrium  growth  methods.  Effectively,  this  is  similar  to 
increasing  the  nitrogen  chemical  potential  //N.  Although 
the  formation  energy  of  N0  decreases  as  one  raises  /(N,  the 
formation  energy  of  the  (N2)Q  complex  decreases  twice  as 
fast.  In  thermal  equilibrium,  this  would  end  up  increasing 
the  concentration  of  the  donor  (N2)G  complex  instead  of 
the  acceptor  N0.  However,  because  the  formation  of  (N2)Q 
from  No  +  (N0)o  involves  the  breakup  of  one  SDM,  the 
diffusion  of  N„  and  the  formation  of  another  SDM,  one 
might  be  able  to  adjust  the  growth  temperature  to  be  low 
enough  that  (N2)Q  is  suppressed. 

The  calculated  frequencies  for  (N2)o+  and  (NC)o,  which 
are  the  most  stable  form  for  these  complexes  in  p-type 
ZnO,  are  2108  and  1995  cm-1,  respectively.  These  frequen¬ 
cies  fall  within  200  cm-1  of  the  experimentally  observed 
values  [2].  Moreover,  recent  XPS  experiments  have  found 
strong  signals  of  N-N  and  C-N  bonds  with  similar 
characteristics  to  those  of  free  diatomic  molecules  [4]. 
These  observations  strongly  suggest  the  existence  of  (N2)Q 
and  (NC)0  in  these  samples. 


4.  Conclusion 

First  principles  total  energy  calculations  show  that  H,  C, 
and  N  can  passivate  N0  acceptors  in  ZnO.  H  can  passivate 
N0  and  form  N0-H  complexes  with  the  binding  energy  of 
approximately  1  eV.  C  and  excess  N  can  passivate  N0  and 
form  low  energy  (NC)0  and  (N2)0  substitutional  diatomic 
molecules  (SDM)  on  the  oxygen  site  in  ZnO.  The 
calculated  vibrational  frequencies  are  in  qualitative  agree¬ 
ment  with  recent  IR  experiments.  Our  calculations  also 
show  that  both  (NC)0  and  (N2)0  are  donors  in  p-type 
ZnO.  These  results  should  shed  new  light  on  efforts  to 
improve  the  fabrication  of  high  quality  p-type  ZnO  by 
nitrogen. 
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Abstract 

First-principles  calculations  show  that  the  hydrogen  configurations  in  GaAsN  depend  on  how  hydrogen  is  introduced  into  the  sample. 
Since  proton  and  neutral  H  have  different  ground  states,  the  proton  injected  into  the  sample  by  H-irradiation  follows  a  unique  energy 
pathway  to  form  a  charged  dihydride,  instead  of  the  charge-neutral  H2  monohydride.  The  subsequent  charge  neutralization  causes  the 
spontaneous  canting  of  the  dihydrides.  The  resulting  canted  N-2H  structure  explains  the  recent  puzzling  1R  observation,  the  recoveries 
of  the  GaAs  band  gap  and  lattice  parameter,  and  the  dihydride  symmetry  determined  by  the  XANES  experiment.  It  may  also  have  broad 
implications  for  ion  implantation  studies  in  other  solids. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Large  size-mismatched  dilute  alloys  such  as  GaAsN  [1] 
and  ZnSeO  [2,3]  are  promising  alloy  semiconductors  with 
unique  physical  properties.  For  example,  with  only  an 
atomic  percent  of  N  incorporated  into  GaAs,  the  band  gap 
decreases  by  several  hundred  meV  [4].  On  the  other  hand, 
H  irradiation  causes  a  nearly  full  recovery  of  the  GaAs 
band  gap  [5]  and  lattice  parameter  [5]  from  those  of 
GaAsN.  The  formed  H  complexes  after  H  irradiation 
should  be  charge  neutral  due  to  their  large  concentration 
( >  102"cm-3).  First-principles  total  energy  calculations 
suggested  that  the  lowest-energy  charge-neutral  H  complex 
is  the  a-FF(N)  complex  with  characteristic  Ga-H  and  N-H 
monohydride  bonds  (Fig.  la),  and  that  this  complex  can 
explain  the  above  experimental  observations  [6,7].  How- 
ever,  a  recent  infrared  experiment  [8]  on  H-irradiated 
samples  did  not  observe  any  Ga-H  modes.  Instead,  two 
N-H  stretch  modes  were  observed  at  3195  and  2967  cm-1, 
as  shown  in  Table  1.  Hence,  both  modes  here  belong  to  the 
N-H  bonds.  In  D-irradiated  samples,  the  two  stretch 
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modes  shift  to  2216  and  2376cm-1,  respectively.  In  H/D 
co-irradiated  samples,  however,  two  additional  D-stretch 
modes,  at  2221  and  2366  cm-1  (see  Fig.  2a),  and  one 
additional  H-stretch  mode,  were  observed.  A  second  H 
mode  is  also  expected,  but  was  not  observed,  which  was 
due  possibly  to  contamination.  The  observation  of  the  two 
(instead  of  one)  additional  D  modes  indicates  that  the  H 
complexes  must  contain  two  inequivalent  H  atoms  strongly 
coupled  to  each  other.  A  recent  XANES  study  [9]  further 
reveals  that  the  symmetry  of  the  2  H  complexes  should  be 
C2v-like,  thus  further  undermining  the  validity  of  the  C3v- 
based  H2  model. 

The  disagreement  between  theory  and  experiment  on  H 
configurations  in  GaAsN  suggests  that  the  observed  H- 
complex  may  not  be  the  ground-state  configuration. 
During  H-irradiation,  the  equilibrium  between  the  injected 
protons  and  the  electron  reservoir  of  energy  EF  is  not 
maintained.  Therefore,  one  cannot  assume  that  protons 
will  take  the  theoretically  predicted  ground-state  structure 
for  neutral  H. 

In  contrast  to  the  conventional  view  of  H  irradiation  as 
simply  another  convenient  mean  to  enhance  H  concentra¬ 
tion  [H]  in  the  solid,  in  this  paper,  we  will  show  that  the 
final  H  configurations  can  be  altered  by  how  the  H  is 
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introduced  into  the  solid.  For  H  irradiation  of  GaAsN,  the 
relaxation  pathway  for  proton  must  be  considered  instead 
of  the  neutral  H,  because  H  irradiation  in  large  quantity 
causes  the  accumulation  of  positive  charges  that  cannot  be 
easily  dissipated.  As  a  result,  the  injected  protons  enter  the 
sample  and  get  trapped  in  the  proton-ground-state 
structure,  N-2Hbc,  rather  than  the  neutral-H-ground-state 
structure,  a-H2.  The  subsequent  charge  neutralization 
results  in  the  formation  of  metastable  canted  N-2H 
complex,  which  cannot  be  converted  easily  to  the  ground- 
state  structure,  a-H2(N),  due  to  a  kinetic  barrier.  The 
canted  N-2H  model  explains  all  experimental  observa¬ 
tions,  including  IR  modes  (both  in  frequency  and  relative 
amplitude),  dihydride  symmetry,  and  the  bandgap  and 
lattice  parameter  recovery.  The  good  agreement  between 
theory  and  experiments  establishes  a  new  direction  to 
examine  ion  implantations  in  semiconductors  and  solids  in 
general. 


2.  Method 

Our  calculations  are  based  on  the  density  functional 
theory  within  the  local  density  approximation  (LDA).  The 
electron-ion  interactions  are  described  by  the  ultra-soft 
pseudopotentials  [10],  as  implemented  in  the  VASP  codes 
[11].  The  valence  wavefunctions  are  expanded  in  a  plane- 
wave  basis  with  a  cutoff  energy  of  348  eV.  All  calculations 
were  performed  using  a  64-atom  supercell.  The  vibrational 
frequencies  were  calculated  based  on  the  dynamical  matrix 
approach  [12].  We  used  the  special  k-points  scheme  for 
Brillouin  zone  summation  and  relaxed  all  the  atoms  until 
the  forces  were  less  than  0.02  eV/A. 


3.  Results 

Once  protons  enter  the  sample,  they  are  unlikely  to  be 
neutralized,  because  the  ( +  /0)  level  for  the  diffusing  H  is 
near  the  conduction  band  minimum  of  GaAs,  which  is  well 
above  that  of  GaAsN.  It  therefore  makes  sense  to  study  the 
H+  instead  of  H°  traps  in  H-irradiated  GaAsN. 

The  most  likely  proton  traps  in  GaAsN  are  the  N- 
centered  di-proton  complexes,  among  which  we  found  that 
neither  the  oc-H22+  nor  the  interstitial  H2  +  (Ga)  (namely,  a 
charged  molecule  at  the  tetrahedral  interstitial  site  next  to  a 
Ga)  is  a  good  candidate  for  the  proton  trap,  because  a 
positive  charge  completely  destabilizes  the  respective 
atomic  structure.  Instead,  the  N-2H|c  complex  is  found 
to  be  the  ground-state  for  a  di-proton  complex  (see  Fig. 
lb).  Charge  neutral  a-H2  and  interstitial  H2(Ga)  may 
occasionally  form  by  simultaneously  creating  two  holes 
near  the  valence  band  maximum.  Even  if  this  happens, 
however,  our  calculation  shows  that  the  resulting  i- 
H2°(N)  +  2 h+  and  H2jint  (Ga)  +  2 h  +  +  NAs  are  significantly 
less  stable  than  the  N-2Hgc  (by  0.47  and  0.80  eV/H, 
respectively).  On  the  other  hand,  the  trapping  of  two  H+  to 
form  N-2Fl|c  lowers  the  energy  by  0.48eV/FI  [13].  Each 
trapped  H  forms  a  strong  bond  with  the  N  (see  Fig.  lb)  by 
passing  the  positive  charge  to  a  nearby  threefold-coordi¬ 
nated  Ga  in  the  N-2Fl|c  complex.  This  could  be  an 
important  initial  step  to  separate  holes  from  the  dihydrides 
for  charge  neutralization. 
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Fig  2.  The  various  D  stretch-modes  for  (a)  experiment,  (b)  N-2H2  +  ,  and 
(c)  canted  N-2H.  Solid  lines  indicate  D  stretch-modes  in  fully  deuterated 
samples,  whereas  dashed  lines  indicate  additional  D  stretch-modes  in  H/D 
flux. 


Table  1 

Stretch  and  wag  mode  frequencies  in  units  of  cm-1.  N/A  stands  for  not  available 


Experiment 

ct-Hj 

N-2H2+ 

Canted  N-2H 

Stretch 

wag 

Stretch 

wag 

Stretch 

wag 

Stretch 

wag 

H-modes  (H-flux) 

3195 

2967 

1447 

3140 

944 

3305 

3209 

1328 

3207 

3052 

1417 

D-modes  (D-flux) 

2376 

2216 

1076 

2290 

687 

2431 

2318 

974 

2346 

2215 

1037 

Additional  Modes  (H/D  flux) 

3192 

N/A 

N/A 

None 

3260 

1164 

3204 

3056 

1257 

2366 

2221 

N/A 

2372 

1164 

2336 

2224 

1238 
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Bonapasta  et  al.  proposed  that  the  N-2Hbc  complex 
accounts  for  the  measured  H-modes  [14],  However,  a 
sample  cannot  host  a  few  atomic  percent  of  positively 
charged  H  complexes  and  also  the  N-2H|c  complex  does 
not  agree  with  experimental  observation  of  two  inequi¬ 
valent  N  H  bonds.  Thus,  the  N-2H|J  complex  cannot  be 
the  final  H  configuration  after  H  irradiation.  The  charge 
neutralization  of  the  N-2H|£  must  take  place  as  the 
sample  is  grounded,  but  at  a  later  time.  This  could  happen, 
for  example,  by  capturing  electrons  (injected  from  ground 
by  the  electrostatic  potential  between  sample  and  ground) 
in  the  conduction  band. 

Upon  the  N-2H|</  neutralization,  two  important  things 
happen  spontaneously  to  N-2HIJ:  (1)  the  formation  of  the 
Ga-Ga  bond  in  Fig.  3a  and  (2)  the  canting  of  the  dihydride 
in  Fig.  3b.  The  former  eliminates  electronic  states  from 
within  the  original  GaAs  band  gap  in  agreement  with  the 
gap  recovery  experiment.  The  latter  qualitatively  alters  the 
H/D  vibrational  properties.  Because  of  its  importance,  we 
have  calculated  not  only  the  initial  and  final  positions  of 
canting,  but  also  four  intermediate  ones  along  the  pathway. 
The  smoothness  of  the  curve  in  Fig.  3  suggests  that  the 
results  are  numerically  highly  accurate.  We  note  that  such  a 
canting  occurs  only  when  the  weak  Ga-Ga  bond  forms  due 
to  the  Coulomb  attraction  between  the  electron-rich 
Ga-Ga  bond  and  the  electron-depleted  H  bonded  to  N. 

Table  1  tabulates  the  calculated  H  and  D  frequencies. 
When  compared  with  experiment,  the  dihydride  N-2Hbc 
and  canted  N-2H  models  produce  better  results  for  the  H 
wag  mode,  1328  and  1417  cm-1,  respectively,  than  the 
monohydride  oc-H2  model  does,  944  cm-1.  This  is  another 
indication  that  the  defect  responsible  for  the  IR  signals 
cannot  be  a  monohydride.  Because  the  two  additional  D 
modes  in  the  mixed  H/D  flux  are  key  IR  results,  they  are 
examined  in  greater  details  in  Fig.  2:  The  canting  of  the 
dihydride  leads  to  two  additional  modes  versus  one 
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Fig  3.  Calculated  canting  pathway,  indicating  that  the  uncanted  N-2H  is 
unstable  against  the  canting.  The  structures  of  uncanted  and  canted  N-2H 
complexes  are  shown  in  insets  (a)  and  (b),  respectively. 


additional  mode  of  the  uncanted  N-2H|J.  Only  the 
former  agrees  with  the  IR  experiment.  The  canting  of  the 
dihydride  also  results  in  two  equivalent  minima  in  the 
energy  landscape,  which  become  vibrationally  different 
when  one  H  is  replaced  by  D.  The  barrier  in  Fig.  3  is  only 
50meV.  Hence,  we  expect  that,  at  low  T.  each  N-DH 
complex  should  be  trapped  to  one  of  the  minima  and 
contributes  only  to  one  of  the  additional  D  modes.  At 
higher  T,  however,  there  will  be  no  such  trapping  so  that 
each  N-DH  complex  contributes  to  both  additional  D 
modes,  possibly  in  a  1:1  ratio. 

We  have  also  calculated  the  relative  IR  intensities  (/)  for 
the  D  modes  in  Fig.  2c.  The  ratio  of  IR  intensities 
originated  from  the  same  defect  can  be  approximated  to 
the  first-order  by  the  ratio  of  |dp/d(9,j  [2,15,16],  where  p  is 
the  dipole  moment  of  the  system  and  Q,  is  the  rth  normal¬ 
mode  coordinate.  We  calculated  the  dipole  moment  change 
Ap  for  each  mode  after  taking  same  distance  (\AQ\)  of 
atomic  displacements  along  each  normal-mode  coordinate 
(Q,:).  Then  the  ratio  of  any  two  IR  intensities  here  is  simply 
the  ratio  of  |Ap|2  resulted  from  atomic  displacement  of 
same  \AQ\  in  two  normal-mode  coordinates.  This  yields 
^2215^2224^2336^2346  =  1:0.84:0.44:0.21,  where  the  sub¬ 
scripts  are  the  mode  frequencies  (cm-1).  These  results  can 
be  compared  to  the  measured  IR  intensities,  estimated 
from  reading  Fig.  2  in  Ref.  [8],  hne-hin ^2366^2376 
« 1:0.43:0.33:0.29.  By  reproducing  the  qualitative  experi¬ 
mental  trend,  these  calculations  provide  further  support  to 
the  canted  N-2H  model. 

Once  the  H  atoms  are  trapped  in  a  neutral  canted  N-2H 
complex,  they  may  not  escape  easily  to  form  a-H2,  despite 
that  the  latter  is  0.2eV/H  lower  in  energy.  The  minimum- 
energy  barrier  for  converting  a  canted  N-2H  to  a-H2  is 
calculated  to  be  0.8  eV.  In  addition,  because  of  the  large 
mass  differences  between  H  and  Ga  and  between  H  and  N, 
the  actual  barrier  to  convert  a  canted  N-2H  to  a-H2*  could 
be  significantly  higher  than  the  calculated  0.8  eV  [17].  Note 
that  the  kinetic  energy  carried  by  proton  should  not  help 
the  N-2H  overcome  the  energy  barrier  to  become  a-H2, 
because  it  dissipates  quickly  after  protons  enter  the  sample 
by  a  distance  of  only  a  few  nm  [18],  which  is  two  orders  of 
magnitudes  smaller  than  GaAsN  film  thickness  of  300  nm 
studied  in  Ref.  [8].  Also,  a  GaAs  capping  layer  of  100-nm 
thickness  is  usually  applied  to  eliminate  direct  damages  by 
the  energetic  protons  [19,20].  Hence,  protons  enter  most  of 
GaAsN  layer  by  thermal  diffusion.  Thus,  one  can  safely 
ignore  any  effect  of  energetic  protons  and  local  heating  on 
the  relaxation  of  the  metastable  H  complexes. 

Despite  the  kinetic  barrier,  there  is  a  finite  probability 
that  a-H2(N)  may  form.  However,  even  if  this  happens,  the 
H  at  a-H2(N)  can  be  etched  away  by  diffusing  protons.  Our 
calculations  show  that,  an  H  +  (As),  which  has  roughly  a 
0.2-eV  diffusion  barrier  in  GaAsN,  can  spontaneously 
transform  any  a-H2(N)  into  an  interstitial  H2(Ga)  plus  an 
H  +  (N).  This  reaction  is  not  only  exothermic  with  energy 
lowering  by  0.61  and  0.48  eV  for  H°  and  H  +  ,  respectively, 
but  also  barrierless  other  than  the  diffusion  one.  These, 
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combined  with  the  proposed  proton  relaxation  pathway, 
explain  why  a  metastable  dihydride  C2v  symmetry  was 
observed  by  the  recent  XANES  study,  instead  of  the  global 
minimum-energy  monohydride  C3v  symmetry. 

Finally,  we  studied  the  recovery  of  the  GaAs  lattice 
parameter  upon  H  irradiation.  To  mimic  the  epitaxial 
growth,  we  fixed  the  planar  lattice  parameters  to  that  of 
bulk  GaAs,  but  allowed  the  third  one  along  [0  0  1]  to  relax. 
Assuming  the  lattice  parameter  of  bulk  GaAs  is  100%,  we 
obtain  98.69%  for  NAs,  98.81%  for  HbC(N),  99.59%  for  a- 
H2(N),  and  100.23%  for  H2(Ga).  The  canted  N-2H  may 
take  two  distinct  orientations  with  respect  to  the  [001] 
direction  with  degeneracy  of  2  and  4,  respectively.  The 
corresponding  lattice  parameters  are  100.55%  and  99.37%, 
so  the  weighted  average  is  99.76%.  This  result  is  slightly 
better  than  that  of  a-H2,  but  could  be  expected,  because  for 
H2,  two  H  atoms  cut  only  one  strained  N-Ga  bond, 
whereas  for  canted  N-2H,  the  two  H  atoms  cut  two  N-Ga 
bonds,  instead,  leaving  more  rooms  for  the  lattice  recovery. 

4.  Summary 

In  summary,  we  have  developed  a  model  for  ion 
relaxation  following  the  initial  implantation.  The  theory 
is  qualitatively  different  from  the  ordinary  equilibrium 
defect  theory  and  may  have  implications  for  a  broad  range 
of  ion  implantation  experiments.  Application  of  this  theory 
to  GaAsN:H  reveals  a  canted  N-2H  complex,  which 
simultaneously  accounts  for  all  major  experimental  ob¬ 
servations  including  the  puzzling  IR  frequencies  and 
amplitudes,  the  dihydride  symmetry  by  XANES,  the 
recovery  of  the  GaAs  band  gap  by  optical  measurement, 
and  the  recovery  of  the  GaAs  lattice  parameter  by  X-ray 
diffraction. 
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Abstract 

We  studied  the  role  of  impurities  in  nitrogen-doped  ZnO  thin  film  to  understand  the  difficulty  of  producing  p-type  ZnO  via  nitrogen 
doping.  The  ZnO:N  films  were  fabricated  by  low-pressure  metal-organic  chemical  vapor  deposition  (MOCVD)  using  diethylzinc  (DEZ) 
and  nitric  oxide  (NO)  precursors.  Although  very  high  levels  of  nitrogen  incorporation  were  observed  ('"TO21  cm-3),  acceptor 
concentrations  were  typically  low  (1014-1017cm-3).  The  investigation  suggests  that  the  low  carrier  concentrations  are  possibly  due  to 
compensation  and  passivation  effects  by  hydrogen  and  carbon  impurities  unintentionally  incorporated  into  the  films  from  the  metal- 
organic  precursor.  X-ray  photoelectron  spectroscopy  (XPS)  demonstrated  that  carbon  was  a  bulk  impurity  in  MOCVD-grown  films. 
Secondary-ion  mass  spectrometry  (SIMS)  analysis  confirmed  the  presence  of  carbon  and  indicated  that  hydrogen  was  also  a  bulk 
impurity.  The  concentration  of  carbon  contaminant  was  found  to  increase  with  nitrogen  doping.  Both  XPS  and  Fourier  transform 
infrared  spectroscopy  (FT1R)  data  indicated  that  defect  complexes  (CHV,  NHr,  and  NCT)  are  likely  present  in  MOCVD-grown  ZnO 
films.  First-principles  calculations  predict  that  the  N0-H  and  (NC)o  defect  complexes  are  neutral  and  1  +  charge  state;  therefore,  the 
existing  carbon  and  hydrogen  passivate  the  nitrogen  acceptor  species.  Thus,  we  believe  a  low  hole  concentration  in  MOCVD-fabricated 
ZnO:N  films  are  partially  due  to  inadvertent  passivation  by  hydrogen  and  carbon. 

©  2005  Elsevier  B.V.  All  rights  reserved. 

PACS:  71.55.Gs  impurity  and  defect  levels  in  II-VI  semiconductors 
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1.  Introduction 

ZnO  is  a  wide-band  gap  (~3.37  eV)  semiconductor  that 
has  been  used  extensively  in  surface  acoustic  wave  devices 
and  solar  cells  [1,2].  It  has  typically  been  used  as  an  n-type 
semiconductor  because  it  is  straightforward  to  prepare  the 
material  as  n-type  with  or  without  extrinsic  doping.  In 
recent  years,  the  international  research  community  has 
devoted  significant  effort  to  fabricate  p-ZnO  because 
manufacturable  p-type  material  would  allow  applications 
such  as  ultraviolet  (UV)  light-emitting  diodes  and  lasers. 
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To  date,  there  is  substantial  evidence  that  p-type  ZnO  has 
been  achieved  by  a  number  of  groups  [3-10].  However,  it 
has  typically  been  difficult  to  prepare  high-quality  p-type 
ZnO  material  reproducibly;  the  p-type  ZnO  is  usually 
characterized  by  relatively  low-hole  concentration,  low- 
hole  mobility,  and  instability  [3,6,7,9,11].  Especially  under 
the  ultraviolet  photoelectron  illumination,  the  p-type  ZnO 
film  appears  as  n-type  [12], 

Efforts  to  make  p-ZnO  thin  films  have  involved  several 
deposition  methods  (e.g.,  CVD,  sputtering,  spray  deposi¬ 
tion,  diffusion,  oxidation  of  Zn3N2)  and  several  different 
dopants  (e.g.,  nitrogen,  phosphorus,  arsenic,  antimony, 
and  lithium).  A  carrier  concentration  as  high  as 
1.7  x  1019cm-3  has  been  reported  for  phosphorus-doped 
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p-type  ZnO  [13].  In  the  case  of  nitrogen,  co-doping  with 
aluminum,  gallium  and  indium  has  been  reported  with 
varied  degrees  of  success  [4,14—16].  In  many  cases,  the 
window  of  the  deposition  parameters  for  achieving  p-type 
material  was  narrow  [6,14,17].  Post-deposition  annealing 
was  used  in  some  cases  to  realize  p-type  conduction  [3]. 

Nitrogen  has  been  regarded  as  the  most-promising 
acceptor  dopant  because  of  its  low  ionization  energy, 
suitable  ionic  radius,  ease  of  handling,  low  material 
toxicity,  and  abundance.  Much  effort  has  been  devoted 
to  nitrogen-related  doping  and  co-doping,  but  low- 
resistivity  p-type  ZnO:N  has  yet  to  be  obtained.  For 
ZnO:N  film  grown  by  CVD,  the  typical  carrier  concentra¬ 
tions  are  only  1015  1018,  despite  nitrogen  concentrations  in 
the  range  of  1021cm-3  [6,18],  The  reasons  for  the  low 
electrical  activity  of  incorporated  nitrogen  are  still  under 
debate,  but  it  is  well  known  that  films  grown  via  MOCVD 
can  contain  hydrocarbon  impurities  that  originate  from  the 
metal-organic  precursors  and  that  such  impurities  can 
affect  the  action  of  dopants  such  as  nitrogen.  In  this  paper, 
we  report  on  the  results  of  an  in-depth  study  of  carbon  and 
hydrogen  impurities  in  ZnO:N  thin  films  fabricated  by 
MOCVD  and  discuss  the  possibility  of  the  compensation 
of  the  p-type  nitrogen  doping  by  various  hydrogen-  and 
carbon-related  defect  complexes. 

2.  Methods 

Nitrogen-doped  ZnO  thin  films  were  fabricated  by  low- 
pressure  MOCVD.  The  precursors  were  DEZ  (C2H5)2Zn), 
oxygen,  and  dilute  NO  gas  (2  wt%  NO/Ar).  The  carrier  gas 
for  all  the  precursors  was  N2  gas.  The  deposition 
temperature  was  varied  between  200  and  550  °C.  Undoped 
ZnO  samples  were  grown  using  the  reaction  between  DEZ 
and  02  gas.  The  nitrogen-doped  ZnO  samples  were  grown 
using  DEZ  and  NO.  Two  types  of  substrates,  Si  and 
Corning  1737  glass,  were  used.  The  Si  substrate  was 
suitable  for  subsequence  FTIR  analysis,  whereas  the  glass 
substrate  was  suitable  for  other  characterization  techniques 
such  as  X-ray  diffraction  (XRD),  XPS,  SIMS,  and  optical 
and  electronic  property  analyses. 

Atomic  force  microscopy  (AFM)  (Dimension  3100  with 
Nanoscope  Ilia  controller,  using  Si  cantilevers)  was  used  to 
assess  the  topography  of  the  thin  films.  Crystal  properties 
were  assessed  using  XRD  (Sintag  Model  PTS).  A  Cary  5G 
spectrophotometer  with  an  integrating-sphere  detector  was 
used  to  measure  total  transmittance  and  reflectance 
spectrum  in  the  wavelength  ranging  from  350  to  2000  nm. 
Electrical  properties  were  studied  by  Hall  measurement 
(BioRad  Model  HL5500)  and  capacitance-voltage  (C-V) 
analysis  (double-Schottky  surface-contact  mercury  probe, 
MSI  Electronics  Hg-2C-5).  The  nitrogen,  hydrogen,  and 
carbon  concentrations  were  analyzed  using  SIMS  (Cameca 
IMS  5f)  and  XPS  (Physical  Electronics  Model  5600). 
Details  of  the  XPS  analysis  can  be  found  in  a  previous 
publication  [12].  The  N-H  bond  signatures  were  studied  by 
room-temperature  FTIR  transmittance  with  ZnO  films 


grown  on  parallelogram-shaped  single-crystal  silicon  sub¬ 
strates  and  mounted  in  an  attenuated  total-reflectance 
(ATR)  accessory  attached  to  a  Nicolet  Magna-IR  550 
spectrometer  [19]. 

We  performed  first-principles  calculations  based  on 
density  functional  theory  (DFT)  within  the  local  density 
approximation  (LDA)  and  ultrasoft  pseudopotentials;  As 
implemented  in  the  Vienna  ab-initio  simulation  package 
(VASP)  codes.  Zn  3 d  states  were  included  in  the  valence, 
the  cutoff  energy  for  the  plane-wave  expansion  is  300  eV. 
We  used  a  supercell  approach  and  monkhorst-pack  A-point 
mesh  for  brillouin  zone  integrations.  Wurtzite  supercells 
containing  36  atoms  were  used  with  the  A'-point  mesh  of 
2x2x2.  Test  calculations  of  up  to  a  96-atom  cell  have 
shown  that  this  is  adequate  for  well-localized  defects  being 
studied  here.  For  defects  with  non-zero  charge;  a  jellium 
background  is  used  to  prevent  the  divergence  of  the 
madelung  sum.  The  definition  of  formation  energy  of  a 
defect  from  first-principles  calculations  can  be  found 
elsewhere  (e.g..  Refs.  [20,21]). 

3.  Results  and  discussions 

3.1.  General  properties  of  N-doped  ZnO 

Nitrogen  is  a  group-V  element;  hence,  it  should  be 
possible  to  use  as  a  p-type  dopant  by  substituting  it  for  O  in 
ZnO.  An  O-poor  growth  ambient  is  generally  required  for 
p-type  doping  by  nitrogen  [22].  In  our  experiments,  if 
oxygen  is  used  as  an  oxidizer  and  NO  gas  is  used  as  a 
dopant,  the  resulting  ZnO  film  contains  no  detectable 
nitrogen  doping  [23].  Using  NO  gas  for  both  oxidizer  and 
dopant,  nitrogen  concentrations  as  high  as 
2.60  x  1021  cm-3  were  achieved.  This  high  concentration 
of  nitrogen  can  be  achieved  in  a  wide  processing  range. 
Nevertheless,  the  growth  window  for  a  successful  fabrica¬ 
tion  of  p-type  film  is  very  narrow.  The  p-type  ZnO:N  films 
are  achieved  only  in  a  growth  temperature  range  of 
400-440  °C  [6], 

A  typical  topography  of  the  MOCVD-formed  ZnO:N 
films  on  a  glass  substrate  is  illustrated  in  Fig.  1.  The  figure 
clearly  shows  that  the  film  is  polycrystalline  with  the  grain 
size  smaller  than  100  nm,  and  an  average  surface  roughness 
of  around  10-20  nm.  XRD  analysis  indicated  that  the 
ZnO:N  films  maintained  the  Wurtzite  structure  of  undoped 
ZnO.  However,  compared  with  the  strong  (002)  orienta¬ 
tion  preference  in  undoped  ZnO  film,  ZnO:N  deposited  on 
glass  lost  the  (002)  orientation  preference  and  became 
randomly  oriented,  as  shown  in  Fig.  2.  The  transmittances 
of  undoped  ZnO  and  ZnO:N  films  are  shown  in  Fig.  3. 
Because  of  the  wide  fundamental  optical  bandgap, 
undoped  ZnO  is  colorless  and  highly  transparent.  The 
increased  absorption  near  the  band  edge  made  the 
nitrogen-doped  ZnO  film  yellowish. 

Results  of  Hall  probe  analysis  on  several  typical  ZnO 
samples  are  shown  in  Table  1.  Undoped  ZnO  is  n-type  with 
a  carrier  concentration  in  the  range  of  1016-1018cm“3.  The 
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Fig.  1.  The  topography  of  a  0.74-pm-thick  ZnO:N  film  grown  on  a  glass 
substrate  at  400  °C. 


Fig.  2.  X-ray  diffraction  data  of  ZnO  and  ZnO:N  films.  Both  films  were 
deposited  on  glass  substrate  at  400  °C. 

nitrogen-doped  ZnO  samples  are  no  longer  n-type,  but  are 
either  semi-insulating  or  weak  p-type  with  the  carrier 
concentration  of  around  1013-10l8cm-3.  A  double-Schott- 
ky  surface-contact  C— V  technique  was  used  to  confirm  the 
Hall  results.  The  Hg  probe  contact  makes  it  possible  to 
contact  both  p-  and  n-type  samples.  In  this  method,  the 
conductivity  type  of  the  samples  is  determined  from  the 
response  of  the  depletion-layer  width  against  the  polarity 
of  the  bias  applied  to  the  inner  main  electrode.  A  ZnO:N 
sample  measured  by  Hall  with  a  p-type  carrier  concentra¬ 
tion  of  6.6  x  1013cm-3  was  analyzed  by  C-V.  As  seen  in 
Fig.  4,  the  C-V  technique  confirmed  that  this  sample  was 
p-type  (positive  slope  under  reverse  bias),  with  calculated 
carrier  concentration  of  ~10l4cm-3. 


Wavelength  (nm) 

Fig.  3.  Optical  transmittance  of  ZnO  and  ZnO:N  films.  Both  films  were 
deposited  on  glass  substrate  at  400  °C. 


The  p-type  behavior  was  achieved  only  with  the  growth 
temperatures  of  about  40CM150  °C.  The  best  electronic 
properties  achieved  in  this  temperature  range  were  on  the 
sapphire  substrate  with  carrier  concentration  of 
8.36  x  10l7cm-3,  Hall  mobility  of  4.55cm2/Vs,  and 
resistance  of  1.64  Q  cm.  At  growth  temperatures  lower 
than  400  °C,  the  him  deposition  rate  was  low  and  the  him 
was  n-type.  When  deposition  temperatures  were  increased 
to  500  °C,  the  him  reverted  to  n-type.  We  will  discuss  the 
temperature  dependence  of  conducting  type  in  the  later 
sections. 

3.2.  Impurity  and  defect  complexes  analysis 

In  this  work,  we  focus  on  the  effects  of  hydrogen  and 
carbon  impurities.  Hydrogen  is  an  impurity  that  can  be 
present  in  ZnO  made  by  many  fabrication  techniques 
[3,9,23-25].  An  isolated  interstitial  H  is  predicted  to  be 
exclusively  donor  in  ZnO,  i.e.,  it  always  exists  in  the  form 
of  an  H+  and  bonds  with  the  host  O  atom  [26].  Carbon 
impurities  are  of  special  concern  with  MOCVD-grown 
hlms  because  of  the  metal-organic  precursors  that  are  used. 
Previous  work  has  demonstrated  that  this  growth  techni¬ 
que  results  in  bulk  carbon  concentrations  much  higher 
than  other  growth  techniques  such  as  reactive  sputtering 
[12].  In  n-type  ZnO,  the  carbon  could  substitute  on  the  O 
site  and  be  an  acceptor.  A  recent  study  indicated  that 
carbon  could  exist  as  a  donor  in  ZnO  in  several  forms  [27]. 
For  example,  in  semi-insulating  or  p-type  ZnO,  carbon  can 
form  diatomic  species  with  oxygen  or  nitrogen  and  share 
an  oxygen  site. 

The  hydrogen  and  carbon  impurity  levels  in  the  ZnO 
hlms  were  studied  by  SIMS  depth-profiling.  The  measure- 
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Table  1 

Electrical  properties  of  ZnO  films  deposited  on  glass  substrate  at  deposition  temperature  of  400  °C 


Doping  condition 

Sample  ID 

CC  (cm-3) 

fi  (cm2/V  s) 

p  (Ocm) 

Conductivity  type 

Undoped 

L3353.2 

-8.42  x  101S 

3.2 

0.235 

n 

L4084.1 

-8.38  x  1016 

6.3 

11.8 

n 

N-doped 

L3357.2 

2.90  x  1014 

579 

37.2 

p 

L4098.2“ 

8.36  x  1017 

4.55 

1.64 

p 

(For  Hall  measurement,  the  error  bars  are  very  small  for  carrier  concentrations  higher  than  10l5cm~3  and  mobilities  greater  than  1  cm2/V s.  In  this  table, 
the  sample  L3367.2  has  a  large  error  bar.  Carrier  concentration  data  was  measured  in  the  region  between  2.9  x  1014  and  4.8  x  1014ctrC3. 

“Sapphire  substrate. 


Fig.  4.  Measured  C-V  characteristics  of  1/C2  vs.  reverse-bias  voltage  on 
nitrogen-doped  ZnO  film.  The  film  is  1.8  pm  thick  and  deposited  on  a  glass 
substrate. 

ment  on  an  undoped  ZnO  film  fabricated  with  MOCVD 
was  used  as  a  benchmark.  The  hydrogen  and  carbon 
concentration  as  a  function  of  growth  temperature  for 
undoped  ZnO  film  is  shown  in  Fig.  5.  The  hydrogen  and 
carbon  concentrations  are  high  at  low  deposition  tempera¬ 
ture.  As  the  deposition  temperature  increased,  both 
hydrogen  and  carbon  concentrations  decreased  and  likely 
stabilized  around  400  °C.  The  similar  dependence  of  carbon 
and  hydrogen  concentrations  on  deposition  temperature 
indicates  that  carbon  and  hydrogen  possibly  exist  together 
as  defect  complexes. 

For  nitrogen-doped  ZnO  films,  the  dependencies  of  the 
hydrogen  and  carbon  concentration  on  the  growth 
temperature  are  different  from  the  case  of  unintentionally 
doped  ZnO  films.  Fig.  6  shows  that,  in  the  ZnO:N  samples, 
the  hydrogen  and  carbon  concentrations  did  not  niono- 
tonically  decrease  as  deposition  temperature  increased.  The 
samples  grown  at  400  °C  had  the  lowest  hydrogen  and 
carbon  concentrations.  As  the  deposition  temperature 


100  200  300  400  500  600 

Deposition  Temperature  (°C) 


Fig.  5.  A  SIMS  analysis  of  hydrogen  and  carbon  concentration  on  the 
MOCVD-fabricated  ZnO  samples.  The  ZnO  films  were  deposited  at 
400  °C  with  DEZ  and  O,  flow  rate  of  1  and  100  seem,  respectively. 


Fig.  6.  A  SIMS  analysis  of  hydrogen,  carbon,  and  nitrogen  concentration 
on  MOCVD-fabricated  ZnO:N  samples.  The  ZnO:N  films  were  deposited 
at  400  °C  with  DEZ  and  NO  gas  flow  rate  of  1  and  36  seem,  respectively. 
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increased  to  above  400  °C,  hydrogen  and  carbon  concen¬ 
trations  both  increased  again.  If  we  compare  the  hydrogen 
and  carbon  concentrations  between  the  nitrogen-doped 
and  undoped  samples,  we  find  that  the  hydrogen  concen¬ 
trations  are  about  the  same  level.  However,  the  carbon 
concentrations  in  nitrogen-doped  samples  are  about  ten 
times  higher  than  those  in  the  undoped  samples.  A  possible 
explanation  for  the  increased  carbon  impurities  in  the 
ZnO:N  samples  will  be  provided  later  in  this  paper.  In  any 
case,  it  is  clear  that  under  current  film  growth  conditions, 
when  NO  gas  used  as  an  oxidizer  instead  of  02  gas,  the 
carbon  concentrations  in  the  ZnO  film  increase. 

The  FTIR  results  strongly  support  the  hypothesis  that  at 
least  some  of  the  carbon  and  hydrogen  co-exist  as  defect 
complexes  in  the  undoped  ZnO.  Fig.  7  shows  the  FTIR 
absorbance  spectra  in  the  2750-3050  cm-1  region.  Two 
absorbance  spectra  were  taken  from  the  undoped  and 
nitrogen-doped  ZnO  samples  that  were  both  deposited  at 
400  °C.  On  the  undoped  ZnO  film  (dotted  line),  we 
observed  a  series  of  sharp  absorption  peaks  in  the  range 
consistent  with  CH  stretch-mode  frequencies.  These 
specific  peaks  correlate  well  with  the  observed  frequencies 
of  the  CHA  (x  =  1,  2,  and  3)  stretch  modes  [25,28]. 
Combined  this  with  our  impurity  analyses  results,  we 
believe  that  in  our  MOCVD-fabricated  ZnO  films,  a  large 
fraction  of  carbon  is  present  in  the  form  of  CHV  defect 
complexes.  These  complexes  most  likely  come  from  metal- 
organic  precursors. 


FTIR  absorbance  spectra  of  the  ZnO:N  sample  further 
indicated  that  the  NO  gas  affects  the  formation  of  the  CHV 
complexes  in  the  ZnO  film.  Although  SIMS  data  show  the 
increase  in  carbon  concentration  for  the  ZnO:N  sample, 
Fig.  7  shows  that  with  nitrogen  doping  (solid  line),  the 
sharp  IR  absorbance  peaks  associated  with  the  CHX 
complex  are  largely  decreased.  The  new  broad  peak 
emerges  [19]  around  3020  cm”1  (assigned  as  N0-H  peak 
[29]).  This  observation  indicates  that  the  use  of  NO  gas  as  a 
precursor  reduces  the  absorption  of  CHV  radicals  on  the 
surface  of  the  substrate  and/or  that  the  existence  of 
nitrogen  in  the  ZnO  film  strongly  attracts  the  impurity  H 
and  promotes  the  break-up  of  CHA.  defect  complexes.  The 
other  possibility  is  that  with  nitrogen  doping,  a  reaction 
between  nitrogen  and  carbon  occurs,  which  we  discuss 
next. 

There  is  some  correlation  between  carbon  and  nitrogen 
in  the  nitrogen-doped  ZnO  film.  SIMS  results  indicated 
that  the  carbon  concentrations  in  the  ZnO:N  films  are 
generally  higher  than  those  in  undoped  ZnO  films.  Fig.  8 
shows  that  for  the  nitrogen-doped  ZnO  films,  very  strong 
IR  absorption  peaks  exist  in  the  range  of  frequencies 
consistent  with  NC,  NN,  and  CO  modes  [19],  Based  on 
first-principles  calculations,  the  formation  energies  of  the 
NC  and  NN  complexes  are  much  lower  than  for  CO  [27]. 
Therefore,  the  peaks  around  the  wave  number  of 
1840  cm”1  are  quite  possible  due  to  the  NC  and  NN 
complexes.  Our  previous  work  suggests  that  the  NC 


Wavenumber  (cm'1) 


Fig.  7.  FTIR  absorbance  spectra,  showing  the  2750-3050cm_1  region, 
were  taken  from  undoped  (dotted  line)  and  nitrogen-doped  (solid  line) 
ZnO  films.  Both  ZnO  samples  were  deposited  at  400  °C.  The  absorption 
peaks,  which  correlate  well  with  the  frequencies  of  the  CH*  (x  =  1.  2,  and 
3)  stretch  modes  [25,28],  are  clearly  shown  on  the  undoped  ZnO  film,  but 
not  on  the  nitrogen-doped  film. 


Fig.  8.  FTIR  absorbance  spectra,  showing  the  1400-2200 cm-1  region, 
were  taken  from  the  same  samples  as  Fig.  7.  Compared  with  the  undoped 
(dotted  line)  ZnO  film,  the  nitrogen-doped  (solid  line)  ZnO  film  has  very 
strong  absorption  peaks  in  the  frequency  region  possibly  due  to  NC,  NN, 
and  CO  stretch  modes. 
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complex  is  far  more  likely  in  the  MOCVD-formed  ZnO:N 
film,  compared  with  the  NN  defect  complex  [12].  There¬ 
fore,  in  the  ZnO:N  film  fabricated  with  the  DEZ  and  NO 
gas,  the  majority  of  the  IR  absorption  peaks  around  the 
wave  number  of  1840  cm"1  may  be  due  to  the  NC  complex. 

In  order  to  produce  high  quality  p-type  ZnO:N,  carbon 
contamination  must  be  reduced.  However,  identifying  and 
reducing  the  carbon  incorporation  pathways  for  MOCVD 
grown  ZnO:N  is  quite  difficult.  The  decomposition  kinetics  of 
DEZ  under  vacuum  or  inert  atmospheres  are  understood 
fairly  well  [30,31],  but  the  presence  of  02  or  NO  in  the  growth 
chamber  will  significantly  effect  the  decomposition  chemistry. 
Ethyl  (C2H5)  radicals  and  ethylene  (C2H4)  are  expected  to  be 
highly  reactive  with  O,  N,  and  NO,  leading  to  many 
pathways  for  carbon  incorporation  [32-36].  These  reaction 
products  may  then  incorporate  in  the  film  as  the  C’HX  and 
NC  species  evident  in  Figs.  6-8.  The  high  reactivity  of  ethyl 
and  ethylene  with  N  may  also  account  for  the  differences  in 
carbon  incorporation  between  the  doped  and  undoped  films. 
Further  work  on  the  kinetics  of  ZnO:N  growth  and  carbon 
incorporation  is  ongoing. 

3.3.  The  possible  compensation / passivation  effects  by 
hydrogen  and  carbon 

Hydrogen  passivation  appears  to  be  a  very  important 
issue  for  p-type  ZnO  doping.  Isolated  interstitial  hydrogen 
was  predicted  to  be  an  exclusive  donor  (H  +  )  in  ZnO  [37]. 
The  interstitial  H2  molecule  was  predicted  to  be  unstable  in 
ZnO,  except  in  an  extreme  n-type  (Fermi  level  (E¥)  is  very 
close  to  the  conduction-band  minimum).  The  calculated 
formation  energy  (per  H  atom)  of  the  H2  molecule  in  ZnO 
is  shown  in  Fig.  9,  in  comparison  with  the  formation 
energy  of  an  isolated  H+  and  N0-H  complex  in  ZnO. 
Under  a  Zn-rich  conditions,  the  formation  energy  of  the  H2 
molecule  in  ZnO  has  a  higher  formation  energy  than 
isolated  H+  and  N0-H  complexes  in  almost  the  entire 
range  of  the  Fermi  energy. 

The  H2  complexes,  which  have  been  predicted  to  be 
stable  in  GaP:N  and  metastable  in  GaN,  have  high 
formation  energy  in  ZnO  and  are  thus  unlikely  to  exist 
[38].  That  H2  complexes  are  predicted  to  be  unstable  in 
ZnO  is  mainly  because  H-  is  not  stable  in  ZnO  [39].  A  H2 
complex  can  be  viewed  as  a  complex  between  H+  (bonded 
with  O)  and  H-  (bonded  with  Zn). 

For  nitrogen-doped  ZnO,  one  of  the  hydrogen  passiva¬ 
tion  effects  has  been  discussed  in  detail  in  a  previous  paper 
[29].  We  also  investigated  the  possibility  of  forming  the 
N0  H2  complex.  In  Fig.  9,  the  calculated  formation  energy 
of  the  Nq-H2  complex  in  ZnO  is  shown  (solid  line)  in 
comparison  with  the  sum  of  the  formation  energy  of  an 
isolated  H+  and  N0-H  complex  (dotted  line)  in  ZnO. 
Because  the  N0-H  complex  has  a  neutral  charge  state, 
there  is  no  Coulomb  driving  force  for  the  complex  to  bind 
with  another  H+  to  form  a  NQ-H2  complex.  Indeed,  it  is 
more  difficult  to  form  a  NQ-H2  complex  than  an  isolated 
H2  (which  is  also  not  stable)  in  ZnO.  Therefore,  the  main 


Fig.  9.  The  formation  of  the  No-Fh*  complex.  Dashed  lines  represent  the 
formation  energies  of  N0-H  complex  and  an  isolated  H+  defect.  The  solid 
lines  represent  the  formation  energy  of  the  No-H2*  complex.  The  dotted 
lines  show  the  sum  of  the  formation  energies  of  the  No-H  complex  and  an 
isolated  H+  defect.  The  N0-H2*  complex  has  a  total  formation  energy 
higher  than  the  sum  of  individual  defects,  which  means  the  complex  is 
unstable.  The  Zn-rich  condition,  N2  phase  precipitation  limits  (for  N 
chemical  potential),  and  H2  phase  precipitation  limits  (for  H  chemical 
potential)  were  used  in  the  calculations.  The  formation  energy  of  the  H2 
molecule  in  ZnO  is  also  shown.  The  Ep  is  relative  to  the  valence-band 
maximum. 

effect  of  hydrogen  impurities  in  nitrogen-doped  ZnO 
should  be  to  form  N0-H  neutral  complexes,  and  thus  to 
passivate  the  N0  acceptors. 

Front  the  above  discussion,  we  can  see  that  for  an 
undoped  ZnO,  carbon  impurities  most  likely  exist  as  CHV 
defect  complexes.  In  the  case  of  nitrogen  doping,  there  are 
fewer  CHV  complexes  and  other  forms  of  carbon  impurities 
such  as  the  NC  complexes  are  present.  Based  on  first- 
principles  calculations,  we  investigated  several  carbon- 
related  defects  in  ZnO.  Interstitial  carbon  (Q)  was  found  to 
be  a  donor,  whereas  substitution  carbon  (C0)  is  a  neutral 
or  deep  acceptor.  Both  C;  and  C0  have  high  formation 
energies;  therefore,  they  should  not  exist  in  high  concen¬ 
tration.  The  N0  and  C;  strongly  bound  to  form  a  (NC)o 
complex.  The  strong  binding  energy  suggests  that  under 
equilibrium  Zn-rich  conditions,  (NC)0  defect  complexes 
are  likely  to  form.  Among  carbon-related  defects,  the 
(NC)o  complex  has  a  lower  formation  energy  than  Q  and 
C0-  The  predictions  well  support  the  experimental  ob¬ 
servation  that  the  carbon  solubility  increases  in  the 
nitrogen-doped  samples  in  comparison  to  the  undoped 
samples. 

First-principles  calculations  predict  that  at  almost  any 
Fermi  level  (except  very  near  the  conduction-band  mini¬ 
mum),  the  (NC)0  is  a  donor,  stable  in  the  1+  charge  state, 
and  has  a  vibrational  frequency  of  1995  cm-1  [27].  There¬ 
fore,  the  increased  carbon  concentration  in  the  ZnO:N 
films  will  compensate  the  N0  doping  effect. 
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4.  Summary 

We  have  discussed  hydrogen  and  carbon  impurities  and 
their  potential  compensation/passivation  effects  in 
MOCVD-fabricated  and  nitric  oxide-doped  ZnO:N  films. 
Compositional,  optical,  electronic,  and  element  chemical 
states  characterization  were  performed  comparatively  on 
both  undoped  and  nitrogen-doped  samples. 

By  MOCVD  growth  using  NO  doping,  we  can  achieve 
very  high  nitrogen-doping  concentration  in  ZnO  films. 
However,  hydrogen  and  carbon  are  unintentionally  doped 
into  our  MOCVD-fabricated  ZnO:N  films.  The  carbon 
concentrations  increase  greatly  in  nitrogen-doped  samples, 
suggesting  that  the  unintentionally  doped  carbon  impu¬ 
rities  prefer  to  form  defect  complexes  with  nitrogen. 
Because  of  the  strong  bonding  between  hydrogen  and 
carbon  with  nitrogen,  the  (N0-H)0  and  (NC)o+  defect 
complexes  have  quite  low  formation  energies  and  could 
exist  in  large  concentrations.  Our  result  illustrates  that  with 
MOCVD  growth  of  N-doped  ZnO,  strict  control  of  carbon 
and  hydrogen  concentrations  may  be  required  to  achieve 
high  carrier  concentrations  in  p-type  ZnO:N  films. 
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